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Wear and Friction Characteristics of 
Structural Materials in Liquid Sodium 


By W. A. Glaeser* 


ABSTRACT: A review of current friction and wear technology 
in sodium-cooled reactors indicates that excessive wear and 
high friction constitute design problems. Selection of materials 
and determination of safe load levels on the basis of friction- 
or wear-failure criteria is based on data from tests in which 
contact conditions are closely simulated. Current wear theory 
must be applied both in establishing test parameters and in 


applying test data to design. The surface chemistry of 


molybdenum-, chromium-, and tungsten-containing alloys is an 
important factor in sliding contact in sodium environment. 


CRITICAL WEAR SITES 


Fast-flux reactors have inherent wear and friction 
problems associated with the use of liquid-sodium 
coolant. Mechanical failure of devices involving con- 
tacting sliding or stationary surfaces immersed in 
sodium can occur by adhesion or by severe wear. Some 
configurations that have been recognized as potentially 
critical friction and wear sites include fuel-pin space 
separators, fuel-assembly duct pads, control-rod-drive 
roller nuts, core-support-structure positioning surfaces, 
and fuel-handling devices (grappling-, linear-, and 
rotary-motion mechanisms). These mechanisms see a 
variety of contact conditions, including contact pres- 
sures ranging from 100 psi to several thousand psi, 
stationary contact, sliding contact, small amplitude or 
fretting motion, and rolling contact. Their operating 
environment is high-purity sodium (under 5 ppM 
oxygen) at temperatures from 400 to 1350°F. 

To illustrate the type of friction and wear problems 
being encountered in the above mechanisms, | will 
briefly describe two components currently under in- 





*Lubrication Mechanics Division, Battelle Memorial In- 
stitute, Columbus Laboratories, Columbus, Ohio 43201. 


tensive investigation for the Fast Flux Test Facility 
(FFTF) program: (1) fuel-assembly duct pads and 
(2) fuel-pin wire-wrap spacers. 


Fuel-Assembly Duct Pads 


The general design of the FFTF core involves 
bundles of fuel pins contained in stainless-steel vertical 
hexagonal pipes or ducts. The ducts are clamped in a 
close-packed array when assembled in the core. When 
the ducts are in position in the core, liquid sodium 
flows through the ducts extracting heat from the fuel 
pins and also surrounds each duct. To allow room for 
circulation of the sodium, each duct is separated from 
its neighbor by spacers. The spacers are pads integral 
with the ducts and constitute several circumferential 
bulges at even axial spacing on each duct. The ducts 
contact each other via their duct pads or spacers. 
During core assembly or after removal and reinsertion 
of a fuel-assembly duct, the ducts must be compacted 
into a close-packed array by means of a hold-down 
mechanism. The compacting is done in 400°F sodium. 
If excessive friction is experienced between contacting 
duct pads in this environment, the required close 
seating of all ducts in the array will be prevented. 
Clamping forces available limit the maximum tolerable 
friction level. Experiments are currently being carried 
out at several laboratories (Liquid Metal Engineering 
Center and Westinghouse Advanced Reactors Division) 
to develop minimum friction surfaces for duct-pad use 
in high-temperature sodium. 

When the reactor is in operation, duct pads will be 
in contact under a bearing pressure—from hold-down 
forces—of several thousand psi at a temperature of 
about 1150°F. The duct pads can remain under these 
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conditions for as long as a month. Adhesion of the 
contacting surfaces after service in the reactor could 
prevent removal of a duct. Therefore selection of 
materials requires that tendency for self-welding in 
liquid sodium be known. 

The duct pads rub up and down against each other 
during reactor operation owing to unequal expansion 
of core components and vibrations caused by flowing 
sodium. Wear of the duct pads must be limited to 
maintain design spacing and to keep wear debris 
circulating in the liquid sodium to a minimum. It is 
necessary to know the load capacity (beyond which 
severe wear or galling occurs) of spacer materials to 
ensure against this type of failure. 


Fuel-Pin Wire Wrap 


Individual fuel pins in a bundle are to be kept 
separated from each other and the duct surface by a 
wire-wrap configuration. This design involves use of a 
stainless-steel wire wrapped around each pin in a 
long-pitch helix. The resulting contact conditions are 
shown in Fig. 1. Wire wrap—duct contact is approxi- 
mately a cylinder-flat configuration and the wire 
wrap- pin contact is approximately a crossed-cylinder 
contact condition. In either case, concentrated contact 
conditions exist. During reactor operation, these wire- 
wrap contacts will be subjected to small amplitude 
sliding motion as a result of uneven expansion and 
contraction of individual pins. Vibration of the pins 
induced by sodium flow will result in a combined 


CLADDING 







=F VIEL 


~CONTACT POINT 


— FUEL PIN 


——WIRE WRAP 


Fig. 1 Contact conditions in wire-wrap separators for fuel 
pins. 
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impact and scrubbing motion—a particularly severe 
condition for fretting. This presents a difficult high- 
temperature wear problem (about 1350°F) at the wire 
wrap—fuel pin contacts. Penetration of the fuel-pin 
cladding by the wire wrap during fretting and impact 
must be allowed for in selecting cladding thickness. 
This is particularly difficult because wire wrap—pin 
contact stresses are not established (fuel-pin bowing 
can produce unknown variations in pin-to-pin loading). 
Therefore penetration rates cannot be predicted with 
sufficient precision. Efficient heat transfer requires 
minimizing cladding thickness and not much of the 
thickness allowance can be set aside for wear. Ultimate 
concern is to be able to reliably predict how much 
cladding must be allowed to prevent penetration to the 
fuel. 

Aside from conventional friction and wear con- 
siderations in design, sliding contact in high- 
temperature, high-purity sodium presents special 
problems of its own which cannot be solved by 
application of classic wear principles. For one thing, 
compatibility with high-temperature sodium limits the 
selection of contact materials to materials like 316 
stainless steel, a material well known for its tendency 
to gall and cold weld. 

In reviewing the current state of liquid-sodium 
wear and friction technology, one finds that progress is 
being made in developing approaches to designing for 
wear in liquid sodium. At the same time, there are still 
many pitfalls and limitations to this field. This paper is 
intended to point out progress and pitfalls. 

As a basis for analyzing wear and friction problems 
in sodium, one starts with current wear theory. The 
modifications and special considerations that sodium 
environment imposes can then be added. 


WEAR THEORY 


Over the past 20 years or so, a general wear 
equation has been developed based largely on empirical 
observations: 


y=—— (1 
Pm ) 
or 
KwL 
d= (2) 
Pn 


where V = wear volume 
N = load 
L = sliding distance 





WEAR AND FRICTION CHARACTERISTICS OF 


Pm = pressure at which plastic penetration occurs, 
i.e., hardness (Brinell, Diamond, Pyramid, 
Vickers) of the softer member of the sliding 
couple 

d= apparent depth of wear (assuming even wear 
over contact area) 

w=apparent bearing pressure = load/apparent 
contact area 

K = wear coefficient 


This equation, in one form or another, has been 
developed from different approaches by Archard,'*? 
Kragelskii,*** and Rabinowicz.° 

Before the limitations of this equation are dis- 
cussed, it would be appropriate to define some of the 
essential parameters involved in the wear process since 
there is often much confusion about what is meant by 
“wear.” 

There are several modes of wear which can occur 
and must be recognized as having differing effects on 
the performance of mechanical devices. Mild abrasive 
or adhesive wear generally occurs under lubricated 
conditions, and it is the condition under which the 
general wear equation is valid—that is, wear is a linear 
function of sliding distance and proportional to load.* 
In this form of wear, rubbing surfaces contact through 
a few high points or asperities (these could be 
embedded or trapped hard particles) and either shear 
off as particles or scratch the mating surface. Asperity 
shearoff is adhesive wear, and asperity scratching or 
ploughing is abrasive wear. In mild wear, the difference 
between adhesive and abrasive is not significant as far 
as the general equation is concerned, and the final 
effect is a gradual even polishing of the surface by 
removal of extremely fine particles (in the micron 
range). 

Severe wear is a second mode which occurs most 
often in unlubricated or poorly lubricated sliding 
contact or when lubrication breaks down because of 
“overload.” Severe wear is characterized by tearing or 
deep scratching of the surface and often accompanies 
formation of growing adhesive junctions or coalescence 
of several growing asperity junctions. This mode of 
wear does not follow the general wear equation, and 
the end result depends on whether abrasive or adhesive 


wear is Occurring. 





*Thus wear can be reliably predicted on the basis of 
appropriate wear tests. 

+For instance, in severe abrasive wear, the surface damage 
can be eliminated by making the surface harder than the 


4 


abrasive medium*"* for adhesive wear; on the other hand, this 


principle does not necessarily apply. 
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Fretting wear occurs under very low amplitude 
oscillation between contacting surfaces and results in 
localized removal of surface material (pitting). It is the 
result of continuous plastic deformation of localized 
contact zones and removal and oxidation or corrosion 
of minute particles, presumably by a fatigue process. 
Fretting can occur under corrosive conditions (fretting 
corrosion) or in the absence of a corroding agent 
(fretting wear). Fretting reduces the fatigue strength of 
a structural member if the maximum outer fiber 
bending stress occurs in the region of fretting contact. 
It will cause localized wear or pitting—and can cause 
premature penetration of fuel-pin cladding. Fretting 
does not follow the general wear equation but is 
sensitive to material combination, vibration amplitude, 
and load-displacement mode. 

There are several means of expressing wear 
quantitatively. 


Wear coefficient: 


fi dP» 


m 


VP, e 
K = 1422 i ie... r4z2—— ©) 


— 5 <a 

where V = wear volume, in. 
F ) 
P, = hardness, kg/mm* 





N = load, lb 
L = sliding distance, in. 
d = depth of wear, in. 
P= bearing pressure, psi: M/apparent area of 
contact 
Specific wear: 
V in? 1 in. 
ate” eereere Pe ty ae (4) 
NL (lb) (in.) P(L) (psi) (in.) 
Wear rate: 
Vins +3 dmil 7 
—— or 10°—— 5) 
ae L in. ( 


Note that these expressions give wear as either 
weight or volume loss or penetration. It is important to 
make this distinction when considering wear in engi- 
neering design. The designer first of all must consider 
what the failure mode will be for a given mechanism. 
For instance, for the wire-wrap fuel-pin spacer, wear- 
through or notching of the cladding would constitute 
wear failure. In this case, the engineer is most 
interested in wear in terms of penetration. In general, 
penetration, or loss of critical dimension by wear is the 
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most prevalent failure mode. From Eq. 1, it is evident 
that the general wear equation is based on volume wear 
or weight loss. Equation 2 is a modification of Eq. 1 to 
provide wear in terms of penetration; however, it is 
true only if wear or deep scratches over the apparent 
contact area are evenly distributed (not the case in 
fretting scoring, for instance). Thus, in selected design 
applications, wear values in terms of linear dimensional 
loss (opening of clearance and penetration of surface) 
can be obtained by dividing the wear coefficient, K, by 
the apparent area of contact. 

The linearity of the general wear equation holds 
only for contact geometries where volume loss is 
proportional to linear loss—as in a flat-on-flat con- 
figuration. However, for line contact or point contact, 
linear wear is not proportional to volume wear. A plot 
of depth of penetration vs. volume wear is shown in 
Fig.2 for cylinder vs. flat and ball vs. flat contact 
conditions (popular specimen configurations for wear 
testing). For the line contact and point contact, 
penetration is roughly proportional to volume wear or 
weight loss for the first few percent penetration (/). 
Beyond that, however, wear coefficients based on 
volume wear rates would deviate more and more from 
true penetration (A) and thus predict much larger 
values than, in fact, occur. The parameter w, or width 
of wear scar, commonly used as a measure of wear in 
wear testing (Four Ball and Timken), is a poor 
indicator of volume wear, especially at small scar 
widths. Thus wear coefficients or specific wear values 
cannot be used interchangeably with those based on 
depth of penetration or width of wear scar where 
point- or line-contact conditions exist. If an existing K 
value is to be used, it must first be expressed in volume 
wear with the general wear equation and converted to 
depth of penetration with a chart such as that in Fig. 2. 

The general wear equation indicates that wear is 
proportional to total sliding distance. If a wearing 
system complies with this characteristic, then total 
wear can be predicted or extrapolated from relatively 
short-term wear tests designed to establish the slope of 
the curve or wear rate. However, when this approach is 
being applied to wear problems, several things must be 
kept in mind: (1)as discussed previously, the wear 
equation is based on volume wear rate, and contact 
geometry must be considered when wear data from a 
given configuration are being applied to a design having 
a different configuration; (2) environment has a power- 
ful influence on wear, and therefore data must be 
obtained under the same surface chemistry conditions 
as anticipated in the application; and (3) wear is rapid 
and nonlinear in the early or wear-in stage, and a wear 
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Fig. 2 Relation between depth of penetration and volume 
wear. 


test must be run until linear wear is established before 
determination of the wear rate, or, if only one wear 
measurement can be made, the test should be suf- 
ficiently long to make the effect of wear in the 
resulting wear-rate determination negligible. This con- 
cept is illustrated in Fig. 3. 





TOTAL WEAR 














SLIDING DISTANCE ———> 


Fig. 3. Effect of rapid initial wear on wear characteristics. 


Wear rate will change drastically when load is 
increased beyond a critical level. Experience shows that 
most sliding systems exhibit a critical load above which 
wear becomes extremely high and surface damage 
occurs (severe wear regime). The critical load character- 
istics is shown schematically in Fig.4. The bearing 
stress (P.) at which gross yielding occurs usually 
coincides with a rapidly rising wear coefficient. The 
critical bearing stress (P.) can be estimated from 
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Fig. 4 Relation between contact pressure and wear coeffi- 
cient. 


indentation hardness values (Brinell, Vickers, and 
Knoop). 

Welsh® has shown that for medium-carbon-steel 
(0.10 to 0.90 wt.% carbon) unlubricated sliding in air 
environment, two abrupt transitions in wear rates can 
occur as the load is increased. These transitions are 
shown in Fig. 5. This shows that with this material, 
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Fig.5 Wear mode transitions for medium-carbon steel as 
influenced by applied load. 
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after the transition to severe wear, if the load is 
increased sufficiently, a sudden decrease in wear rate 
occurs. The shape of the curve in Fig.5 reflects 
particular material properties and surface heating ef- 
fects. In mild steel, mild wear prevails as long as a 
protective oxide is supported on the surface. Increasing 
the load breaks up the oxide film since surface heating 
weakens support of the film. Heavy loads produce 
sufficient plastic flow to transform the surface to a 
hard ‘friction phase’—-probably martensite—which 
has sufficient penetration resistance to support a 
protective oxide film. This characteristic may not hold 
for materials other than mild steel (Welsh found no 
transition for high carbon steel). The surface chemistry 
with sodium as the environment imposes some special 
characteristics for such a sliding system and will be 
discussed later in more detail. 

When wear data are being used for design of 
mechanical components, the two concepts described 
above might be considered: (1) estimation of total 
wear over a required life and (2) maximum load 
capacity without exceeding the severe wear limit. Both 
of these approaches require wear test data, and the 
method of testing should be based on the principles 
discussed previously. For instance, because of the 
effect of contact geometry, material hardness (in- 
fluenced by temperature), and environment, testing 
should simulate actual contact conditions as closely as 
possible. The term P,, in Eqs. 1 and 2 varies with 
temperature. Surface chemical effects are influenced 
by temperature and not covered in Eqs.1 and 2. If 
load capacity or the critical load at which wear rate 
changes drastically is to be determined, the lower-load 
range should be explored, as well as the upper range, to 
establish whether more than one critical load level 
exists. Reduction of design load on the basis of 
heavy-load tests might not be sufficient because one 
could get into an even higher wear rate at an 
intermediate load. Austenitic stainless steel will trans- 
form at the surface under heavy-sliding contact stress. 
The resulting “friction phase’—presumably — un- 
tempered martensite—acts in the same way Welsh 
described for mild steel. 


FRICTION THEORY 


Modern friction theory is based on the concept of 
true area of contact. That is, when surfaces contact 
under loads that do not produce gross yielding, actual 
surface contact occurs only at a few minute high points 
or asperities. The friction force, then, is the force 
required to shear the asperity junctions (shearing 


REACTOR TECHNOLOGY, Vol. 15, No. 1, Spring 1972 





6 WEAR AND FRICTION CHARACTERISTICS OF MATERIALS IN SODIUM 


produces wear debris). The force will be controlled by 
the shear strength of the adhesive bond and the yield 
strength of the weaker of the two contacting materials. 
Bowden and Tabor’ have proposed that the size of the 
real contact area is proportional to the load and 
inversely proportional to the penetration hardness of 
the softer material in the contact pair; that is, 


W 
A “| (6) 


where A = area of real contact, in.” 
W = load, lb 
P = indentation hardness, psi 


Shearing of these contact areas produces the friction 
force, 


(7) 


where F = friction force, lb 
= coefficient of friction 
S=shear strength or cleavage strength of the 
bond or weaker material in the contact, psi 


The above relation shows that the friction coefficient is 
independent of load and apparent area of contact. It is 
also independent of sliding velocity so long as frictional 
heating does not change the mechanical properties of 
the surface. Since the shear strength of many materials 
does not exceed about one-half of their indentation 
hardness, friction coefficients above 0.5 should not be 
expected. In actual fact, friction coefficients of 1.0 and 
higher are commonly reported for dry, clean, metallic 
surfaces. The simple concept, however, does not take 
into account the plastic flow that occurs during sliding 
of the contact surfaces, resulting in growth of the real 
contact area. Ploughing of hard constituents in one 
surface through softer phases in a mating surface will 
also contribute to the total friction force. 

The shear strength of the contact junctions, S, in 
Eq. 7 is a complex factor and very sensitive to surface 
conditions. In high vacuum, with truly clean surface 
conditions, S approaches the shear strength of the 
softer material in the sliding pair. When the surface is 
“contaminated” with oxides and/or adsorbed gases or 
vapors, the friction coefficient can drop several orders 
of magnitude. Just a few molecular layers of an organic 
material on a metal surface can account for such 
changes in friction behavior. The use of published 
friction coefficients has always been unreliable because 
of unknown surface conditions present during a test. 
This has been especially true in liquid-sodium tech- 
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nology, where wide differences exist in published data 
on friction coefficients for similar materials. 

Even the structure of an adsorbed monolayer can 
influence friction. This was demonstrated by Zisman® 
in a series of elegant experiments in which he caused a 
ball to slide over a glass plate with adsorbed mono- 
layers of various pure hydrocarbons laid down under 
conditions of extreme cleanliness. He showed that 
friction decreased measurably when straight-chain 
molecules with increasingly longer chain lengths were 
used. A minimum friction level was attained when the 
chain length of the molecule adsorbed was that at 
which close packing occurs. Oxides differ in “lubricat- 
ing”’ ability. The lower oxides of iron (Fe304 or FeO) 
which exist in partial vacuum conditions produce lower 
friction than does the Fe,03 found in normal atmo- 
sphere. 

It is possible to enter the high-friction regime of 
clean surfaces when contact conditions (heavy load— 
high speed, for instance) produce frictional heating, 
softening of the surface layers, and breakthrough of 
protective films. This principle is used in the friction 
welding process for joining metals of different 
chemistry. 


Sodium Environment 


The principles of wear and friction discussed 
previously are modified to some extent in sodium 
environment. The wear coefficient, K, in Eq. 1 will be 
modified. Resistance to penetration, P,, will also be 
reduced with increased temperature. Essentially the 
differences are of the same type that are noted when 
any change is made in the environment surrounding 
sliding surfaces. The high-purity sodium environment 
expected in a reactor should provide a more constant 
environmental effect than a normal air atmosphere 
does. However, as will be shown later, in a high- 
temperature sodium environment, surface effects are 
more sensitive to slight shifts in temperature and 
oxygen level. 

Surface reactions are known to take place between 
sodium and materials immersed in the sodium. In 
general, the sodium acts as a reducing agent to clean 
off most films generated during manufacture. As has 
been pointed out, these complex surface contan inants 
play an important role in suppressing adhesicn and 
reducing friction and wear. However, when certain 
metallic elements (notably chromium, tungsten, and 
molybdenum) and a source of oxygen are present in 
sodium, an oxidation process can occur in which a 
sodium-metal-oxide film is formed. This reaction prod- 
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uct has been found capable of influencing wear and re- 
ducing friction in elevated-temperature sodium. The 
process has been likened to boundary lubrication by 
conventional oil lubricants where the oil supplies a 
source of continuous replenishment of the solid 
boundary film as it is worn away. Sodium molybdate, 
formed in the reaction between sodium oxide and 
molybdenum, is a known effective solid lubricant. 

Kissel, Glaeser, and Allen? and, later, Roberts'° 
found that friction of alloys containing molybdenum, 
tungsten, or chromium can be lower in oxygen- 
containing sodium than in an inert-gaseous environ- 
ment. Although the reactions producing beneficial 
films on surfaces immersed in sodium are complex and 
not completely understood, two key reactions can be 
considered as determining whether or not the surface 
chemical effects are probable: 


M+0O, —> M,.0, + AF (8) 
and 


4Na + O, > 2Na, 0 — 135 kcal per mole (9) 
If the free energy (AF) of Eq. 8 is positive or a smaller 
negative number than that of Eq. 9, the sodium oxide 
is more stable, and the probability is that the metal 
oxide will be reduced by the sodium. If the metal- 
oxide reaction shows a larger negative AF, then the 
chances for stable reactions between oxygen, sodium, 
and the metal surface are favorable. 

Roberts'° suggested that below 700 to 800° F, the 
molybdenum, tungsten, and chromium oxides are more 
stable and hence the boundary lubrication mechanism 
is probable for alloys containing these elements. 
However, above 800°F the only possible candidate is 
chromium oxide. This has been demonstrated in 
friction experiments conducted in oxygen-containing 
sodium with two nickel-based alloys having similar 
hot-hardness properties, one containing chromium and 
the other having no chromium. The chromium- 
containing alloys’! performed with lower friction and 
less surface damage than the nonchromium alloy at 
temperatures to 1200°F. 

The reaction process that produces beneficial 
boundary lubrication is apparently sensitive to the 
combined effect of oxygen concentration and sodium 
temperature. Recent studies at the Liquid Metal 
Engineering Center (LMEC)'?"'? using austenitic 
stainless steel as the wear components have shown that 
chromium additions can be beneficial to various 
limiting maximum temperatures, depending on the 
oxygen content of the sodium. This work assumes that 
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the Na—Cr—O compound significantly influences fric- 
tion, and wear data* have been established relating the 
free energy of NaxCryO, to the free energy of oxygen 
dissolved in sodium. The latter assumes, as stated in the 
LMEC work, “that the oxygen activity of the sodium 
with an oxygen concentration of c is approximately 
equal to c/c.,¢ and that the experimental curve 
selected for oxygen saturation in sodium (C.at) as a 
function of temperature is correct.” 

The chart shown in Fig. 6 was developed by LMEC 
from a cross-plot of free-energy curves. It shows two 





Zone of no protection 






Stable boundary film 





TEMPERATURE, °F 











OXYGEN, ppM 


Fig. 6 Combined effect of sodium oxygen level and tempera- 
ture on wear-resistant film formation for 304 stainless 
steel.!)>1? 


regions—one where sodium getters oxygen from the 
stainless steel, and hence no beneficial film forms, and 
the other where a stable boundary film forms and 
provides protection. The chart indicates clearly the 
importance of having enough oxygen present, partic- 
ularly at the higher temperatures, to form a stable, 
protective film. More recent data from LMEC'* 
indicate that, at 5 ppM oxygen, beneficial films may 
form on chromium-containing alloys somewhere 
between 1100 and 1200°F. This depends on the 
chromium content of the surfaces (reduced concentra- 
tion pushes the critical temperature toward 1100°F). 
The effects of oxygen levels lower than 5 ppM and 
reaction kinetics are yet to be determined. This work 
also indicates that the change of performance with 
oxygen level may be even more sharp than shown in 





*Friction is assumed to be a sensitive indicator of the 
presence or absence of the complex oxide on the stainless-steel 
surface. 
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Fig. 6. These data serve to illustrate the importance of 
oxygen-level measurements and the absolute need for 
accuracy and reliability when defining sodium environ- 
ment for wear and friction effects. The difference 
between 3 and 5 ppM oxygen can be important. 

The relation shown in Fig. 6 can also be useful for 
analysis of test data. However, because of the com- 
plexity of reactions occurring at the sodium —stainless 
steel interface, it is not possible to define accurately 
the chemical reactions or predict the final surface 
chemistry. Extensive corrosion tests at MSA Research 
Corp.'* have resulted in a visible gray scale forming at 
1200°F on 316 stainless steel in sodium containing 
oxygen “‘in excess of 20 ppM,” and this scale has been 
identified as sodium chromite, Na,Cr,O,. At high 
temperature (1400°F), a scale of similar chemical 
content but of different phase was formed. 

Austenitic stainless steels in the hot regions of 
flowing-sodium circuits tend to lose chromium and 
nickel more rapidly than they lose iron. The resulting 
composition change in the surface produces a thin 
layer of ferrite on the order of several tenths of a mil in 
depth. The composition of this phase can go as high as 
90 wt.% iron. Flow velocity and oxygen content of the 
sodium appear to influence the selective attack. The 
same effect has been seen in Inconel 718 and the 
Stellites. Above 1350°F the leaching of nickel in 
nickel-based alloys occurs at a higher rate than in 
austenitic stainless steels.’ © 

The softer surface layer formed by selective re- 
moval of elements from an alloy exposed to flowing 
sodium will have a higher sensitivity to galling and 
adhesion than the bulk alloy during surface contact (P, 
in Fig. 4 is lowered). Therefore it is possible that, after 
exposure to flowing high-temperature sodium in the 
reactor, core components (selected on the basis of 
friction, wear, and adhesion characteristics for normal 
surface composition) will fail because of lowered 
resistance to adhesion and plastic flow. 

The effect of change in surface composition of 
alloys used in surface-contact applications in sodium 
may be influenced by wear and by the availability of 
contact surfaces to the flowing sodium. If the ferrite 
phase is worn away as fast as it is formed (and this is a 
distinct possibility and dependent on the kinetics of 
the process), very little influence on adhesion and 
galling characteristics of the alloy is expected. Further- 
more, most surface contacts—especially static flat-on- 
flat or pin-in-hole configurations—will inherently 
shield the contact areas from flowing sodium, and 
there is little likelihood of leaching in that condition. 
The most vulnerable configurations will be those in 
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which the contact surfaces are separated for long 
periods in flowing high-temperature sodium and are 
then brought back into sliding contact. 

The mass-transfer process responsible for the above 
effects can also result in deposits that could influence 
wear and self-welding characteristics of material com- 
binations. It is possible that these deposits could 
reduce wear and self-welding; however, if the deposits 
are loose and much harder than the surfaces them- 
selves, abrasive wear could result. No specific experi- 
mental data have been found in this general area. 

Wear and friction behavior, then, is definitely 
influenced by the above surface chemical effects. This 
also explains the wide divergence in wear data which 
has been the general experience of the past 10 to 15 
years. 

Using these effects deliberately for wear or friction 
reduction in reactor design is not practical as yet. This 
is because of the sensitivity of reactions to oxygen 
content, operating temperature, temperature gradients, 
and the possibility that surface reactions can proceed 
at too high a rate or result in undesirable surface phase 
changes (formation of ferrite). In addition, it is not 
expected that all complex metal oxides will be good 
boundary lubricants. These effects must be considered, 
however, in the interpretation of wear and friction data 
and in their application to component design. 


USING WEAR THEORY 
IN REACTOR DESIGN 


It is quite evident that a designer cannot select 
materials or adjust load conditions in a_high- 
temperature sodium environment on the basis of 
known material properties. Material selection can be 
governed to some extent on published wear and 
friction data in which the sodium environment has 
been fully defined (oxygen level, temperature, flowing 
vs. stagnant, and other contaminants), and the princi- 
ples concerning contact geometry, wear regime, etc., 
discussed earlier are adhered to. Perhaps the most 
useful application of the principles discussed here has 
been that for defining test conditions to obtain data 
for critical components. Since it is too expensive and 
time consuming to test full-scale reactor components,* 
the practical approach has been to develop test 
equipment in which the contact condition can be 
simulated closely (contact area, geometry, load, and 
sliding velocity) and the sodium environment is closely 





*Some tests have been carried out in Germany on fuel-rod 
bundles and spacer grids in sodium.'” 
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controlled and monitored. The use of the vanadium 
wire technique for measuring oxygen level has allowed 
wear experiments to be carried out at established 
oxygen levels close to 3 ppM. 

In the past 3 years, data have become available on a 
number of materials for wear and friction applications 
in liquid sodium. The data are more reliable than that 
developed in the late 50s and early 60s because the 
importance of sodium-environment control has been 
recognized. 

Some of the materials currently being evaluated are 
listed in Table 1. 

Surface-treatment techniques are currently under 
investigation for improving wear resistance of 
corrosion-resistant high-temperature alloys. Surface 
treatment is also being evaluated for reduction of 
friction levels. (Friction coefficients can range between 
0.5 to 1.0 in high-purity liquid sodium.) If the 
resistance to penetration is increased at temperature 
Pm in Eqs. 1, 2, and 7, reduction in wear and friction 
levels is anticipated. Essentially this means increasing 
the hot hardness of the surface material. This can be 
done by coating surfaces with hard carbides or oxides. 
In the past, this approach has not been too successful 
because of the tendency for hard coatings to fracture 
and spall off surfaces under vibratory contact. Ad- 
vanced coatings techniques, however, are developing 
better, tougher hard-surface coatings that warrant 
serious consideration. These techniques include chemi- 
cal vapor deposition, spark hardening, sputtering, and 
plasma spraying of carbides and selected metals. 


SUMMARY 


Wear and friction of mechanical components in the 
core region of liquid-metal fast breeder reactors 
(LMFBRs) warrant special consideration during design. 
The high-temperature liquid-sodium environment limits 
material selection to materials with limited resistance 
to adhesion and galling. Selection of materials and load 
limits currently is based on test programs involving test 
conditions simulating projected conditions in the reac- 
tor. The basic principles of friction and wear can be 
applied in setting up tests and interpreting data as long 
as the special conditions imposed by the reactor 
sodium environment are recognized. Boundary lubrica- 
tion in liquid sodium is influenced crucially by the 
combination of sodium temperature, oxygen content, 
and time of exposure. French and British researchers 
have both found that exposure time in 600°C 
(1112°F) sodium with less than 10 ppM oxygen 
content changes wear characteristics of chromium- 
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Table 1 Materials Under Evaluation for 1200°F 
Bearing Surfaces in LMFBR-Type Reactors 





Material Comments 





Has similar corrosion resistance to 
sodium as 304 or 316 stainless 
steel; at the same time, has less 
tendency to gall and self-weld. 
Over four times as hard as stain- 
less steel at 1200°F. Has shown 
continuing promise as a wear- 
resistant material among nickel- 
base alloys. 

An austenitic matrix containing 
titanium carbide particles. Much 
improved galling resistance over 
austenitic stainless steel. 

Not as good wear resistance as Inconel 
718 but shows ability to provide 
double-oxide boundary lubrication 
in sodium. Shows promise at light 
loads. Tendency for self-welding. 

LT-2 An aluminum oxide—tungsten— 
chromium cermet compatible 
with sodium. Shows good 
wear and antiseizure resis- 
tance, 

Cobalt-base alloys are very limited 
for reactor use, especially where 
wear debris may be liberated. 
Show good wear resistance in 
sodium. 

This alloy shows good resistance to 
self-welding but poor wear 
characteristics in sodium. 
Sodium molybdate reaction will 
provide good boundary lubrica- 
tion at temperatures below 
800°F. Sodium molybdate de- 
composes at higher temperatures. 

Although it shows excellent com- 
patibility with high-tempera- 
ture sodium, this material is 
notorious for its galling and 
seizure tendencies. However, 
this material is widely used in 
such reactor core components as 
fuel-element cladding and does 
experience contact conditions. 
Under light loads, it shows 
sensitivity to boundary lubrica- 
tion by complex oxide reaction. 


Inconel 718 


Ferro TiC 


Hastelloy C 


Stellite 6 


TZM molybdenum 


316 stainless steel, 
20% cold worked 





containing alloys. Recent experiments at LMEC. 
Canoga Park, Calif., have indicated that boundary 
lubrication effects can occur with as little as 3 ppM 
oxygen content in sodium in the appropriate tempera- 


ture range. 
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The transition from mild wear to severe wear is 
perhaps of more significance in reactor mechanisms 
than wear rates. Mild polishing-type wear or light 
scoring, for instance, can be tolerated in the fuel-rod 
separators (wire wrap) where contact with the stainless- 
steel cladding is concerned. Severe wear would damage 
or penetrate the cladding. 

Surface chemistry is equally important in frictional 
behavior in sodium. Friction levels can vary by an 
order of magnitude if surface roughness, reaction films, 
or contact geometry are altered. Hot hardness is 
another important factor influencing friction. Friction 
generally increases as a metal alloy is softened by 
increasing temperature. Thermal softening of substrate 
metal can also lead to breakup of hard coatings 
supported by it. 
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NEPA Impact on Nuclear Power Plants 


By Glenn Ellis* 


The National Environmental Policy Act of 1969 
(NEPA) enunciates a broad national policy on the 
natural environment and provides a mechanism to 
implement the policy.’ The NEPA legislation could be 
called a kind of bill of rights for the natural environ- 
ment; it puts the environment on an equal basis with 
the nation’s economic and technical priorities. New 
responsibilities for “encouraging productive harmony 
between man and his environment” are assigned to all 
the federal agencies and to the newly created Council 
on Environmental Quality. The basis for all NEPA 
decisions and actions is a detailed environmental report 
initiated by the applicant which must accompany “all 


federal actions” that “significantly affect the quality of 


the environment.” Federal actions affected include the 
AEC licensing of nuclear power reactors. NEPA be- 
came effective on Jan. 1, 1970. 


CALVERT CLIFFS DECISION 


On July 23, 1971, a U. S. Court of Appeals 
decision on the Calvert Cliffs nuclear power plant in 
the state of Maryland effectively broadened AEC’s 
NEPA responsibilities.2 The decision was basically that 
the AEC construed too narrowly its responsibilities 
under NEPA. The Court held that NEPA requires 
independent substantive review of a// environmental 
matters and makes the AEC directly responsible for 
evaluating the fofa/ environmental impact of nuclear 
power. The AEC had been relying heavily on various 
federal, state, and locally accepted environmental 
standards to evaluate certain environmental effects, 
such as thermal effects. The Court directed instead that 
the environmental costs in each individual case must be 





*Science Services Branch, Office of Information Services, 


U.S. Atomic Energy Commission. 


weighed against the expected benefits. Alternatives 
that affect the balancing of values in the analyses of 
costs vs. benefits must be taken into account. The 
Court also directed that AEC review all nuclear 
facility-licensing actions since Jan. 1, 1970, the effec- 
tive date of NEPA. The Court further required a NEPA 
review of construction permits that were issued before 
NEPA started, in cases where an operating license had 
not as yet been issued. For this review to be as 
effective as possible, the Commission was asked to 
consider a temporary halt in construction pending its 
review and the backfitting of technological innova- 
tions. 

Atomic Energy Commission Chairman James R. 
Schlesinger calls the Calvert Cliffs decision a “water- 
shed event” in the history of nuclear power.* The 
decision has had a major impact on the Commission’s 
program. Every nuclear power reactor must undergo a 
NEPA review except four which had full-term operat- 
ing licenses as of the date NEPA became effective. 
Altogether some 64 licensing applications for 95 
nuclear power plants are directly involved. Less af- 
fected are some five nuclear power plants for which 
operating licenses were issued after NEPA became 
effective and 10 provisional operating licenses issued 
before NEPA. 


NEPA REVIEW PROCEDURE 


The AEC responded to the Calvert Cliffs decision 
with a revised policy and procedure for implementing 
the broadened scope of NEPA, contained in the 
revision to Appendix D of 10 CFR 50, published in the 
Federal Register on Sept.9, 1971, and effective on 
that date.* The additional NEPA review needed for a 
license or permit depends on the status of the plant 
licensing and of the NEPA action already taken. The 
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status varies considerably around the nation for the 
100-odd nuclear power plants being planned and 
constructed. To accommodate the varying need, the 
Appendix specifies the procedure for NEPA review of 
three different status categories of nuclear power 
plants. The procedures actually vary little between the 
categories, the main difference being the more expedi- 
tious handling provided for NEPA reviews of plants 
either in the hearing process or nearing the time for 
hearings; a number of plants in this category are almost 
completed and ready to operate. The Appendix also 
provides authority to permit fuel loading and low- 
power operation of plants that are ready pending their 
NEPA approvals. 

Listed below are the three different status cate- 
gories of permits, licenses, or cases, and the approxi- 
mate number in each category that will require NEPA 
reviews; some of the cases involve more than one 
nuclear power plant: 





Number requiring 
NEPA review 


Status category 
in Appendix D 





Permits or licenses issued between 
Jan. 1, 1970, and Sept. 9, 1971: 

Construction permits 11 

Operating licenses 
Construction permits issued before 

Jan. 1, 1970, for which either an 

operating license or the notice of 

opportunity for hearing on an 

operating license had not been 

issued as of Sept. 9, 1971 ~30 
Cases in the process of hearings or 

for which hearings were pending 

as of Sept. 9, 1971 ~30 





The procedure basic to all the NEPA reviews—and 
generally typical of the procedure to implement other 
NEPA actions—can be briefly outlined: The review 
starts when the applicant submits the supplemental 
information AEC requires to implement the broadened 
NEPA procedures. This information supplements an 
initially filed Applicant’s Environmental Report 
(AER), an AER being necessary to initiate any NEPA 
action. Licensees were asked to submit their supple- 
mental AER by Nov. 9, 1971. They had been asked to 
show cause by Oct. 19, 1971, why their permit or 
license should not be suspended pending NEPA review. 
A hearing opportunity will be given an interested party 
upon the AEC decision to suspend a permit or license. 

The applicant’s supplemental AER must consider 
the environmental impact and any adverse effects of 
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the proposal to construct or operate a nuclear power - 
plant. Among other AER details, the applicant pre- 
pares a cost—benefit analysis of a quantitative nature 
which considers and balances the effects and available 
alternatives. The applicant also assesses the compliance 
of his facility with the various environmental-quality 
standards, including thermal and other water-quality 
standards. With the applicant’s report as the basis, the 
AEC prepares a supplemental-draft Detailed Environ- 
mental Statement (DES) of a scope comparable to the 
AER. Then the AEC invites comments on the DES and 
the AER by circulating them to other government 
agencies and to appropriate government officials with 
environment responsibilities on various levels and by 
making them available to the public. The AEC con- 
siders these comments in preparing a Final Detailed 
Statement (FDS) that, among other environmental 
considerations, includes a final cost—benefit analysis of 
a quantitative nature. The AEC also certifies whether 
the application should be approved, denied, or modi- 
fied to protect the environment. Copies of the AER, 
FDS, and comments received thereon are sent to the 
Council on Environmental Quality. A public notice 
from AEC appears in the Federal Register upon receipt 
of the applicant’s AER and upon completion of each 
of the two DES and FDS statements. The notice says 
the statements are available, requests comments, and 
specifies that comments will be publicly available when 
received and that comments will be considered. 

A hearing is then held before the Atomic Safety 
and Licensing Board (ASLB), at which time the AER 
and FDS are offered in evidence. Any interested party 
may take a position and offer evidence. The ASLB 
determines whether the NEPA requirements have been 
met, whether AEC’s review has been adequate, and 
decides any matters in controversy. Finally, the ASLB 
independently considers the balance among all the 
conflicting factors and determines whether the permit 
or license should be granted, denied, or modified to 
protect the environment. 


NEPA WORK LOAD 


A staggering work load faces the nuclear industry, 
the AEC, and other government agencies in implement- 
ing the broadened NEPA responsibilities. Both operat- 
ing licenses and construction permits require prior 
NEPA approvals. Many plants are either finished or just 
about completed, and each applicant is in or close to a 
hearing for an operating license. Some 11 of the new 
nuclear power plants have a combined capacity of 7.5 
million kW and vary from 87 to 100% of being 
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complete. Other plants are ready to be built, and the 
applicants are seeking construction permits. 


IMPACT OF DELAYS 


A comfortable healthy environment presumes an 
adequate uninterrupted supply of electric power to run 
the devices and machines upon which our technological 
society largely depends. Any significant delays in 
nuclear power-plant startup, for NEPA approval or any 
other reason, would reduce the electricity-generating 
capacity below the levels of peak demands with 
consequent electric-power shortages and frequent in- 
terruptions of power. The new nuclear power plants 
are all necessary just to keep up with the ever- 
increasing demand for more electricity. They have been 
scheduled in response to genuine need.° Virtually 
every one of these new nuclear units would come into 
scheduled service at about the latest time possible to 
maintain a reasonable margin of reserve capacity for 
peak load conditions. The problem is nationwide. 
Virtually no area of the country has any uncommitted 
surplus generating capacity at times of peak loads 
sufficient to supply significant amounts of firm power 
to shortage areas. William McGuire, Chairman of the 
National Electric Reliability Conference, representing 
the nine reliability regions that comprise the nation’s 
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electric-power systems, warns: “The impact on electric 
power could be catastrophic.”© John N. Nassikas, 
Chairman of the Federal Power Commission, cautions 
that delays would have “serious consequences adverse 
to the public interest, including detriment to the 
environment.”® Any significant delays in nuclear 
power-plant startups would thus introduce an adverse 
environmental impact that must be considered as part 
of the NEPA implementation. 
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AMERICAN NUCLEAR SOCIETY— 
CRITICAL REVIEWS 


The Atomic Energy Commission has contracted with the American Nuclear 
Society to prepare for publication on a regular basis detailed Critical Review 
articles written by experts selected from the ANS membership and reviewed 
by an advisory committee of top leaders and scientists in the field. In this 
issue the articles are on pages 16--75. 


Members of the Critical Review Advisory Committee are: 


Sidney Siegel, Chairman Atomics International 

William Chittenden Sargent & Lundy Engineers 
Frank G. Dawson Battelle—Northwest 

Don E. Ferguson Oak Ridge National Laboratory 
F. J. Leitz WADCO Corporation 

Paul Lottes Argonne National Laboratory 
Peter Murray Westinghouse Electric Corporation 
David Okrent University of Arizona 

Herbert Parker Battelle—Northwest 

A. M. Platt Battelle—Northwest 

Joseph Prestele New York Consolidated Edison 
W. C. Redman Argonne National Laboratory 
K. L. Rohde Idaho Nuclear Corporation 
Charles Stevenson Argonne National Laboratory 
Bertram Wolfe Battelle—Northwest 


ANS Critical Review Editor, Norman H. Jacobson 


ANS Officers: 
President, John W. Landis 
Vice President/President-Elect, James R. Lilienthal 
Treasurer, J. Ernest Wilkins, Jr. 


Executive Secretary, Octave J. Du Temple 


Comments on the articles should be communicated directly to the ANS 
Critical Review Advisory Committee or Editor, 244 E. Ogden Avenue, 
Hinsdale, Ill. 60521. 
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ANS GUIDELINES FOR PREPARING CRITICAL REVIEWS 
FOR REACTOR TECHNOLOGY 


One interpretation of the words ‘’Critical Review’’ emphasizes the word critical. That type of 
article would have as its purpose the discussion of a single subject, which has been uncertain and 
perhaps controversial, in considerable depth. An example might be a critical review of the values of 
alpha for plutonium, assessing all the work done and arriving at a best current estimate. This paper 
is primarily addressed to specialists on the particular subject, and serves as an authoritative source 
for the information they use. 

In the second interpretation of a Critical Review, the emphasis is on the word review, in the 
sense of survey. Such a paper would be broader and probably more descriptive in scope. A paper 
on ‘Solubility of Metallic Elements in Liquid Sodium” is an example of this category. Such a 
paper would be addressed to a much broader group of readers and written in a fashion that would 
be of interest to the majority of reactor technologists. It would provide enough of an introduction 
that the specialist from another field could immediately appreciate why the subject is important. 

The criteria for these two types of articles may not be very different: 

1. The paper should be based on a thorough coverage of relevant work on the subject from 
many sources. A review based on the work of only one individual or group is better suited for 
publication in one of the regular society journals. The review should be critically selective, 
reporting only the most valid results and indicating why some prior results have a questionable 
status, The review should call attention to significant gaps where more work is required in the 
subject of the article. 

2. The paper should be timely, on a subject of active current interest. 

3. The scope of articles acceptable as Critical Reviews includes all the subject areas identified 
for Reactor Technology (Economics, Physics, Mechanics, Construction, Fuel Elements, Fuel 
Cycles, Fluid and Thermal Technology, Fuel Processing, Components, Operating Performance), as 
well as Materials (including Source and Special Nuclear), Environmental Effects, and Effluent 
Management. Not desired are (1) reports of original research proposed for first publication and 
(2) review articles directed toward the specialist in the field of nuclear safety (which is covered by 
the AEC’s bimonthly review Nuclear Safety). 

4. The paper should be organized so that it is of immediate practical use to the readers. Such 
organization requires: attention to consistent use of units, presentation of important data in 
summary curves or tables, and a fully adequate bibliography to the original literature. Details on 
size, honorarium, style, etc., can be obtained from ANS. 
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Fast Breeder Fuel 


By P. Murray* 


Abstract: Fuel-pin irradiations in test reactors have demon- 
strated the high-performance capabilities of the mixed oxide- 
stainless steel system, Fission-gas release and fuel swelling can 
be accommodated, and the benefits arising from radiation- 
induced effects of swelling and creep in stainless-steel cladding 
give confidence that high-burnup targets for liquid-metal- 
cooled fast breeder reactor (LMFBR) fuel can be achieved. 
Fuel—cladding reaction has emerged as a potential problem 
that may limit peak cladding temperatures to which mixed- 
oxide fuel pins are designed and operated, Further detailed 
results on these limits are needed, particularly in the high-burn.- 
up (100,000 MWd/metric ton) range. Reactions caused by 
fuel—sodium interactions under dynamic coolant conditions do 


not appear to cause major operating problems, Definition of 


fuel-pin thermal-performance models requires additional ex- 
perimental data to minimize the uncertainties. 

Stainless-steel swelling and irradiation-induced creep have, 
in effect, focused attention on the fuel assembly and behavior 
of the fuel-pin spacer system and the flow-duct system, These 
are vital areas in which more precise materials-irradiation data 
and interpretation for design are necessary, Present uncer- 
tainties can be overcome by design, given more reliable 
irradiation data to the high fluences required for LMFBR 
application, Accommodation of swelling does, however, intro- 


duce an economic penalty that highlights the importance of 


longer-term metallurgical solutions to the swelling problem, 


In the United States, although the breeder reactor 
program has been vigorously supported for many years, 
recent environmental, energy-, and fuel-supply prob- 
lems and related studies have further reinforced the 
importance of expeditiously completing the develop- 
ment of the liquid-metal-cooled fast breeder reactor 
(LMFBR). There is today clear consensus among 
cognizant agencies of the federal government, utilities, 
and the industrial community as to the need for 
proceeding vigorously with the development of the 
breeder reactor as a major answer to the nation’s 





*Westinghouse Electric Corporation, Advanced Reactors 


Division, Waltz Mill Site, P.O. Box 158, Madison, Pa. 15663. 
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energy supply and related environmental problems. 
The President, in his Energy Message of June 4, 1971, 
emphasized the need to develop the breeder reactor 
and stressed the importance to the nation that the 
successful demonstration of the LMFBR be achieved 
by 1980. 

All the major industrial countries of the world are 
developing LMFBRs on a national priority basis, and 
each has committed significant resources to test facili- 
ties and demonstration plants. The Soviet Union is well 
along in construction of the world’s two largest 
LMFBRs of 350- and 600-MW(e) capacity, and the 
United Kingdom and France are scheduled to begin 
operation of 250-MW(e) plants in 1972 and 1973, 
respectively. West Germany is designing the SNR 
demonstration plant, and Japan has announced plans 
to construct a demonstration plant for operation in 
this decade. West Germany and France have also 
recently announced a joint plan to build two 
1000-MW(e) LMFBRs, the first to be located in 
France, and the second, later, in Germany. 

Approximately $500 million annually is being 
devoted worldwide to fast breeder development and 
the LMFBR projects shown in Table 1. 

The advantages of the fast breeder have been well 
documented in recent years.’ The key to its eco- 
nomic success lies in obtaining low fuel-cycle costs. 
Fuel-fabrication costs must be minimized since they 
are estimated to account for about 60% of the 
fuel-cycle costs. Burnup in the fast breeder must be 
maximized since fuel-cycle cost is inversely propor- 
tional to burnup. Many estimates of the components of 
the fuel-cycle-cost targets have been made, for exam- 
ple, in the Edison Electric Institute (EEI) fast breeder 
study.> In all the estimates the importance of fuel 
performance is clearly brought out.* In_ particular, 
burnup targets of 100,000 MWd/metric ton need to be 
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Table 1 LMFBR Projects 

















Pool or 
Name Country MW(t) MW(e) loop Operation 
Operable 
BR-5 Russia 5* Loop 1959 
Dounreay Great Britain 72 15 Loop 1959 
EBR-II United States 62.5 20 Pool 1963 
Fermi United States 200+ 67 Loop 1963 
Rapsodie France 40 Loop 1967 
SEFOR United States 20 Loop 1969 
BR-60 (BOR) Russia 60 12 Loop 1970 
Under Construction 
BN-350 Russia 1000 + 150 Loop 1972 
PFR Great Britain 600 250 Pool 1972 
Phenix France 600 250 Pool 1973 
FFTF United States 400 Loop 1974 
BN-600 Russia 1500 600 Pool 1976 
Planned 
KNK II West Germany 58 20 Loop 1973 
JOYO Japan 100 Loop 1974 
PEC Italy 140 Modified pool 1975 
SNR West Germany § 730 300 Loop 1976 
MONJU Japan 750 300 Loop 1976 
Demo 1 United States 750-1250 300—500 Not decided 1977 
Demo 2 United States 750-1250 300—500 Not decided 1979 
CFR United Kingdom SPS 1250 Not decided 1979 
1000 MW(e) France 2500 1000 Not decided 1980 
Demo 3 United States 750-1250 300—500 Not decided 1981 
1000 MW(e) Germany 2500 1000 Not decided 1982 
*To be increased to 10 MW(t). 
+Operated in 1970. 
+ Dual-purpose desalting and electric power. 
§ In cooperation with Belgium and the Netherlands. 
achieved for economic success. The magnitude of the Table 2 Relative Performance Criteria 
fuel problems that need to be solved can be seen from for LMFBR Fuels and Thermal Reactor Fuels 
the comparison of the relative performance criteria for 
LMFBR fuel and thermal reactor fuel given in Table 2. Required improvement 
Sats ; a? in fast breeder 
To meet these objectives, the continued operation 
‘ : a reactor fuel as com- 
of the Experimental Breeder Reactor II (EBR-II) and pared with thermal 
expeditious completion of the Fast Flux Test Facility Performance criterion reactor fuel 
(FFTF) are essential to the success of the LMFBR 
program. The LMFBR program needs not only the fast Burnup, MWd/metric ton 5 to 10 times 


Fuel power density, 
MW/ton 3 times 
Core power density, 


fluxes and testing capabilities that EBR-IIl now pro- 
vides but also the even higher fast fluxes and greater 


testing capabilities, particularly that of testing in closed kW/liter 10 times 
loops which the FFTF will provide. Vitally important Total neutron flux, 

irradiation testing of a variety of fuels, reactor control neutrons/(cm? )(sec) 30 times 
materials, and structural materials in a controlled and Neutron exposure or fluence, 

. ; ; ; neutrons/cm? 50 times 
instrumented fast-neutron flux which approximates Theensad thick. °Fhene $ os tae 


that required in economic fast breeder reactors will be 
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provided by the FFTF. In fact, the more recently 
discovered phenomena of stainless-steel swelling and 
irradiation creep reinforce the need for testing at the 
comparable fluxes and fluences. 

The purpose of this review is to concentrate on 
those factors which influence the irradiation perfor- 
mance of the fuel, fuel pin, and fuel assembly and to 
summarize the information available from the test 
reactors and irradiation programs in relation to the 
objectives defined for the LMFBR fuel. Attention is 
concentrated on the mixed oxide—stainless steel sys- 
tem since this system has been chosen for FFTF, for 
the LMFBR demonstration plants in the U. S. program, 
and also for the LMFBR plants being constructed 
overseas. 


FAST BREEDER FUEL-DESIGN 
PARAMETERS 


It was early recognized that the high power density 
necessary for fast breeder application, together with 
the use of sodium as coolant, required the use of 
subdivided fuel as compared with light-water fuel. 
Therefore fuel-pin diameters of ~0.25-in. diameter 
were considered at an early stage. Christensen® has 
pointed out that LMFBRs are, in fact, fuel temperature 
limited because of the constraints arising from the 
increased unit powers and the higher burnup. The 
primary factors determining the maximum allowable 
heat rating are thermal conductivity, solidus of the 
mixed-oxide fuel, radial power distribution in the fuel, 
and surface temperature of the fuel. On the basis of 
zero fuel melting for both cases, these factors and their 
effects result in a penalty of ~30% reduction in 
maximum allowable heat rating for an LMFBR relative 
to a light-water reactor, as shown in Table 3. 

Burnup effects, although substantial, include com- 
pensating factors. Restructuring of the fuel and its 
effects on thermal conductivity and the variation of 
the fuel—cladding gap AT with irradiation history are 
being measured® by in-reactor instrumented experi- 
ments in the EBR-II. 

A further effect recognized early was the high 
fission-gas release from mixed-oxide fuel irradiated to 
burnups’’® in excess of 25,000 MWd/metric ton. The 
interna. pressure generated within the pin and the 
resulting hoop stress in the cladding were considered 
initially to be the main failure mechanism at high 
burnup. 

As a result, most designers incorporated sufficient 
plenum space in the pin to limit the internal gas 
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Table 3. Comparison of Water-Cooled Power 
Reactors and LMFBRs in Terms of Factors That 
Affect Their Relative Limiting Heat Ratings 





Approximate percent 
changes in LMFBR 
allowable heat rating 
compared with 
water-cooled 


Differentiating factor power reactor, kW/ft 





Increased coolant temperature in 

the LMFBR aS 
Addition of large fractions of 

plutonium to LMFBR fuels 


Effect on thermal conductivity 8 
Effect on melting temperature 5 
Lack of flux depression in fast flux —§ 


Increased porosity for swelling 
accommodation in the LMFBR 


Effect of thermal conductivity —-12 
Restructuring effects +25 

Increased fuel—cladding AT because 
of higher heat fluxes in the LMFBR -10 


High-burnup effects in the LMFBR 
On thermal conductivity and melting 


temperature —10 
On fuel—cladding AT because of 
gap closure +10 
Net reduction in allowable heat 
rating for the LMFBR oxide fuel 3g 





pressure to reasonable levels (~800 psi) at end of life. 
The provision of a plenum has therefore been accepted 
universally, although the location at the top or bottom 
of the pin has varied in different designs. Thus, in U. S. 
designs, considerations of ejection of fission-gas bub- 
bles following failure and their effect on heat transfer 
in the fuel assembly have led to the choice of a top 
plenum. In the U.K. Prototype Fast Reactor (PFR), 
the plenum has been located in the bottom and cooler 
part of the pin to obtain the advantage of a smaller 
plenum length. The plenum location is also influenced 
by such nonfuel considerations as reactor-vessel height. 

In actual fact, the early irradiation tests? carried 
out with reasonable volumes of plenum space incorpo- 
rated in the pin showed that the swelling of the fuel 
and its interaction with the cladding was more likely to 
be a source of failure. Accordingly, the use of 
lower-density fuel to accommodate fuel swelling be- 
came both a major design consideration and a principal 
objective in the experimental programs. The definitive 
results obtained on this aspect of fuel performance in 
the EBR-II and their interpretation will be discussed 
later in this review. The method of incorporation of 
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porosity in the fuel has also resulted in variations and 
options. Low-bulk-density (85 to 90% of theoretical) 
solid pellets are the preferred choice in most designs 
leading to smear-density values (when the effect of the 
fuel—cladding gap is included) of 80 to 85%. Annular 
pellets are preferred in the PFR fuel because of the 
improvement in thermal performance, but there is 
some question as to the effectiveness of the central 
hole (as compared to distributed porosity in the 
low-bulk-density solid pellet) in accommodating fuel 
swelling. Consideration of cladding material has re- 
sulted in the choice of 316 stainless steel with a 
thickness of 15 to 16 mils in most cases. The peak 
hot-spot cladding temperatures arise not only from 
consideration of pin performance but also from plant 
optimization studies. In the United States, values in the 
range of 1250°F have been selected.!°°!! 

In summary, pin designs have centered on small 
diameters (~0.25 in. in diameter), provision of plenum 
space to minimize internal fission-gas pressure, and the 
use of low-smeared-density fuel to accommodate fuel 
swelling. Trade-off studies carried out by different 
organizations have resulted in some variation between 
the parameters actually selected. In general, these 
variations are small, and the correctness of the particu- 
lar choice will be determined only after statistical 
testing in the FFTF or demonstration reactors. 

In most fuel-assembly designs, the basic unit is a 
hexagonal duct, ~5 in. across the flats, containing 217 
pins. The spacer system supporting the pins is either 
wire wrap or grid, as will be discussed in detail later. 
Thus, in the FFTF design,’? the main component of 
the driver fuel assembly is the flow-duct region 
containing a close-packed hexagonal wire-wrapped 
array of 217 fuel pins, each having a diameter of 
0.230 in. The helical wire spacers provide radial sup- 
port for the pins and ensure space between them for 
coolant flow. The external structural member of the 
fuel assembly is the coolant-flow duct. Load pads are 
provided on the outsides of the ducts at elevations 
above the active core and at the tops of the ducts. 
These load pads separate adjacent ducts and provide 
bearing surfaces for restraining loads and interactions 
between ducts. The reactor core then consists of a 
hexagonal array of these vertical assemblies, the posi- 
tions and types of which are shown in Fig. 1. 

In the Westinghouse demonstration plant,'? the 
0.25-in. diameter fuel pins have a 16-mil wall thickness 
of 20% cold-worked 316 stainless steel and are 104 in. 
long, including 42 in. of fission-product-gas plenum. 
Spacing is by honeycomb grid in a 217-pin hexagonal 
duct. Average and peak design burnups are 75,000 and 
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103,000 MWd/metric ton, respectively, representing a 
peak linear power rating of 15.1 kW/ft and a maximum 
fuel temperature of 4840°F. 

The Atomics International (Al) demonstration 
plant'* has a fuel-pin diameter of 0.27 in., 80% smear 
density, 85% pellet density, with 271 pins wire 
wrapped in the hexagonal duct. 

In a similar way to the appreciation of the effects 
of fission-gas release and swelling of fuel in the design 
of fuel pins, the advent of irradiation-induced swelling 
and creep in stainless-steel cladding and core structural 
materials has had considerable impact. The implica- 
tions with respect to design and performance have been 
the major task over the past 3 years. Three principal 
areas were identified: 

1. The effect on fuel-pin performance. 

2. The effect on irradiation testing programs. 

3. The effect on the spacer support systems and the 
duct. 


Items 1 and 3 will be discussed in detail in later 
sections of this review. Analysis of the implications on 
irradiation testing of fuel pins and fuel assemblies in 
present test reactors led to the early realization of the 
importance of fluence-to-burnup ratio. Most previous 
irradiation tests had concentrated on obtaining high 
burnup in the fuel in the shortest possible time by 
overenrichment, leading to a low fluence-to-burnup 
ratio. The new irradiation effects in the stainless steel 
at the higher fluxes and fluences for the LMFBR meant 
that, for integral experiments, much higher fluence-to- 
burnup ratios were needed. In Table 4 the fluence-to- 
burnup ratios are compared for the FFTF fuel-perfor- 
mance parameters and the irradiation experience for 
highly enriched fuel. 

As a result of these and similar evaluations, extra 
emphasis was placed on programs involving irradiation 
of nonfueled cladding and structural materials to high 
fluence. The fuel-irradiation programs were reoriented 
to take account of this effect as far as possible. True 
matching of the effects of irradiation and burnup on 
fuel, cladding, and structural materials can only come 
from the integral driver-assembly behavior in FFTF. 


IRRADIATION EXPERIENCE IN 
TEST REACTORS 


Recognition of the magnitude of the problems that 
had to be solved for fast breeder fuel has resulted in 
concentration on a wide variety of irradiation tests on 
fuel pins, fuel assemblies, cladding, and structural 
materials in fast test reactors over the past 10 years. 
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Note: Reactor Capability, 
6 Closed Loops 


Special-purpose test, 2 


Open test assemblies, 3 


instrumentation, 1 


Oscillator, 1 


W BE) SD €9 


General-purpose closed loops, 2 


Open test assembly with proximity 


Alternate control-rod positions, 2 


P. Murray 


Drivers, 28 inner (Rows 1—4) 
and 47 outer (Rows 5-6) 


Safety rods, 3 


Peripheral control rods, 8 
Fixed shim rods, 7 


Reflectors, 93 


CR In-core shim/scram rods, 4 


Fig. 1 FFTF core map. 


The principal fast test reactors used are the EBR-II in 
the United States, the Dounreay Fast Reactor (DFR) 
in the United Kingdom, the Rapsodie in France, and 
the BR-5 in the USSR. Significant results have been 
achieved over the past 5 years, and these are best 
described under the individual test reactor. The limita- 
tions of existing fast test reactors, particularly with 
respect to flux, also led to the commitment to the 
FFTF in the United States, as outlined earlier. In 
addition, provision has been made in the PFR and 


REACTOR TECHNOLOGY, Vol. 15, No. 1, Spring 1972 


Phenix demonstration plants for the inclusion of test 
assemblies that are primarily variants of the reference 
fuel in these reactors. 


Experimental Breeder Reactor II 


The EBR-II was originally proposed and designed 
to demonstrate the engineering and operating feasi- 
bility of the fast reactor for central-station power-plant 
application. By the time the plant was in operation, the 
long-range national emphasis had shifted to much 
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Table 4 Comparison of FFTF Fuel-Performance Parameters 
and Irradiation Experience on Mixed-Oxide Fuel Pins 








Cladding Fast flux, Fast fluence, Fluence/burnup 
temp. Power level, Burnup, 10° neutrons/(cm? )(sec) 107? neutrons/cm? ratio, 
Region range, °F kW/ft at.% (E > 0.1 MeV) (E > 0.1 MeV) 10? /at.% burnup 
FFTF Fuel-Performance Parameters 

Inlet 660-715 5-8 4-6 1.6—2.9 4.5—7.8 1.1-1.3 
Core 805-970 12.4 8.6 4.5 4255 1.4 

Outlet 930-1150 5-8 4-6 1.6-2.9 4,5-7.8 1.1-1.3 

Irradiation Experience 

Inlet 660-715 9-15 3-—7* 1-2 2-6* 0.5—0.9 
Core 805-970 8-15 3-7 1-2 2-6 0.5—0.9 
Outlet 930-1150 7-15 3-7 1-2 2-6 0.5—0.9 





*Burnup and fluence in the EBR-II have now reached 11 at.% and 8 x 107? neutrons/cm?, as will be discussed later. 


larger plants using dilute ceramic fuel, and the mission 
of the EBR-II was changed to that of a primary 
irradiation facility for the LMFBR program. Because 
the EBR-II was not designed as an irradiation facility, 
the transition required changes in several original 
concepts of operation as follows: 

e An increase in overall core size to accommodate 
the irradiation experiments. 

e Changes in fuel enrichment and blanket composi- 
tion. 

e A shift of operating philosophy from that pri- 
marily of an engineering test facility to that of a 
high-priority neutron producer. 

e The design of a variety of irradiation vehicles for 
use by experimenters. 


e The development of methods to monitor the 
changing nuclear conditions of a heterogeneous reac- 
tor. 

e The establishment of operational reviews and 
procedures that would allow for the acceptance of the 
increased risk that must accompany any experimental 
program and still not compromise the availability of 
the reactor in its mission as an irradiation facility. 


The successful operation of the EBR-II as an 
irradiation facility with a high availability is evident 
from the record of its operating history and from the 
number and variety of experimental irradiations 
already accomplished, as well as those presently under- 
going tests in the reactor (Table 5). A summary of the 


Table 5 EBR-II Conversion to an Irradiation-Test Capability 





1968 1969 1970 1971 





Plant-Operation Improvements 


Calendar days of power operation 
At 45 MW(t) 
At 50 MW(t) 
At 62.5 MW(t) 

Plant factor, % 


Irradiation Experiments 


Total experiments since operation began 
Encapsulated fuel 
Unencapsulated fuel 
Cladding and structural materials 
Absorber materials 


74 

55 129 183 
9 39 143 
Bs 38 61 39 
367 440 589 675 
72 368 747 1611 
248 312 350 557 
7 11 15 38 
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Table 6 Number of Fuel Pins Irradiated in EBR-II* 
Up to 5 at.% 5 to 10 at.% Beyond 10 at.% Total 
Encap- Unencap- Encap- Unencap- Encap- Unencap- Encap- Unencap- 

Fuel sulated sulated sulated sulated sulated sulated sulated sulated 
Oxide 147 559 127 241 39+ 0 313 800 
Carbide 49 0 qs 0 8 0 128 0 
Nitride 18 0 0 0 0 0 18 0 
Metal + 

cermet 70 196 4 38 0 0 74 234 
Total 284 755 202 279 47 0 533 1034 





Number of Capsules with Structural Materials Irradiated in EBR-II* 





2 to 4 x 107? neutrons/cm? 
(E > 0.1 MeV) 


Up to 2 x 10? neutrons/cm? 
(E > 0.1 MeV) 


4 to 6 x 10*? neutrons/cm? 


Beyond 6 x 107? neutrons/cm? 


(E > 0.1 MeV) (E > 0.1 MeV) 





262 143 


32 37 





*Status as of December 1971. 


+Peak cladding fluences up to 0.95 x 10? * neutrons/cm? (E > 0.1 MeV) have been achieved. 


number of fuel pins irradiated in the EBR-II and of the 
burnup levels achieved is shown in Table 6. 


The highlights of the EBR-II irradiation results may 
be detailed as follows: 


1.A significant number of encapsulated mixed- 
oxide fuel pins have been examined'!* recently after 
irradiation to burnups of over 100,000 MWd/metric 
ton in the EBR-II. Results from these experiments have 
significantly increased confidence that UO,—PuO, fuel 
pins and subassemblies can be designed to achieve 
required economic and performance targets in a safe, 
reliable manner. To date, most of these irradiations 
have been carried out at cladding temperatures 100 to 
150°F lower than the optimum required for the 
LMFBR. Pins now in the reactor are being irradiated in 
the higher temperature range. High-burnup, mixed- 
oxide fuel pins have successfully attained target expo- 
sures of from 80,000 to over 100,000 MWd/ metric ton 
{Nuclear Materials & Equipment Corp. (NUMEC) 
Series C pins and General Electric Series F2 pins]. 
Postirradiation examinations indicate that the maxi- 
mum diameter increases of 316 and 347 stainless-steel- 
clad fuel pins ranged from 2.2 to 5.6% at peak fuel 
burnups to ~12at.% and total fluences up to 
~1.0 x 107? neutrons/cm?. Initial smeared fuel den- 
sity has a major effect on the diametral increases of 
these fuel pins. Figure 2 shows the peak fuel-pin AD/D 
data as a function of peak burnup and smeared fuel 
density. This figure illustrates the major change in 
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diameter observed as smeared fuel density increases 
and is clear evidence of the effect of initial porosity in 
accommodating fuel swelling. 

2.Interim and terminal postirradiation examina- 
tions of assemblies of fuel pins irradiated in the 
EBR-II, in contact with the EBR-II primary sodium 
coolant, to exposures over 50,000 MWd/metric ton 
have been completed. No failures occurred. Measure- 
ments on fuel pins clad with 304, 304L, and 316 
stainless steels at ~30,000 and ~50,000 MWd/metric 
ton indicate that the 316-clad pins swell less (~50%) 
than the 304- or 304L-clad pins. 

















32 
ue | | | | | | | | | 
a 
6 6 High burnup — 
a 
Ww /_—— — 
oO 
2 
7 
rs} 4h a 
ais 
<x = = 
cc 
be 
= 2 

= ° =! 
2 é Medium 
- ° 0. burnu 
. o of ons 
. eo a wee 
A ie «: ae 
5 80 84 88 92 96 100 
Ww 


FUEL SMEARED DENSITY, % 


Fig. 2 Correlation between fuel-pin diameter change, burnup, 
and fuel smeared density. 
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3. Fuel—cladding reactions are being studied exten- 
sively, as will be discussed later in this review. The 
magnitude of the reaction at the inside diameter of the 
cladding is not significant at cladding inside-diameter 
temperatures <1100°F (0.5 to 1.0 mil); however, at 
higher temperatures, >1250°F, extensive reactions 
have been observed (2.0 to 4.0 mils). 

4.A considerable program is now under way on 
FFTF driver-pin and wire-wrapped subassemblies.' ° 
Three 37-pin EBR-II subassemblies, PNL-3, -4, and -5, 
were designed with fluence-to-burnup ratio as a major 
parameter. The PNL-3 subassembly has a cladding 
fluence-to-fuel burnup ratio typical of a mixed-oxide- 
fueled fast reactor (FFTF) obtained through the use of 


25% PuO,—75% *2%UO, fuel. The goal exposure of 
a p 


this test in EBR-II was predicated on cladding swelling 
to 1.1% AD/D and subsequent mechanical interaction 
between the fuel-pin bundle and the outer coolant 
duct, rather than on fuel burnup. The burnup calcu- 
lated to correspond to 1.1% cladding AD/D is 
30,000 MWd/metric ton, and the cladding fluence is 
5.8 x 107? neutrons/cm? (E > 0.1 MeV). 

The PNL-5 subassembly has satisfactorily reached 
initial goal burnup of 50,000 MWd/metric ton. Follow- 
ing disassembly and inspection of the pins, selected 
pins will be destructively examined. The remaining pins 
will be incorporated into fail-safe Mark-E subassemblies 
as part of a run-to-failure test series. 

The plan is to irradiate one Mark-E-19 subassembly 
until a pin failure occurs. The pins in another Mark- 
E-19 subassembly will be subjected to interim examina- 
tion at incremental burnups to observe the onset of 
fuel—cladding mechanical interaction and incipient 
failure. 

5. During the FFTF conceptual design phase, 20% 
cold-worked 316 stainless steel was selected as the 
reference fuel-pin cladding, and a set of reference fuel 
parameters was established. Accordingly, the PNL-6, 
-7, and -8 experiments test a second iteration of the 
FFTF fuel-pin design. In this series, emphasis is placed 
on testing fuel that has been statistically sampled for 
such characteristics as pellet diameter and density. 

The latest iteration of the FFTF driver-fuel-pin 
design is embodied in three 61-pin subassemblies, 
PNL-9, -10, and -11. The PNL-9, -10, and -11 fuel-pin 
design features are: 0.230-in.-OD by 0.200-in.-ID 20% 
cold-worked 316 stainless-steel cladding; planar fuel 
smeared density 85.5% of theoretical; 0.194-in.- 
diameter dished end pellets; and 0.006-in. fuel- 
cladding diametral gap. 

6.In addition to the fuel-pin-design test subas- 
semblies described, irradiation of an instrumented 
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subassembly, PNL-17, designed to test the thermal— 
hydraulic calculations used in the design of noninstru- 
mented subassemblies, was recently completed in the 
EBR-II. 


7. The first FFTF  fuel-development-program 
TREAT transient test of a fast-flux-irradiated, stain- 
less-steel-clad mixed-oxide fuel pin demonstrated an 
unexpected ductility, ~/4,% AD/D, without failure of 
the cladding. A second transient test of a fast-flux- 
irradiated pin was conducted at a greater energy input 
to yield gross failure of the pin. Preliminary postir- 
radiation evaluation confirms that the upper energy 
limit for pin-cladding integrity under transient over- 
power conditions is consistent with design safety 
criteria being used for FFTF driver fuel. 

TREAT tests of preirradiated, prototypic Fast Test 
Reactor (FTR)-length fuel-column pins are also being 
conducted. 

8. Detailed information is now being obtained on 
the kinetics of in-reactor restructuring of the fuel!” 
and on the spatial distribution of porosity in mixed- 
oxide fuel during irradiation.’ ® 

9. As indicated, reactions have been observed at the 
fuel—cladding interface in fuel pins irradiated in both 
fast- and thermal-flux reactors. The austenitic stain- 
less-steel-clad wall adjacent to the fuel reacts with 
fission products and fuel. Constituents of the cladding 
are transported into the fuel matrix, and intergranular 
attack of the cladding is observed. Components of the 
reaction phase have been identified by electron-probe 
microanalysis and include uranium and plutonium 
from the fuel; cladding alloy elements; impurities from 
the fuel or cladding, such as silicon; and a number of 
fission products—cesium, iodine, barium, tellurium, 
palladium, and molybdenum. Strontium and zirconium 
have also been observed in the fuel near the cladding 
interface. Cesium, iodine, tellurium, and silicon have 
been detected in the grain boundaries of the cladding. 

The extent of cladding penetration has been 
measured on metallographic sections from irradiated 
fuel pins, and a summary of General Electric results’ ? 
is shown in Fig. 3. A correlation of increased penetra- 
tion with increasing temperature is clearly evident. In 
addition, there is the indication that the least reaction 
is observed with lower oxygen-to-metal (O/M) ratios.?° 

Cladding penetration to the extent of about 20% of 
the 0.015-in. wall thickness has been observed in pins 
operating at peak LMFBR temperatures with the 
penetration resulting from the effects of chemical 
reactions of fuel, fission products, and cladding, and 
from intergranular attack of cladding. 
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Fig. 3 Effect of temperature on the fuel—cladding reaction. 


The presently available fuel—cladding reaction data 
have the following characteristics: 

a. With the exception of the F12 pins (shown in 
Fig. 3), which may have operated with fuel center 
melting, all data with cladding temperatures >1200°F 
were obtained from thermal-flux tests with fuel having 
O/M ratios of 2.00. 

b. Insufficient data are available to assess quantita- 
tively the effect of cladding type, fuel O/M ratio, 
fuel-impurity content, fuel-fabrication technique, and 
cladding stress; but there is a clear trend indicating 
that the reaction is reduced the lower the O/M ratio. 

c. The high-temperature tests were all run for 
~10,000 hr to average fuel burnups of 
~140,000 MWd/metric ton. Thus end-of-life data are 
available for only a limited time and burnup range. It is 
quite probable that the reaction rate decreases with 
time. 

d. In addition to differences in fission yields, the 
thermal-flux depression causes a large fraction of the 
fission events to occur adjacent to the cladding. This 
provides an unrealistic concentration of fission prod- 
ucts adjacent to the cladding. 


A considerable program is now under way in the 
EBR.-II to obtain more quantitative information and to 
determine the limiting conditions imposed by fuel— 
cladding chemical interaction. 
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testing pins in trefoil experiments? and for 77-pin 
grid-type subassemblies irradiated in the core center. 
The reactor is downflow and cooled with NaK, and 
burnups of ~75,000 MWd/ton have been achieved on 
pins particularly under subassembly conditions at peak 
hot-spot cladding temperatures up to 1290°F. Cover- 
gas entrainment in the coolant has led to some 
problems of overheating, causing the failure of about 
50 pins out of approximately 800 irradiated to date. 
The experience gained not only has led to confidence 
in the steady-state performance of fuel at high cladding 
temperatures for the PFR initial target burnup but also 
has provided experimental evidence that fuel-failure 
propagation is highly unlikely under these irradiation 
conditions. 

Considerable data on the performance of mixed- 
oxide pins have recently been reported???* with 
particular emphasis on the swelling of the stainless-steel 
cladding. Data are presented on solution-treated and 
cold-worked M316 and M316L, aged Nimonic PE16, 
and FV548. Figures 4 and 5 illustrate the results on 
M316 and M316L cladding. These results show that the 
cladding-volume changes account for most of the 
observed diametral changes in both solution-treated 
and cold-worked materials. In a later section of this 
article, a discussion of stress effects on stainless-steel 
swelling suggests that some of the volume change may 
be stress induced. 

Limited results are available on six of the 77-pin 
grid-type subassemblies irradiated in the DFR. Figure 6 
shows the total diametral strain for cold-worked M316 
cladding from the Mark 1, 1B, 2, 2A, 2B, and 2C 
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Fig. 4 Relation of total diametral strain to predicted strains 
attributable to swelling for solution-treated M316 and M316L 
steels, 


subassemblies. Figure 6 also shows data to 10 at.% 
burnup obtained in trefoil experiments. The small 
diametral strains reported for Nimonic PE16 are of 
particular interest, as will be discussed later. 

The extensive experimental evidence from the 
United Kingdom with peak fluence up to ~8 x 107? 
neutrons/cm? shows: 

1. Diametral strains and cladding-volume changes 
are smaller in cold-worked M316 cladding than in 
solution-treated material. 

2.The presence of finely dispersed gamma-prime 
particles in a precipitation hardening alloy, such as 
Nimonic PE16, provides a means of suppressing swell- 
ing. 


Irradiation-creep experiments on helical springs 
conducted in the DFR have yielded important informa- 
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Fig.5 Relation of total diametral strain to predicted strains 
attributable to swelling for cold-worked M316 and M316L 
steels, 


tion.25 Data are available on M316, Nimonic PE16, 
Nimonic 80A, nickel, and molybdenum. The results 
obtained at ~480°F indicate that irradiation creep is 
linear with respect to stress and fluence. The Nimonic 
PE16 springs exhibited the lowest level of irradiation 
creep. This result, coupled with the low swelling 
reported earlier for this alloy, has important theoretical 
and technological implications. 


Rapsodie 


Since Rapsodie went to power in March 1967, a 
surveillance program has been conducted to establish 
the behavior and performance of mixed-oxide fuel 
elements. A milestone was reached in the French 
sodium-cooled fast reactor program when the initial 
target burnup of 35,000 MWd/metric ton was achieved 
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Fig. 6 Variation of pin-diameter change with burnup. 


in the central assemblies of the core. As burnup 
progressed, selected standard fuel assemblies were 
discharged from the core during periodic reactor 
shutdowns. 


The standard assembly contains a hexagonal bundle 
of 37 fuel pins. The reference pin, 48.5 cm long, has a 
34-cm column of sintered pellets 96% dense with the 
composition 26 wt.% PuO,—74 wt.% UO, (60% en- 
riched), and 1.98 stoichiometry. Type 316 stainless- 
steel-clad tubing, 5.80mm in inside diameter and 
6.70 mm in outside diameter, has an integral helical 
fin, and the initial diametral fuel—cladding gap?°’*” is 
0.23mm. The schedule of examination and _fuel- 
irradiation parameters are presented in Table 7. 


Detailed examination and precise nondestructive 
measurement of complete assemblies and of all fuel 
pins did not indicate any sign of cladding failure, 
bowing distortion, or weight or length change.”*® 
Continuous measurements of pin diameter at 90° 
orientation from the fin have shown no diameter 
increase within selection criteria of cladding tubing 
batches, Whatever the irradiation level, epithermal- 
neutron radiographs of pins showed permanence of 
interpellet spacings and appearance of two crack 
patterns in the fuel column. Histograms of fuel-column 
length vs. burnup indicate maximum elongation of 
0.5%; however, interpellet spacings and transverse 
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fractures make any determination of swelling rate 
difficult. 

Postirradiation examination of metallographic sec- 
tions indicated a general trend in fuel-microstructure 
evolution. 

1. Operation at power, 335 W/cm: No gap closure 
or grain growth (except small elongated grains at the 
edge of radial cracks near the center). 

2. Operation at linear power between 335 and 
375 W/cm: Gap closure, outer ring of as-fabricated 
structure surrounding equiaxed growth zone (porosity 
at grain boundaries) during steady-state operation. 


When the fuel sticks to the cladding, a new 
circumferential crack pattern is superimposed on the 
radial cracks. 

Fission-gas determination performed on a limited 
number of pins in each assembly indicates a drastic 
increase from 10 to 55% release with burnup. 

Pressure burst tests have also been performed on 
cladding tubings after fuel removal, and a 2 to 3% 
residual ductility was measured after a fluence of 107? 
neutrons/cm?. 

Although the reactor sodium was never contami- 
nated, activity bursts of xenon in argon cover gas 
occurred during reactor operation, and minor cladding 
failures have occurred. At the burnup levels to date, 
fission-gas release and cladding properties, rather than 
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lable 7 Operational Conditions and Results 





Reactor cycle number 





1 2 3 A 5 

Nominal reactor power, MW 20 20 20 24 24 
Local peak linear power 

(Row 01), W/cm 335 375 435 400 400 
Cumulated peak burnup level, 

MWd/ton 150 2,300 10,400 20,000 31,200 
Radial gap closure No Partial Full Fuel stuck to clad- 

ding 

Fission-gas release, % 1 10 50-55 55-60 





fuel swelling, are considered to be the major problems 
in the performance of Rapsodie driver fuel. 


Nondriver Fuel Testing in Rapsodie? ° 


The objective of other fuels testing in Rapsodie is 
to resolve various design questions for the Phenix fuel 
element. Among several important differences between 
Rapsodie and Phenix fuel, these tests are intended to 
investigate the following: 

1. Burnup levels of over 80,000 MWd/metric ton 
(Rapsodie fuel pins designed for 35,000 MWd/metric 
ton). 

2. Smeared densities in the range 75 to 80% as 
compared with 89% for Rapsodie driver fuel. 

3. Linear power ratings about 10% higher than 
Rapsodie fuel. 

4.Hot-spot temperatures of 1290°F or higher 
(compared with 1200°F for Rapsodie). 

5. Lower axial blanket integrated into fuel pin and 
the fission-gas plenum at the bottom of fuel pin. 

6. Wire-wrap spacer instead of the integral fin. 


In addition, there will be some testing of annular 
rather than solid pellets, and some cladding-material 
variations will be examined, particularly Hastelloy X. 
Two parameters that cannot be tested are the longer 
height of the Phenix subassembly (850 vs. 340 mm) 
and the larger size of the assembly. 

The two major efforts to date for investigating 
these variables consist of the Aurore and Orphee tests. 
The Aurore program includes six assemblies that are to 
be examined in pairs at burnup levels of ~40,000, 
60,000, and 80.000 MWd/metric ton. The diameter of 
the Aurore pins has been increased to 7.15 mm 
(compared with 6.6 for Phenix) to achieve the desired 
linear power and thermal conditions. This yields linear 
power levels of about 480 W/cm, maximum nominal 
cladding temperatures of about 1200°F, and hot-spot 


temperatures of about 1330°F. Because of the larger 
pin diameter, the subassembly includes 27 pins ar- 
ranged in three groups of nine (dividing the hexagonal 
subassembly into three equal regions). 

The Aurore experiments serve to provide informa- 
tion on most of the aspects needed for the Phenix 
subassembly. The major variations are the shorter -pin 
length, larger pin diameter, and separation of the lower 
blanket from the fuel-pin assembly. 

The Orphee tests, with higher fuel enrichments, 
accomplish all the conditions that are achieved in 
Aurore, but, in addition, the fuel-pin diameter is 
retained at 6.6mm and the lower axial blanket is 
integrated into the fuel pin with plenum space at the 
bottom. Very few details of the Orphee tests are 
available, but it is believed that these include refractory 
cladding and possibly other fuel-pin design variations. 

Examination of some Aurore subassemblies has 
already taken place. Behavior of the fuel and cladding 
to a burnup level of about 60,000 MWd/ metric ton was 
reported to be extremely good, and levels up to 
90,000 MWd/metric ton have been reached recently. 


Cladding- and Structural-Material 
Testing 


Testing of cladding can be divided into two broad 
categories: 

1. Experiments which are placed in the core and 
which are basically fissile driver elements with a 
portion of the assembly replaced by cladding samples. 

2. Experiments which replace some of the pins in 
the fertile, or radial, blanket subassemblies and which 
are at the periphery of the fast region of the core. 


In addition, there are elements for investigations of 
different procedures for fixing wire wraps to the fuel 
pins, elements in which the cladding experiment is 
removable so that it can remain in the core for periods 
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in excess of driver-fuel lifetimes, and finally there are 
capsules located in the reflector for structural-material 
irradiation. 

The first class of cladding experiments includes 
cladding tensile specimens located in flowing sodium, 
experiments to investigate cladding-creep characteris- 
tics, cladding specimens with an insulating gas space to 
permit control of cladding temperature, tubular clad- 
ding specimens with high internal pressures (and 
contained in special, pressure-safe capsules), and clad- 
ding specimens with a small (fuel) reheater section to 
permit higher-than-normal temperature levels. Al- 
though many of the cladding tests have purposes 
specific to Orphee or Aurore testing, most attempt to 
provide data on a broad range of properties such as 
uniaxial and biaxial creep, embrittlement, swelling, and 
on the importance of fabrication variations (annealing, 
cold working, and grain size). 

In general, it can be said that fuels and cladding 
experience from Rapsodie to date has been very good, 
with results showing reasonable agreement with the 
experience of other workers. To date, a large number 
of pins have been satisfactorily irradiated to high 
burnup. French experience has yielded extensive de- 
tailed results concerning mixed-oxide fuel behavior to 
burnup levels of over 60,000 MWd/metric ton. 


Fortissimo® ° 


The Fortissimo project is an outgrowth of the 
program for the second core loading of Rapsodie, with 
the added objective of providing as great an increase in 
the flux level as possible, without major modifications 
of the system. Although it is desirable to increase the 
ratio of flux to power as much as possible, criticality 
requirements severely limit any reductions in enrich- 
ment or volume (which would serve to increase the 
flux-to-power ratio). As a result, an increase in power 
level from 24 to 37 MW (approximately ) was accepted 
to increase the flux from its present level of 
1.9x 10'5 to 3x 10!* neutrons/(cm? (sec). The in- 
crease in power level is accommodated by increasing 
the sodium flow from 890 to about 1100 m*/hr, and 
by increasing the AT through the core from 190 to 
250°F (allowing a hot-spot cladding temperature of 
1290°F vs. the present 1200°F). 


Driver Fuel-Element Modifications 


The major design change concerns the driver fuel 
element, which was modified to double its power 
output without any increase in linear power rating. 
This is achieved by using a smaller diameter fuel pin. 
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Addition of another ring of fuel pins per subassembly 
increases the number of pins from 37 to 61, with the 
pin diameter reduced from 6.7 to 5.1 mm. Linear 
power rating is maintained at about 400 W/cm, which 
is very near the current levels in Rapsodie. Cladding is 
316 stainless steel with a wire-wrap spacer. Gas-expan- 
sion space is provided at the top and bottom of the 
pin. The mixed-oxide fuel pellets are solid, sintered 
pellets with a density of about 92% of theoretical and a 
diameter of 4.2 mm. Reducing the number of driver 
subassemblies provides space for materials capsules, 
with the reactivity decrease being compensated by a 
modest increase in enrichment (30% plutonium and 
85% 735U, compared to the present 26% plutonium 
and 60% ?3*U). 

Hot-spot temperature and flow and pressure-drop 
limitations result in a reduction in core height from 34 
to 32 cm. Pressure-drop limitations and potential cavi- 
tation at the subassembly outlet have also resulted in 
the elimination of orificing of the subassemblies, 
providing a higher flow rate through all the driver fuel 
elements. 

The most notable feature of the now completed 
Fortissimo project was the decision to increase the flux 
at the expense of obtaining experience on a practical 
driver fuel-element design (which more closely approxi- 
mates Phenix-type fuel). The French decision was 
related to the fact that their experience, as elsewhere, 
shows that cladding problems (swelling, creep, and loss 
of ductility) are likely to be more difficult to resolve 
than fuel problems (fuel swelling and restructuring). 
An even more compelling reason for the French 
decision is the Phenix schedule, which requires early 
firm decisions on fuel and cladding design specifica- 
tions. Fortissimo will reduce the time required for 
acquisition of fuel information from about | year to 
9 months and will reduce the duration of cladding 
experiments (to the desired levels of integrated flux for 
initial Phenix conditions) from about 2% years to 
about 1% years. 


BR-5 


The BR-5 [5 MW(t)] has been the principal USSR 
fast-flux testing facility and operated from startup in 
1960 with a PuO, core until May 1965 when a 
235U-enriched uranium carbide core was installed. 
During this campaign, 37 of the PuO, subassemblies 
failed, and the experience obtained constitutes the 
largest to date on failed fast breeder oxide fuel. 

Examination of the 37 failed oxide assemblies after 
6 at.% burnup revealed only one or two failed pins per 
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assembly, and one of the highest burnup oxide 
subassemblies has been examined in detail.?! This 
subassembly had been irradiated for 5 years in the 
BR-S to a burnup of 6.1 at.% and a fluence of 
2.4 x 107? neutrons/em? at a maximum cladding 
temperature of 1080°F. At the start of the irradiation, 
the calculated oxide center temperature was 2550°F; 
at the end of the run, 2100°F. 

The fuel pin in the BR-5 PuO, assembly (19 
wire-wrapped pins) has a 5-mm-OD stainless-steel tube 
with a wall thickness of 0.4 mm. The tube is loaded 
with pellets of sintered plutonium dioxide with a 
density 88 to 93% of theoretical. In the highest burnup 
assembly, the initial gap between the pellets and 
cladding was 0.03 to 0.15 mm in 4 pins and 0.04 to 
0.07 mm in the remaining 15 pins. The height of the 
fuel pin is 280 mm. In the upper part of the fuel pin, 
the plenum had a height of 1 to 4mm to compensate 
only for axial thermal expansion of the pellets. There 
was little provision, therefore, for accommodation of 
the internal pressure arising from fission-gas release, 
and thus the failure of the fuel pins is believed to arise 
from the excessive pressure stress on the cladding. All 
the fuel pins are wire wrapped and assembled into a 
hexagonal tube of stainless steel with a wall thickness 
of 0.5 mm and a diameter of 26 mm. 

No external changes were observed on inspection 
of the highest burnup assembly. After cutting of the 
hexagonal tube, it was established visually that the fuel 
pins had a bright surface, and the spacing wires were 
unchanged. Only one pin in the assembly was found to 
be damaged—the central one. The damage occurred 
along the entire length of the pin, and cracks in the 
cladding were distributed along opposite faces. 

Fission-gas release from the fuel was 80%, and, 
because of the small plenum, it was established by 
calculation that under operating conditions the pres- 
sure inside the fuel element was 200 to 250 atm, which 
corresponds to a very high stress on the cladding. 
Under these conditions of high internal stress, it is 
remarkable that there were only a relatively small 
number of failures. 

The fuel pellets in undamaged pins mainly showed 
radial cracks, and there was no gap between the fuel 
and the cladding. Pellets in the central part of the 
damaged pin had considerable ovality, and extrusion of 
the fuel into the cracks in the cladding had occurred, 
confirming that plastic flow of the plutonium dioxide 
had occurred during irradiation. 

Microstructural analysis shows that the PuOQ, pel- 
lets from the center of the undamaged pins contained a 
large number of pores whose shape, size, and distribu- 
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tion varied according to the transition from the 
periphery to the center of the pellet. In the peripheral 
layer of the pellet, minute pores were distributed inside 
the grain, with larger pores along the grain boundary. 
Near the center of the pellet, the pore size inside the 
grain was increased; the number of pores was reduced 
but became larger around the boundary of the grain. 
At certain places in the center part of the fuel element, 
the pores around the grain boundaries were joined in a 
continuous lattice. Quantitative metallographic analysis 
showed that in the inner part of the pellet, with a 
diameter of 2 mm, the pores constitute 20 to 25%, and 
the porosity of the peripheral layer is 10%. The average 
porosity throughout the whole section of the pellet 
was 14%. In a PuO, pellet from the damaged pin, at a 
depth of 0.5 mm from the surface, large pores were 
observed with a diameter up to 50 um. 

If, as assumed, the oxide fuel swells to the extent 
of 1.7% per 1 at.% burnup due to the formation of 
solid fission products, then at a burnup of 6% the 
increase of volume should be ~10%. The initial 
porosity of the fuel is ~10%, and for the condition of 
compensation for the increase of volume due to 
porosity, the diameter of the fuel pellet should not be 
increased. However, according to the data the average 
porosity of the pellet was not reduced, but increased 
by 2 to 4%, which indicates swelling of the fuel due to 
the formation of fission-gas bubbles. Swelling of the 
fuel in the damaged pin was 2.8% per | at.% burnup of 
heavy atoms, i.e., it turned out to be considerably 
greater than in the undamaged pin because of the 
unrestrained swelling in the failed region. 

On the basis of the results of this investigation of 
the highest burnup assembly, the following conclusions 
can be drawn: 

1. The fuel pins of an assembly with stainless-steel 
cladding and fuel in the form of PuO, pellets, after 
operating at burnups up to 6.1 at.%, were found to be 
in a satisfactory state, except for the central pin, which 
had two long cracks on opposite faces of the cladding 
along the entire length of the fuel element. 

2. The quantity of gaseous fission gases released is 
about 80% of the amount formed, and the high 
internal gas pressure resulting from the small plenum 
volume was responsible for the failure. 

3. Fission-gas swelling during irradiation of PuO, at 
a temperature of 1650 to 2550°F is observed. The 
maximum increase of volume of the fuel in the 
damaged pin was about 2.8% per at.% burnup of heavy 
atoms. 

4. A tentative estimate of the stressed state of 
undamaged pins and the increase of their diameters 
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(~2 to 3%) confirms that these pins were close to 
failure. 


Failed-fuel propagation had not occurred, and this 
experimental experience, although obtained on PuOQ, 
fuel as opposed to mixed-oxide fuel, indicates fuel-pin- 
failure propagation is unlikely with fast breeder fuel. In 
spite of the severe failure in the one pin, there was no 
fuel—coolant reaction, and the sealing effect of the 
plastic PuO, resulted in no entry of fuel into the 
coolant stream. 


Total Irradiation Experience 


The total experience to date from EBR-II, DFR, 
Rapsodie, and BR-5 is summarized? in Table 8. The 


Table 8 Oxide Irradiations in Fast Test Reactors 








No. of pins No. of pins 
Country Reactor irradiated failed 
USSR BR-5 4600 120-150 
United Kingdom DFR 800 50 
France Rapsodie 8000 l 
United States EBR-II 1400 10 





failure rate has been about 1.4%. If one excludes such 
anomalies as the early BR-5 failures, which were caused 
by too small a fission-gas plenum, and eliminates the 
DFR failures, which arose from cover-gas entrainment 
in the coolant, the failure incidence becomes very small 
(<0.1%). This low incidence shows that, despite the 
wide range of parametric variations, fairly rudimentary 
models of mixed-oxide fuel behavior permit successful 
design of fuel pins. Thus, despite the large increase in 
the number of fuel pins subjected to testing, the 
number of failures has steadily decreased each year. 
The mixed oxide/stainless steel system has there- 
fore been shown to be a reliable one for the LMFBR, 
although more detailed work is needed to define the 
performance limits under the high-fluence conditions, 
higher temperatures, and the correct burnup-to-fluence 
ratio for the LMFBR. Many features may contribute to 
the excellent endurance shown by the test results, such 


as the plasticity of the fuel and the accommodation of 


fuel swelling by the incorporation of initial porosity; in 
addition, there is the strong indication of the advantage 
conferred by stainless-steel swelling in delaying the 
point at which fuel—cladding interaction occurs to very 
high burnup. Creep under irradiation can also play an 
important role. The detailed interplay of the main 
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phenomena involved is discussed later in the sections — 
on analytical models of fuel-pin and fuel-assembly 
behavior. 


EFFECTS OF FAST-FLUX 
IRRADIATION ON CLADDING 
AND STRUCTURAL MATERIALS 


Stainless-Steel Swelling 


Cawthorne and Fulton?? first reported that voids 
or empty cavities were formed in stainless steels 
irradiated in the DFR to doses in excess of 10? 
neutrons/cm?. During development work on fuel ele- 
ments for fast reactor applications, electron-micro- 
scope examination by the thin-foil technique was 
carried out on samples of stainless steel irradiated in 
the DFR, in the form of cladding on experimental fuel 
elements. The steel had a composition falling within 
the AISI 316 stainless-steel specification. A common 
feature of all the specimens examined after irradiation 
to total neutron doses >10?? neutrons/cm? in the 
temperature range 750 to 1130°F was the presence of 
spherical or polyhedral cavities ranging in size from the 
smallest observable to 1500 A. 


Bubbles had been seen before in irradiated stainless 
steel,?* and their origin had been ascribed to the 
segregation of helium generated by (n,a) reactions 
from !°B when the flux is chiefly thermal, but also 
from nuclear reactions from the constituents of stain- 
less steel in fast fluxes. In the fast reactor irradiations, 
however, the void population was too great to be due 
to helium-gas bubbles. 


Interactions with fast neutrons result in the dis- 
placement of atoms from their equilibrium lattice sites, 
creating “interstitials” and ‘“‘vacancies.”’ Interstitials are 
displaced atoms that are forced into nonequilibrium 
sites within the solid, and vacancies are the atomic 
holes that remain within the solid after the atoms are 
displaced. At temperatures relevant to the economical 
operation of fast reactors, both of these defects can 
migrate throughout the solid as a result of thermal 
activation. The interstitial migrates by jumping from 
one site to another, whereas the vacancy migrates as a 
result of a neighboring atom jumping into the vacant 
hole. At any point, large numbers of both interstitials 
and vacancies are freely migrating at random through- 
out the irradiated components. Such migrating defects 
can suffer a variety of fates: an interstitial and a 
vacancy can annihilate one another by recombining; 
both interstitials and vacancies can be captured and 
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lost at extended defects within the solid, such as 
dislocations or grain boundaries; interstitials can cluster 
together to form plates of extra atoms within the solid; 
and, finally, vacancies can cluster together to form 
either plates of missing atoms or large holes, such as 
voids, Plates of extra or missing atoms give rise to 
circumferential strains?* within the solid, which are 
observed as rings or loops. Therefore plates of extra 
atoms are called interstitial loops, and plates of missing 
atoms are called vacancy loops. 

When steels or nickel-based alloys are irradiated at 
room temperature, the lattice expansion due to inter- 
stitial loops and the contraction due to vacancy loops 
essentially counterbalance one another. However, 
during irradiation at elevated temperatures, e.g., 
930°F, a substantial fraction of vacancies agglomerate 
to form three-dimensional cavities rather than loops. 
Such cavities produce vanishingly small strains within 
the solid which can no longer balance the expansion 
produced by interstitial loops. Since agglomerating 
interstitials can only form loops, the result is a net 
expansion of the solid. 

Current theories of void formation suggest that 
trace quantities of gaseous elements, either occurring 
naturally within the solid or created during irradiation 
by transmutation, are necessary to nucleate three- 
dimensional voids. One such source is helium atoms, 
which are continually created within the reactor 
components as alpha particles via nuclear reactions. 
Helium atoms are highly insoluble within solids and 
rapidly collect together to form minute pockets of gas. 
It is these pockets of gas which are thought to act as 
perfect nuclei around which the migrating vacancies 
grow into voids. 

The voids therefore develop through the condensa- 
tion of vacancies, which are generated prolifically 
during irradiation in a fast-neutron flux, on the helium 
nuclei produced by (n,a) reactions. This conclusion 
was supported by the behavior of the voids on 
annealing in the Cawthorne—Fulton experiments. 
Specimens containing voids of average diameter 210 A 
and occupying 1.5% of the total volume were annealed 
for 1 hr at 1290 and 1650°F. Although the 1290°F 
treatment left the voids relatively unchanged, they 
were almost completely removed by the 1650°F 
treatment, leaving a residue of porosity confined 
mostly to the grain boundaries which could reasonably 
be equated with the helium content. 

In typical irradiation experiments, an axial tem- 
perature gradient is maintained in fuel pins that are 
exposed in flowing sodium; one such pin was used to 
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investigate the effect of irradiation temperature on 
void population. Cladding specimens were taken at 
intervals of 3in. along a 22-in. fuel pin clad in 
solution-treated 316 stainless steel after irradiation in 
the temperature range 520 to 1040°F. Voids were seen 
only in the three hottest specimens, the average void 
size increasing with irradiation temperature. The total 
void volume peaked at 950°F. 

Since Cawthome and Fulton’s initial results, fur- 
ther papers have been published** reporting observa- 
tions of voids in other metals as well as in various 
stainless steels. The volume changes associated with the 
presence of voids is of major technological significance, 
and thus it is of practical importance to achieve an 
understanding of the factors affecting both the forma- 
tion of voids and their growth. 


Theories of Swelling 


The two theories of irradiation-induced swelling 
which have received the most attention center on the 
homogeneous-nucleation-and-growth model proposed 
by Harkness and Li®® and the diffusion-controlled 
kinetic model of pore growth proposed by Bullough, 
Eyre, and Perrin.? 7 

The Harkness—Li model has been employed with 
some success to describe the swelling profile of two 
flats of a safety thimble irradiated in the EBR-II to a 
calculated fluence of 1 x 107° neutrons/cm?. The 
model of Bullough et al. has successfully predicted void 
growth under irradiation and its dependence on such 
parameters as the damage production rate and tempera- 
ture.3® 

Both theories require the solution of coupled 
differential equations that describe the steady-state 
conditions of a defect production rate balanced by the 
rate of loss of point defects through recombination and 
migration to microstructural “sinks,” such as voids, 
dislocations, grain boundaries, and precipitate particles 
under various boundary conditions. Since the micro- 
structure changes continually during irradiation, it is 
necessary to solve the conservation equations in each 
time increment by numerical procedures. 


Consideration is given to the vacancy and inter- 
stitial concentrations under steady-state conditions 
during irradiation. Account is taken of the loss of 
vacancies and interstitials by recombination and by loss 
to dislocations or other sinks. It is assumed that, in 
addition to random diffusion process, the interstitials 
undergo a preferential drift migration to dislocations 
due to the long-range interaction between the stress 
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field of the dislocation and the misfit strength of the 
interstitials; there is considerable evidence that this 
interaction is negligible in the case of a vacancy. Thus, 
by taking account of this preferred drift of the 
interstitials to the dislocations and the effects of 
recombination, the condition necessary for void 
growth is deduced. 

Both theories assume that one defect sink, i.e., the 
dislocation, has a sufficiently large preferential attrac- 


tion for interstitial atoms to permit the existence of 


the excess concentration of vacancies required for 
growth of the voids. The complex dislocation arrange- 
ment that then influences the kinetics of void forma- 
tion is approximated in both cases by a regular 
geometric array or a continuous distribution of point- 
defect sinks. One major difference is that the Harkness 
and Li model recognizes the existence of other sinks, 
such as grain boundaries and precipitate particles, 
whereas the Bullough, Eyre, and Perrin model assumes 
the growth process to be governed by the effective 
dislocation-sink density alone. The validity of this 
assumption must be tested, for it implies that grain 
refinement and precipitation reactions will have little 
influence on void formation. This conflicts somewhat 
with the observations that alloys with finely dispersed 
coherent precipitates (i.e., Nimonic PE16) swell less 
than similar alloys and that incoherent precipitates in 
various alloys are frequently associated with very large 
voids.” 3 

The only metallurgical processing variable that 
could have a significant influence is the amount of cold 
work introduced into the material, which determines 
the dislocation density. It is known experimentally, as 
will be discussed later, that 20% cold-worked 316 
stainless steel shows considerably less swelling than 
solution-treated material. The Bullough model predicts 
a linear variation of volume increase with dose for 
dislocation densities representative of solution-treated 
material but a dependence of AV/V on the square of 
the dose for cold-worked material. Thus, although the 
cold-worked material may show considerably less 
swelling initially, the rapid increase of swelling with 
increasing dose may ultimately result in larger volume 
increases than would have occurred with solution- 
treated material. In contrast, the Harkness and Li 
model predicts that swelling is actually enhanced by 
increasing the dislocation density at low irradiation 
temperatures. At high temperatures, e.g., T > 930°F, 
the addition of sinks through cold work, grain refine- 
ment, and precipitate dispersion reduces the swelling 
by decreasing the vacancy supersaturation to very low 
levels. 
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Although both approaches explain qualitative fea- 
tures observed experimentally, particularly the varia- 
tion with temperature and the influence of cold work, 
further theoretical development is needed since: 

1. Neither theory takes into account the possible 
effects of heterogeneous and continuous void nuclea- 
tion, as suggested by some recent experimental obser- 
vations.? ? 

2. Saturation of the swelling produced by irradia- 
tion with increasing time is intuitively expected. As 
currently formulated, neither model contains explicit 
upper bounds on saturation, although both groups have 
recently proposed saturation mechanisms.” 4°44! 

3. There is some evidence, as will be detailed later, 
that stress-induced void growth may govern swelling of 
fuel cladding at high irradiation temperatures. At 
present, neither theory adequately accounts for the 
presence of a hydrostatic stress during irradiation. 


Analysis of Stainless-Steel 
Swelling Data 


Simultaneously with the development of theories 
and models for the stainless-steel swelling phenomenon 
in fast fluxes, analysis of experimental swelling data as 
they became available from fast test reactors was 
carried out as a priority task in many laboratories. The 
main aims were: (1) to determine the relations with 
fluence and temperature so as to be able to predict the 
swelling at the higher fluences appropriate to the FFTF 
and the LMFBR and (2) to start to develop design 
equations incorporating these parameters as rapidly as 
possible. Early attempts*? to fit the data to equations 
of the form 


SE = A(gtyB exp (— C/RT) (1) 
where AV/V is swelling, @t is fluence (EF > 0.1 MeV), 
and A, B, and C are constants, were successful in that 
these equations (1) gave the first indication, prior to 
theoretical studies, of the variation of swelling with 
temperature and (2) showed that the swelling—fluence 
relation appeared to follow a power law. 

In view of the uncertainties in flux and tempera- 
ture in the data available, it is not unnatural that 
estimates of the constants A, B, and C by different 
laboratories vary. For example, Boltax and Claudson*? 
are providing a continuing review of stainless-steel 
swelling data in support of FFTF core design studies. 
The swelling equations recommended in August 1970 
are as follows: 
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For solution-treated 316 stainless steel: 


AV 


y = (gt) 2-05-27 /0+78/07 ee 40)107?°] 


% 





exp (32.6 — 5100/T—0.015T) (2) 
For 20% cold-worked 316 stainless steel: 


AT =9.0 + 1035 (or)! > [4.028 — 3.712 x 10? 


x (T — 273) + 1.0145 x 10% (T — 273)? 
=—7F59 10° G79 | (3) 


where %(AV/V) = stainless-steel swelling, vol. % 
gt =fast fluence (£>0.1 MeV), neu- 


trons/cm? 
T = absolute temperature, °K 
6 =T — 623 


Figures 7 and 8 illustrate the temperature and 
fluence dependency of swelling from these equations. 
Equations 2 and 3 are empirical, and there is no 
particular significance to the forms of the equations. 
The temperature-dependent fluence exponent indi- 
cated in Eq.2 is equal to 1.54 and 1.82 at 400°C 
(750°F) and 500°C (930°F), respectively. For reactor 
design studies, these equations are usually used with 
limits of +60 and —35% of the nominal values. In the 
analysis of high burnup EBR-II fuel pins, which will be 
discussed later, the swelling equation used for solu- 
tion-treated 347 stainless steel is 20% higher than that 
used for 316 stainless steel. This position on 347 
stainless steel is based on limited immersion density 
data available at a fast fluence near 4 x 107? neu- 
trons/cm?. 
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Fig. 7 Effect of temperature on the swelling of 316 stainless 
steel, 
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Fig. 8 Effect of fluence on the swelling of 316 stainless steel. 


Similarly, early empirical expressions developed by 
Oldberg et al.** (utilizing analytical techniques for 
nonlinear equations) described the swelling behavior of 
304 and 316 stainless steels as a function of fluence 
and temperature of irradiation. As an example, when 
the change in volume was considered as the dependent 
variable, the expression that fitted the available data 
was found to be 


= = 7.53 x 10°(F)!-%? 


x (T — 360)°- 95 (775 — T)®-7!! = (4) 
with a correlation coefficient of 0.938. When the 


fluence was considered the dependent variable, the 
expression that satisfied experimental findings was 


| cb 
S 


—=3.14x 10%(F)!*8! 


— 


x (T — 360)°-©8? (775 — T)9-577 = (5) 


with a correlation coefficient of 0.845, where AV/V is 
the change in volume in percent, F is the fluence in 
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107? neutrons/em? (E>0.1 MeV), and T is the 
temperature of irradiation in degrees centigrade. Evalu- 
ation of this equation beyond the limits of the input 
data was considered a first-order estimate at best. 

When the Boltax and Claudson ground rules are 
applied [(1) use only those data from specimens 
irradiated within +10 in. of the core midplane, and (2) 
adjust each data point by subtracting 0.06% from the 
measured swelling value], the value of B determined 
from a least-squares fit of the solution-treated 304 and 
316 stainless-steel data is ~1.7 with temperature and 
fluence as the independent variables. Bloom and 
Stiegler?® have examined these data using somewhat 
different ground rules, as summarized in Table 9. An 
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and 


(47 z 0.06)  (61)!-°5 (7) 
At 460 to 470°C, a least-squares fit gives 


« (¢t)'-°* and (Ar 0.06) « (gt)':?4 (8) 


This analysis suggests that the fluence dependency of 
the swelling is not constant but is at least a function of 
irradiation temperature (Eq. 2 also exhibited this fea- 
ture). This is substantiated by electron-microscopy 


Table 9 Analysis of Swelling Data from the EBR-II 





No.of — Expression Value of constants 











data used 
Ground rules points in fit* A B c 

All data 93 I exon” 0.78 1180 
85 II 21°%:107> 1.06 2120 
Only data from +10 in. of core 55 I 22 AO? 1.54 4680 
midplane 3D II DSKAG° Lz 5040 
Only data from +10 in. of core 17 ] 6.9x10'? 0.64 5080 
and only irradiation temperatures 17 Il 25x10"* 0:71 5470 


of 370 to 380°C and 460 to 470°C 





*1, AV/V = A(ot)8 exp (—C/RT). Il, AV/V @ — 0.06 = A(ot)8 exp (—C/RT). 


examination of the values of B listed in Table 9 shows 
that it is very sensitive to the data considered. Bloom 
and Stiegler believe that the best data available 
correspond to those specimens irradiated near the 
bottom or the top of the EBR-II core, where the 
temperature is known to be 370 to 380°C (700 to 
715°F) or 460 to 470°C (860 to 880°F), respectively, 
in addition to the criterion of +10 in. from the reactor 
midplane where the flux and neutron spectrum are best 
known. When these ground rules are applied, the value 
of B is found to be 0.64 to 0.71, depending on whether 
or not the data are adjusted by 0.06%. These values are 
in good agreement with the values obtained when all 
data are considered. Bloom and Stiegler have extended 
the analysis further and extracted and analyzed sepa- 
rately the data for solution-treated specimens irra- 
diated at 370 to 380°C and 460 to 470°C. At 370 to 
380°C, a least-squares fit gives 


AV 0.84 


Vv = (6) 
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observations of void concentrations and sizes as func- 
tions of irradiation temperature and fluence. Concen- 
trations at 370 to 380°C are always higher than those 
at 460 to 470°C, but, in the neighborhood of a fluence 
of 107? neutrons/cm?, the concentrations at 460 to 
470°C are increasing at a faster rate than at 370 to 
380°C. Thus it appears that, at these temperatures in 
stainless steel, voids grow very quickly to a limiting size 
above which their growth rate is very small. In this case 
the swelling rate approximately reflects the nucleation 
rate. 

All the above equations for 316 stainless steel were 
developed primarily using the more abundant data on 
304 stainless steel, which was considered a similar 
material, Recent data*® show that 304 and 316 
stainless steels have somewhat different swelling be- 
haviors, and individual equations are now being devel- 
oped by various laboratories. The solution-treated 304 
stainless-steel swelling appears to be about twice that 
of solution-treated 316 stainless steel for essentially 
identical irradiation conditions. In addition, recent 
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unstressed cold-worked 316. stainless-steel swelling 
data*® exhibit significantly lower density changes than 
predicted by Eq. 3, which was drawn from data taken 
from material known to be under stress. The difference 
between the two sets of data may reflect either the 
state of stress or material variability, or both. 

British researchers correlate their swelling data 
somewhat differently from their American  col- 
leagues.*7 They extract a linear relation after plotting 
their data with dose above a threshold value. This 
incorporates the observation that swelling in stainless 
steels does not appear until fluences equivalent to 
about six displacements per atom [1.3 x 10?? neu- 
trons/cm? (E > 0.1 MeV)] have been attained. 


These comparisons illustrate the difficulties of 


analysis and extrapolation of the present data to the 
high fluences needed for the LMFBR application. As a 
result, three approaches have been identified and 
implemented. 

1. Concentration on obtaining accurate data from 
the national program in EBR-II. 

2. Continuous evaluation of data as they arise from 
the test-reactor programs. 

3. Use of conservatism at this point in time for 
design purposes, as indicated in the discussion on the 
Boltax—Claudson analysis. 


Effects of Stress on Stainless-Steel Swelling 


Analysis of the irradiation behavior of stainless 
steel also indicates that stress can have a significant 
effect on the swelling of stainless steel over the 
temperature range of interest in fast reactors. Recent 
studies by Boltax et al.4® identified three possible 
mechanisms for stress effects on stainless-steel swelling, 
including: 

1. The growth of grain-boundary cavities or helium 
bubbles under hydrostatic tension. 

2. The hydrostatic component of irradiation creep. 

3. The growth of voids within grains under hydro- 
static tension. 


The growth of grain-boundary cavities or helium 
bubbles under hydrostatic tension has been reported 
by numerous investigators.?* Stress-rupture tests at 
temperatures of 700°C (1290°F) on stainless steel 
containing 40 ppM helium indicate that the major 
mode of deformation is due to the growth of grain- 
boundary cavities.4° This high-temperature deforma- 
tion mechanism is very important in determining 
density changes, creep rates, and strain limits in 
stainless-steel cladding stressed at temperatures above 
600°C (1110°F). 
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Theoretical studies*® of the effects of stress on 
voids within grains, continuing the work of Harrison*° 
and of Gulden and Kaae,°! indicate that, under high 
fast-flux irradiation [~10'* neutrons/(cm? )(sec) with 
E> 0.1 MeV], void growth can occur at ~0.57;, due 
to stress-motivated vacancy flow to voids. Hydrostatic 
tension produces steeper concentration gradients of 
vacancies around voids which provide an extra driving 
force for void growth. The equations that describe this 
phenomenon are given by: 








4 
AV =~ ar°N (9) 
dA V 2 _dr 4 3 aN 
it = 4rr Na t3™ d (10) 


where AV = volume change arising from void nuclea- 
tion and growth 
r = void radius 
N = number of voids 
t = time 


Experimental work by Straalsund and Bates>? 
indicated that, for stress-free swelling, r is relatively 
independent of fluence (or time) and that swelling is 
principally controlled by nucleation. If we assume that 
nucleation is not affected by stress, then Eq. 10 can be 
rewritten as 


dAV 2n7 dr 4 3 aNgr 
mal =é4 N é — eee "2 ae 
er a a as ay) 


and 














AV AV dr A 
er) ear 4 Or) 
; V eS . V seat ,* V SF 


(12) 
dt Rep dt 

where (AV/V)gsp is stress-free swelling, Ngp is stress- 

free number of voids, and Rg p is stress-free radius of 

voids. 


The equation for the rate of stress-enhanced void 


growth, derived from the work of Gulden and Kaae,°*? 
is given by 


dr_ DEQoy 


= 3 
dt RspkT 7) 


Thus the total swelling is given by the sum of 
stress-free swelling and stress-enhanced swelling as 
follows: 


REACTOR TECHNOLOGY, Vol. 15, No. 1, Spring 1972 








36 FAST BREEDER FUEL P. Murray 
V V , The application of these phenomena, together with 
‘dee 37? & DE Qoy MP 0 ios cine Sts = : 
a V /sp VAsr (14) irradiation creep and irradiation embrittlement, is 
"4 7 dt pe ae reviewed in later sections. 
dt R¢pkT 


where DE = effective vacancy diffusion coefficient 
Q = atomic volume 
oy, = hydrostatic tension 
k = Boltzmann’s constant 


Figure 9 illustrates the magnitude of stress-en- 
hanced swelling from this model in solution-treated 


Accelerator Experiments 
on Void Formation 


Prior to 1970 the only factor known to influence 
the swelling behavior of cladding and structural mate- 
rials was the incorporation of cold work in the 
structure, which, as shown previously, led to a signifi- 
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Fig. 9 Summary of irradiation phenomena that can produce volume changes in solution-treated 316 stainless steel. 


316 stainless steel at a fluence of 5 x 10?? neu- 
trons/em* (E> 0.1 MeV) and a constant hydrostatic 
stress of 10,000 psi. The effect of stress-enhanced 
swelling is small below 550°C and is significant 
between 600 and 700°C. For comparative purposes the 
effects of the volume change arising from irradiation 
creep are included in Fig.9 to show the relation 
between the irradiation effects that can produce 
volume changes in stainless steel. 

Recent theoretical investigations®* on the effects 
of stress on swelling in a triaxial state of stress 
supported qualitatively the conclusions drawn from the 
model outlined above. These investigations have also 
shown that stress may enhance swelling by two 
mechanisms: one mechanism is the enhanced diffusion 
of point defects under stress, and the other mechanism 
is the alteration of the efficiency of the various sinks 
for point defects. 
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cant reduction in swelling. More recently, investiga- 
tions®4°55 of nickel-based alloys containing a fine 
gamma-prime precipitate have led to further advances 
in establishing the directions by which swelling be- 
havior can be reduced by alteration and control of 
microstructure. One such gamma-prime alloy is 
Nimonic PE16, and the available information on its 
swelling characteristics is derived from ion-bombard- 
ment and fuel-pin tests. The data show that the 
swelling rate for this material is low and significantly 
less than 20% cold-worked 316 stainless steel. 

The ion-bombardment tests provide a quick assess- 
ment of the relative swelling behavior of various 
materials over a range of fluences, and equivalent 
LMFBR fluences of 3 x 107? neutrons/cm? can be 
achieved in a few days. A set of experiments has been 
carefully conducted on solution-treated stainless steel, 
20% cold-worked 316 stainless steel, and solution- 
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treated and aged Nimonic PE16, leading to a good 
comparison of candidate LMFBR materials. Before 
these results are discussed, the methods employed for 
these tests and their accuracy will be briefly outlined. 

The tests were conducted by bombarding thin foils 


of materials with 20-MeV C2* ions. For simulation of 


reactor conditions, particularly void-nucleation con- 
ditions, the foils were prebombarded with helium to 
yield ~10 at. ppM of helium in the samples. Volu- 
metric swelling was determined by electron microscopy 
involving the measurement of the mean void size and 
the void concentration. The latter involves measure- 
ment of the foil thickness by stereomicroscopy or 
step-trace measurement, and these techniques limit the 
overall accuracy to about +20%. The technique thus 
provides a reliable means of screening and comparing 
the swelling behavior of different materials and alloys. 

In Fig. 10 the swelling of 20% cold-worked and 
solution-treated 316 stainless steel is compared to that 
of Nimonic PE16 in the temperature range 500 to 
700°C (930 to 1290°F) at an exposure of 40 displace- 
ments/atom (dpa). This exposure is equivalent to 
~9 x 107? neutrons/cm? (EF > 0.1 MeV) in the FFTF, 
and these data can be used to establish the temperature 
dependency of swelling for Nimonic PE16. The results 
clearly demonstrate that swelling of Nimonic PE16 is 
substantially less than that of the stainless steels. Thus, 
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Fig. 10 Comparison of 316 stainless steel and Nimonic PE16 
swelling. 
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at 500°C, swelling of the Nimonic alloy is less than 
10% of that of cold-worked stainless steel, and, at 
600°C (1110°F), it is approximately 30% of that of 
cold-worked stainless steel. It is important to note that 
the high-flux levels used in ion-bombardment tests shift 
the temperature of peak swelling by approximately 
+125°C. Thus the peak for solution-treated 316 stain- 
less steel is found at 600°C in these tests rather than at 
450 to 475°C, as observed from test-reactor results. 
This effect has been predicted from theory and does 
not influence the relative comparison of alloys, but it 
must be recognized when considering applications of 
the data to design. 

The swelling of the same alloys is compared in 
Fig. 11 as a function of fluence at a constant tempera- 
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Fig. 11 Effect of fluence on the swelling of Nimonic PE16 
and 316 stainless steel at 525°C (977° F). 


ture of 525°C (977°F) and can be used to establish the 
fluence dependency of swelling. The conversion of 
fluence to total flux for the FFTF and the LMFBR is 
also shown. These results demonstrate the following 
key points: 

1.The gamma-prime Nimonic PE16 alloy swells 
substantially less than the stainless-steel alloys at doses 
>20 dpa. 
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2. At 100 dpa [2.3 x 1073 neutrons/cm? (E > 0.1 


MeV) for the LMFBR] and 525°C, swelling of 


Nimonic PE16 is <0.4%. 

3. Assuming the same dose dependence at 600°C 
(1110°F) as for 525°C, the peak swelling of this alloy 
is expected to be <3% at 600°C and 100 dpa. 

4, The fluence dependence of swelling for Nimonic 
PE16 is approximately ($f) *, compared to (¢r)' ** for 
stainless steels. Because of the different slopes, the 
curves for the stainless-steel and Nimonic alloys inter- 
sect at low doses. Thus test-reactor data in the range 
400 to 500°C obtained near or below the intersection 
points would show that stainless steels swell equal to or 
less than Nimonic PE16 and result in a misleading 
conclusion. Therefore, establishing the swelling advan- 
tages of this alloy requires in-reactor testing at fluences 
>5 x 1072. neutrons/em? (E> 0.1 MeV) for this 
temperature range. At higher temperatures, 600°C, the 
lower swelling of Nimonic PE16 would be observed at 
lower fluences. 


Nimonic PE16 contains equal amounts of titanium 
and aluminum at concentrations between 0.9 and 
1.5 wt.%. With a suitable heat treatment, precipitation 
of these elements occurs in the form of the gamma- 
prime phase, which has the composition Ni3(TiAl). 
The crystal structure of gamma prime is face-centered 
cubic (fcc), and its lattice parameter is only 0.03% 
greater than that of the austenite matrix. Thus the 
misfit between the spherical gamma-prime particle and 
the matrix is extremely small and results in a low- 
energy, fully coherent precipitate—matrix interface 
which, in turn, is responsible for the high thermal 
stability of the gamma-prime dispersion. Typically the 
gamma-prime particles are ~100A in diameter and 
spaced ~200 A apart, the concentration being in the 
range 10!° particles/cm?. 

Although the precise reason for the swelling resis- 
tance of this alloy is not fully understood at present, 
several possible explanations may be advanced based 
on the microstructure. These are: 

1.The small misfit associated with the gamma- 
prime particle locally compresses the austenite lattice. 
There is an elastic interaction between this strain field 
and vacancies which results in vacancies being tempo- 
rarily trapped, thus increasing the chances of their 
annihilation by interstitials. 

2. Vacancy trapping may occur inside the gamma- 
prime particle, i.e., there is a binding energy between 
vacancies and gamma-prime particles. Vacancies are 
thus delayed in their migration to voids and other 
sinks, and the rate of recombination is increased. 
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3. The austenite matrix contains titanium and 
aluminum in solution, both of which may be expected 
to interact with vacancies. The effect of such a 
solute—vacancy binding energy would be to reduce 
vacancy mobility and again promote recombination 
with interstitials. 

4. The high precipitate concentration ensures that 
interaction between the growing interstitial loops and 
gamma-prime particles occurs at low fluences. The 
resulting restriction in the ability of loops to expand 
reduces their preferential attraction for interstitials. 
This preferential attraction is generally considered the 
primary cause of swelling since it results in excess 
vacancies in the system. Thus, with this reduced 
attraction, recombination and loss of interstitials at 
nonpreferential sinks is enhanced, and swelling is 
reduced. 


Careful electron microscopy and correlation with 
microstructural changes might distinguish between 
these possibilities. Until such studies are performed, it 
is assumed that each of these mechanisms may be 
inhibiting swelling in Nimonic PE16 to some degree. 

In summary, the ion-bombardment results at flu- 
ences equivalent to and greater than expected in the 
LMFBR demonstrate that gamma-prime alloys, such as 
Nimonic PE16, swell substantially less than cold- 
worked 316 stainless steel or solution-treated 316 
stainless steel. The results also show that the fluence 
dependence of swelling is considerably less than that of 
stainless steels. The trends need to be confirmed by 
test-reactor data, and the implications and potential 
benefits of low-swelling alloys are discussed later in this 
review. 


Ductility of Fast-Flux-Irradiated Stainless Steel 


The major available data on the ductility of 
stainless steel irradiated in fast test reactors are from 
short-time tensile tests. Significant information on 
ductility limits has been obtained from pressurized- 
tube experiments®® and also from analysis of fuel-pin 
test data.°’ If one extrapolates the strain at failure in 
the experiments of Bell and coworkers to the low 
strain rates expected in a fuel pin [10° in./(in.)(hr)], 
the strain at failure is about 0.7% for cold-worked 316 
stainless steel. Analysis*’ of the helium-bonded 
mixed-carbide pins irradiated by United Nuclear 
Corporation in the EBR-II given in Table 10 also 
indicates a limiting strain in the cladding (arising from 
the swelling of the carbide fuel) of about 0.8 to 1% for 
solution-treated 316 stainless steel. These thermal- 
creep strains at failure were obtained by subtracting 
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Table 10 Results of Helium-Bonded Mixed-Carbide Pins Irradiated by 


United Nuclear Corporation in the EBR-II 





Pin No. 81 82 83 84 85 86 87 89B 
Pellet density, 7% 98 98 98 98 98 98 98 98 
Fuel—cladding diametral 

gap, mils 8 J 6 10 10 6 2 
Cladding type 316 316 Incoloy 800 316 316 Incoloy 800 Incoloy 800 316 
Cladding outside 

diameter, in. 0.306 0.306 0.302 0.306 0.306 0.302 0.302 0.306 
Cladding thickness, 

mils 24 24 22 24 24 22 22 24 
Burnup, 10*MWd/ 

metric ton 57 ny 57 73 73 73 40 40 
Peak fluence, 

107? neutrons/cm? 3.8 3.8 3.8 4.9 4.9 4.9 i 1.1 
Max. cladding OD 

temp., °C 595 595 600 590 590 590 625 570 
Heat rating, kW/ft 21.6 21.9 21.9 19.8 19.8 19.8 28.8 30.3 
Cladding condition Failed Failed Failed Failed Failed Intact Intact Failed 
Max. OD increase, 

mils ns. 4.4 10.4 7.0 B 5 3.2 BS 
AD/D, % (total) 1.74 1.44 oa ie ee 1.14 


3.45 








the cladding swelling and irradiation-creep contribu- 
tions from the total diametral change. Slightly higher 
strains to failure (1 to 2%) have been obtained in the 
General Electric F2-series mixed-oxide pins irradiated 
in the EBR-II. 

The available information therefore indicates that, 
at the low strain rates corresponding to fuel—cladding 
interaction, the thermal-creep strain at failure is in the 
range of 1%. It should be noted that the present data 
are limited to fluences of ~7 x 10?? neutrons/cm?. 


Irradiation Creep of Stainless Steel 


Irradiation-creep data on cold-worked and solu- 
from 
fast-spectrum helical springs irradiated in the DFR at 


tion-treated 316 stainless steel are available? * 


250°C and uniaxial tensile specimens irradiated in 
thermal spectra.5 °°? Damage-function analysis of the 
thermal- and fast-spectrum data indicated that ir- 
radiation-creep strain is the superposition of two 
components: a hydrostatic component that represents 
a volume increase due to the application of a stress and 
a deviatoric or shear component.°***? The equations 
for the shear and hydrostatic portions of irradiation 
creep are given separately below. 
For solution-treated 316 stainless steel: 


<=0.18x 10° af] — exp (for/1.47 x 10?! )] 


ol 9! 


+ 1.87 x 103° fot (15) 


SH =2.7 x 10% afl — exp (—fét/1.47 x 10?')] 


+28.1x 10°° for (16) 


For 20% cold-worked 316 stainless steel: 


|| oI 


= (0.78 x 10° afl — exp (— fot/6.37 x 107°)] 


~ 
_ 


+ 1.86x 10°° for (17) 


=H = 11.73 x 10° afl — exp (— for/6.37 x 10?°)] 
H 


+279 x 103° for (18) 
where €/@ = ratio of equivalent strain to equivalent 
stress 
€1,/0,; = ratio of hydrostatic strain to hydrostatic 
stress 
ot = total fluence 
f =a neutron-energy-spectrum correction 
factor which is 1.0 for the FFTF mixed- 
oxide core and 1.75 for the EBR-II for 
rows 2 to 4 at the core center line 
a = exp (1.405 — 0.00277) 
Recently additional fast-spectrum data for solu- 
tion-treated 304L stainless-steel helium pressurized 


tubes in the EBR-II have been reported,°°® and this 
permitted the modeling of fast-spectrum irradiation- 
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creep data without the use of damage-function analysis 
to incorporate thermal-reactor data.°' Within the 
range of the data (250 to 380°C), the creep rate was 
athermal, and there was no hydrostatic creep compo- 
nent. Furthermore, the creep rate was linear with 
stress, in good agreement with previous results. The 
form of the equation is 


€=Aokot BoE¢t 
where F is the mean neutron energy, which is predicted 
to be 0.43 MeV for the FFTF core,°*® and A and B are 


materials constants: 


For solution-treated 316 stainless steel, 
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Corrosion and Mass Transfer 


The corrosion rates have been measured in loop 
systems by a number of investigators®?*®* with respect 
to the effects of temperature, oxygen content of the 
sodium, and the sodium velocity, and correlations have 
been derived relating corrosion rate to these factors.° * 
Typical results are given in Table 11, and the impor- 
tance of temperature as the main factor influencing the 
corrosion rate is well established. 

The compositional changes of the sodium-exposed 
surface have also been investigated in detail®® and 
metallic mass-transfer rates established. The average 
metallic mass-transfer rates observed so far are gener- 
ally acceptable for the LMFBR fuel requirements. 

Interstitial elements, principally carbon and nitro- 
gen, are among the most mobile elements in high- 
temperature sodium systems, and relatively minor 


0.99 x 1078 oo variations in the concentratio 1ese elements have a 
Aa ‘ if 6 >22,000 psi ariations in t 1e concel tration of tl S 
oO pronounced effect on the mechanical properties of the 
alloys used in the system. A typical carbon-gradient 
B=3.9x 107° 


For 20% cold-worked 316 stainless steel, 
A=1.7x 1073 


B=2.7x 107° 


EFFECTS IN SODIUM 


Insofar as the fuel-assembly components (cladding 
and structural materials) are concerned, the effects in 
sodium are mainly those of corrosion and mass 
transfer. In addition, the reactions between sodium and 
fuel following fuel-pin failure are being investigated 
extensively .°? 


analysis for 316 stainless steel exposed to 1100°F 
sodium for 2600hr is shown in Fig. 12. In_ this 
example, the surface carbon concentration has de- 
creased from ~600 to 400 ppM, and the depth of 
decarburization is about 10 mils. Comparative re- 
sults°” on different cladding materials are given in 
Table 12, which shows that the stabilized stainless-steel 
grades, such as 347 stainless steel, do not lose carbon. 

The effects on stress-rupture properties®’ follow- 
ing loss of carbon in 316 stainless steel are shown in 
Fig. 13. Therefore 316 stainless steel shows a loss of 
strength under interstitial transfer conditions, and, for 
this reason, the mechanical properties of low-carbon 
316 stainless steel are those taken for design applica- 
tions. The stabilized grades, 321 and 347 stainless 
steels, do not lose carbon; in fact, there are indications 


Table 11 Corrosion and Deposition for Primary Systems 


Steady-State Mass Transfer Rates After 2500- to 3000-hr Exposure 





Hot-leg material 316 S.S. 316 S.S. 316 S.S. 3168.8. 316S.S. 
Cold-leg material 304 S.S. 304 S.S. 304 S.S. 304S.S. 304S.S. 
Hot-leg temp., °F 1100 1100 1325 1100 1325 
AT. oF 300 300 300 300 300 
Oxide content, ppM <10 <10 <10 20-25 20-25 
Hot-leg velocity, ft/sec i 20 20 20 i 
Corrosion rate, 

mils/year —0.08 —0.102 —0.80 —0.39 -1.1 
Cold-leg velocity, ft/sec 7 15 15 15 7 
Deposition rate, 

mils/year 0.03 0.0427 0.071 0.145 0.184 
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Fig. 12 Carbon gradient in 316 stainless steel exposed to 
1100°F sodium for 2600 hr. 


Table 12 Bulk Carbon Contents After 4840-hr 
Exposure to 1300°F Flowing Sodium 





Carbon content, % 




















After Change, 
Material Initial exposure %e 
304 S.S. 0.062 9.017 —73 
316 S.S. 0.073 0.038 —48 
347 S.S. 0.051 0.052 +2 
Incoloy 800 0.088 0.098 +11 
100 

Velocity, 17 ft/sec 

Oxygen, <10 ppM 
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Fig. 13. Biaxial stress-rupture properties of 316 stainless steel 
in sodium and argon at 1300°F. 
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that carbon pickup may take place. The decarburiza- 
tion—carburization mechanisms that can occur with 
cladding and structural materials and their effects on 
mechanical properties are therefore being investigated 
extensively. In addition, the interplay between the 
sodium environmental effects on microstructure and 
the irradiation-induced effects, such as swelling and 
creep, are being investigated in the EBR-II program. 


Fuel—Sodium Reactions 


Following the failure of five encapsulated fuel pins 
in the EBR-II which on postirradiation examination 
showed high swellings ranging up to 34% maximum 
AD/D, the reaction between static sodium and mixed- 
oxide fuel has been investigated extensively.°? The 
reaction has been confirmed in both in-pile and 
out-of-pile static experiments, and the compound 
Na3(U,Pu)O, has been identified as the generic for- 
mula for the reaction product. This compound forms 
in excess sodium under conditions for which the 
oxygen activity is sufficient to stabilize the compound. 
Oxygen is the limiting species and may come from the 
sodium phase and/or from the fuel phase. The thresh- 
old oxygen levels in the sodium phase for stabilization 
of the compound are estimated to be in the range 10 to 
100 ppM. The threshold O/M in the fuel phase adjacent 
to the reaction compound is about 1.97 under con- 
ditions where the reaction product is at 1000°C 
(1830°F) and the sodium is at 900°C (1650°F). The 
threshold average O/M of the fuel, however, will be 
somewhat below 1.97 when the oxygen redistribution 
along the temperature gradient is taken into account. 
The threshold oxygen level for stabilization of the 
compound depends on the sodium temperature, the 
reaction product temperature, and the temperature 
gradient in the fuel phase. 

Swelling due to the compound formation arises 
from the entry of sodium into the fuel matrix along 
grain boundaries, which expands the fuel and the 
cladding in the vicinity of the defect. The swelling 
increases with increasing oxygen available in the fuel. 
This extra oxygen arises from redistribution of oxygen 
because of the temperature gradient and from that 
released by burnup of the fuel. The degree of diametral 
change appears to be larger than expected due to 
uniform swelling and is attributed to localized reaction 
near the point of failure. 

The rate of diametral swelling follows a 1/4-power 
time dependence which suggests a blocking mechanism 
wherein a diffusion-controlled process through a reac- 
tion product layer is rate limiting. Swelling occurs 
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rapidly at first and then slows down. The detailed 
results to date are in static sodium, and oxygen 
removal in flowing sodium can limit the swelling to 
values much lower than those observed in these 
experiments. The observations on failed fuel pins in 
flowing sodium (such as the BR-S failures discussed 
earlier, in which the reaction was not observed with 


PuO, fuel) are evidence for the dominant effect of 


oxygen removal under flowing conditions. Further 
work is needed to confirm this in mixed-oxide fuel- 
pin-failure tests under flowing-sodium conditions. 


ANALYTICAL METHODS TO DESIGN 
AND PREDICT PERFORMANCE 
OF FAST REACTOR FUEL PINS 


In the present state of fuel development, consider- 
able reliance is placed on irradiation tests and their 
interpretation. More information is required on fuel 
pins, subassemblies, and materials operating at the 
conditions expected for LMFBR designs. Progress in 
screening various fuel designs and establishing the 
effects of power ratings, temperatures, and burnup is 
limited owing to: (1) lack of high fast-flux test 
environments, (2) the length of irradiation time re- 
quired to obtain significant results, and (3) limited 
usefulness of thermal-flux tests. 

In parallel with the design screening and the 
irradiation experiments, an analytical approach is 
required to establish guidelines to assist in the defini- 
tion of irradiation programs and for the development 
of optimum fuel designs. When the fuel pin and 
assembly are treated as an interdisciplinary system to 
describe fuel performance during normal and off- 
normal operation, the following benefits result: 

1. Guidance and direction for irradiation testing. 

2. Extrapolation and interpolation from known to 
unknown regions of performance. 

3. Optimization of fuel designs. 


The recognition of these advantages had led to 
concentration on the development of fuel-performance 
codes over the past 3 years, particularly in the U.S. 
program.® § 

In more detail, the fuel-element-performance 
models are being used in several ways. The models 
force the materials experts to make quantitative 
predictions on fuel and cladding performance. This 
indicates areas where experimental and theoretical 
work are most needed. Comparison of model predic- 
tions with each new batch of fuel-element-performance 
data may indicate inadequacies in the medels used but 
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also indicates the degree of agreement in the perfor- 
mance of different fuel elements. Correlating the large 
quantity of performance data on fuel from the EBR-II 
program is the most important near-term benefit. As 
the models develop and experience is gained, they are 
being used to determine fuel designs and to predict the 
behavior of elements in the LMFBR. The most 
challenging objective is the use of the codes to predict 
the effect of various transients and off-normal opera- 
tion on fuel behavior. 

The models developed to describe the operation of 
fuel elements under normal conditions can be arbitrar- 
ily divided into two categories. In the first, the models 
are more phenomenological and consist of fast-running, 
largely steady-state models that use empirical relations 
to describe fuel behavior; for example, the relation 
between gas release and fuel swelling vs. burnup and 
temperature. Only the more observable characteristics 
of the oxide core and the cladding are modeled. This 
type of model draws heavily on the physical aspects of 
the fuel obtained from macroscopic measurements and 
metallography. Since the analytical models are consis- 
tent with the empirical observations, the codes are 
valuable for interpolation and limited extrapolation. 
OLYMPUS (Westinghouse)*® and SWELL (Argonne 
National Laboratory)°® are examples of these models. 
The work at General Electric Company,’7° Atomics 
International,” and Hanford Engineering Develop- 
ment Laboratory’? also belongs to this category. Many 
of these codes have detailed subroutines, and there is 
an overlap with the second category in which more 
realistic models of the various microscopic processes 
that are believed to control the behavior of the fuel 
element are formulated. The main objectives of these 
detailed codes are to relate materials properties and 
physical mechanisms to the widest possible range of 
fuel-element geometries, materials parameters, and 
operating conditions. Because of the firmer physical 
basis for the models, extrapolation over a wider range 
of conditions should eventually be valid. Also, they can 
be used to guide studies of fuel and cladding properties 
since these data are needed as input. CYGRO,7? 
LIFE,”* and FMODEL’® are examples of these codes. 
In conjunction with their development, considerable 
effort has been given to the subroutines of the code 
that analyze a particular mechanism of interest. Fuel- 
swelling and gas-release subroutines, for example, are 
areas being pursued actively by several investigators.” ° 

Analysis of the off-normal behavior of mixed-oxide 
LMFBR fuel is another area of active study since 
knowledge of fuel performance under such conditions 
is necessary for setting realistic design and operating 
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limits and establishing margins for unacceptable dam- 
age or failure. 

Characterization of failure models, pressures pro- 
duced by failure, and movements of fuel and coolant 
are needed for accident analyses.77 Experimental 
guidance is provided by the program in TREAT’” and 
other facilities for simulating LMFBR power excur- 
sions, loss of coolant flow, and short-term effects of 
local abnormalities (such as coolant-channel block- 
ages). 

Relatively long-term effects, such as fission-product 
migration or release, fuel restructuring, and degrada- 
tion of cladding mechanical properties, enter the 
transient analyses only as initial conditions. However, 
related phenomena arise, including, respectively, tran- 
sient pressures produced by fission-gas release upon 
transient heating, internal movement of fuel resulting 
from transient thermal expansion and melting, and 
cladding response to transient loading. 

The central problem in modeling fuel-element 
performance is the determination of the restructuring 
and shape changes of the fuel and cladding. The 
treatment of the shape changes of the cladding is 
relatively straightforward. The behavior of the mixed- 
oxide fuel is more difficult to describe in detail. As a 
result, the main difference between the various codes is 
the degree to which they model the changes taking 
place in the fuel. The more-advanced codes have 
subroutines that treat each of the following processes. 


Cladding Swelling. In mixed-oxide fuel elements 
with low smeared density (~85% of the theoretical 
density) and moderate burnup (~50,000 MWd/metric 
ton), little internal pressure is exerted on the stainless- 
steel cladding, and its dimensional changes are deter- 
mined by fast-neutron-induced swelling. The empirical 
correlation most often used in the United States for 
316 stainless steel has been given previously. 


Cladding Creep. The creep of fuel-element cladding 
is the sum of two parts: thermally activated creep and 
radiation-induced creep. Thermally activated creep 
data are taken from out-of-pile tests. The radiation- 
induced component is superimposed on this and found 
to be dominant at lower stress and temperature,’® as 
discussed previously. 


Cladding Corrosion. Corrosion affects cladding in 
at least three different ways: reduction in thickness 
due to sodium cladding attack;®? loss of strength due 
to removal of interstitial elements, such as carbon;*” 
or fission product attack of the cladding.'® However, 
the empirical data on the kinetics of these reactions are 
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quite limited at the higher sodium coolant tempera- 
tures where the attack is more severe. The data on 
fuel—cladding reaction have been summarized earlier, 
and attention has been drawn to the importance of the 
O/M ratio of the fuel and to the urgent need for more 
data in this area. For conservative design, allowances 
are made in terms of decreased effective cladding 
thickness. As reliable kinetic equations are developed 
for these reactions, they can be incorporated into the 
codes. 


Fuel. The low thermal conductivity of the oxide 
causes steep temperature gradients in the fuel. Conse- 
quently there are great differences in many of the 
properties of the fuel as a function of radial position. 
The changes that occur in the fuel are primarily 
changes in dimension due to swelling, creep, hot 
pressing, thermal expansion, and fuel cracking. The 
rates of all these kinetic processes change greatly with 
temperature. 


Swelling and Gas Release. In the phenomenological 
codes the fuel swelling and gas release due to gaseous 
fission products are calculated using empirical rela- 
tions. For the advanced codes the development of a 
more-detailed swelling and gas-release model is being 
actively pursued. Two approaches are presently being 
adopted. The first, a statistical approach,”? uses the 
Monte Carlo technique to simulate the behavior of gas 
bubbles. The second solves the differential equations of 
the evolution of bubble-size distributions by numerical 
approximations.®® The success of this development 
depends on the ability of these codes to calculate the 
effect of re-solution and various interactions between 
bubbles and microstructural defects. Solid fission 
products also will contribute to fuel swelling, but these 
are less important than gaseous fission products.*? 


Creep. The creep deformation of the fuel is con- 
sidered only in codes where detailed mechanical 
analysis is performed. The creep strain of fuel consists 
of a thermal component and a fission-fragment-induced 
component, similar to cladding. Thermal creep is 
strongly temperature dependent with an activation 
energy®? of ~100 kcal/mole. A review of the relatively 
few data on fission-induced creep** indicates that at 
temperatures below 1000°C (1830°F) the process is 
essentially athermal. Low-temperature fission-induced 
creep rates are many orders of magnitude greater than 
thermal-creep rates. At higher temperatures, fission- 
induced creep shows a similar temperature dependence 
to thermal creep, but at rates approximately an order 
of magnitude higher. 
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rapidly at first and then slows down. The detailed 
results to date are in static sodium, and oxygen 
removal in flowing sodium can limit the swelling to 
values much lower than those observed in these 
experiments. The observations on failed fuel pins in 
flowing sodium (such as the BR-S failures discussed 
earlier, in which the reaction was not observed with 
PuO, fuel) are evidence for the dominant effect of 
oxygen removal under flowing conditions. Further 
work is needed to confirm this in mixed-oxide fuel- 
pin-failure tests under flowing-sodium conditions. 


ANALYTICAL METHODS TO DESIGN 
AND PREDICT PERFORMANCE 
OF FAST REACTOR FUEL PINS 


In the present state of fuel development, consider- 
able reliance is placed on irradiation tests and their 
interpretation. More information is required on fuel 
pins, subassemblies, and materials operating at the 
conditions expected for LMFBR designs. Progress in 
screening various fuel designs and establishing the 
effects of power ratings, temperatures, and burnup is 
limited owing to: (1) lack of high fast-flux test 
environments, (2) the length of irradiation time re- 
quired to obtain significant results, and (3) limited 
usefulness of thermal-flux tests. 

In parallel with the design screening and the 
irradiation experiments, an analytical approach is 
required to establish guidelines to assist in the defini- 
tion of irradiation programs and for the development 
of optimum fuel designs. When the fuel pin and 
assembly are treated as an interdisciplinary system to 
describe fuel performance during normal and off- 
normal operation, the following benefits result: 

1. Guidance and direction for irradiation testing. 

2. Extrapolation and interpolation from known to 
unknown regions of performance. 

3. Optimization of fuel designs. 


The recognition of these advantages had led to 
concentration on the development of fuel-performance 
codes over the past 3 years, particularly in the U.S. 
program.°® 

In more detail, the fuel-element-performance 
models are being used in several ways. The models 
force the materials experts to make quantitative 
predictions on fuel and cladding performance. This 
indicates areas where experimental and theoretical 
work are most needed. Comparison of model predic- 
tions with each new batch of fuel-element-performance 
data may indicate inadequacies in the medels used but 
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also indicates the degree of agreement in the perfor- 
mance of different fuel elements. Correlating the large 
quantity of performance data on fuel from the EBR-II 
program is the most important near-term benefit. As 
the models develop and experience is gained, they are 
being used to determine fuel designs and to predict the 
behavior of elements in the LMFBR. The most 
challenging objective is the use of the codes to predict 
the effect of various transients and off-normal opera- 
tion on fuel behavior. 

The models developed to describe the operation of 
fuel elements under normal conditions can be arbitrar- 
ily divided into two categories. In the first, the models 
are more phenomenological and consist of fast-running, 
largely steady-state models that use empirical relations 
to describe fuel behavior; for example, the relation 
between gas release and fuel swelling vs. burnup and 
temperature. Only the more observable characteristics 
of the oxide core and the cladding are modeled. This 
type of model draws heavily on the physical aspects of 
the fuel obtained from macroscopic measurements and 
metallography. Since the analytical models are consis- 
tent with the empirical observations, the codes are 
valuable for interpolation and limited extrapolation. 
OLYMPUS (Westinghouse)*® and SWELL (Argonne 
National Laboratory)®? are examples of these models. 
The work at General Electric Company,’° Atomics 
International,’?! and Hanford Engineering Develop- 
ment Laboratory’? also belongs to this category. Many 
of these codes have detailed subroutines, and there is 
an overlap with the second category in which more 
realistic models of the various microscopic processes 
that are believed to control the behavior of the fuel 
element are formulated. The main objectives of these 
detailed codes are to relate materials properties and 
physical mechanisms to the widest possible range of 
fuel-element geometries, materials parameters, and 
operating conditions. Because of the firmer physical 
basis for the models, extrapolation over a wider range 
of conditions should eventually be valid. Also, they can 
be used to guide studies of fuel and cladding properties 
since these data are needed as input. CYGRO,7° 
LIFE,”* and FMODEL’® are examples of these codes. 
In conjunction with their development, considerable 
effort has been given to the subroutines of the code 
that analyze a particular mechanism of interest. Fuel- 
swelling and gas-release subroutines, for example, are 
areas being pursued actively by several investigators.’ ° 

Analysis of the off-normal behavior of mixed-oxide 
LMFBR fuel is another area of active study since 
knowledge of fuel performance under such conditions 
is necessary for setting realistic design and operating 
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limits and establishing margins for unacceptable dam- 
age or failure. 

Characterization of failure models, pressures pro- 
duced by failure, and movements of fuel and coolant 
are needed for accident analyses.77 Experimental 
guidance is provided by the program in TREAT”” and 
other facilities for simulating LMFBR power excur- 


sions, loss of coolant flow, and short-term effects of 


local abnormalities (such as coolant-channel block- 
ages). 

Relatively long-term effects, such as fission-product 
migration or release, fuel restructuring, and degrada- 
tion of cladding mechanical properties, enter the 
transient analyses only as initial conditions. However, 
related phenomena arise, including, respectively, tran- 
sient pressures produced by fission-gas release upon 
transient heating, internal movement of fuel resulting 
from transient thermal expansion and melting, and 
cladding response to transient loading. 

The central problem in modeling fuel-element 
performance is the determination of the restructuring 
and shape changes of the fuel and cladding. The 
treatment of the shape changes of the cladding is 
relatively straightforward. The behavior of the mixed- 
oxide fuel is more difficult to describe in detail. As a 
result, the main difference between the various codes is 
the degree to which they model the changes taking 
place in the fuel. The more-advanced codes have 
subroutines that treat each of the following processes. 


Cladding Swelling. In mixed-oxide fuel elements 
with low smeared density (~85% of the theoretical 
density) and moderate burnup (~50,000 MWd/metric 
ton), little internal pressure is exerted on the stainless- 
steel cladding, and its dimensional changes are deter- 
mined by fast-neutron-induced swelling. The empirical 
correlation most often used in the United States for 
316 stainless steel has been given previously. 


Cladding Creep. The creep of fuel-element cladding 
is the sum of two parts: thermally activated creep and 
radiation-induced creep. Thermally activated creep 
data are taken from out-of-pile tests. The radiation- 
induced component is superimposed on this and found 
to be dominant at lower stress and temperature,’® as 
discussed previously. 


Cladding Corrosion. Corrosion affects cladding in 
at least three different ways: reduction in thickness 
due to sodium cladding attack;®? loss of strength due 
to removal of interstitial elements, such as carbon;*” 
or fission product attack of the cladding.'® However, 
the empirical data on the kinetics of these reactions are 
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quite limited at the higher sodium coolant tempera- 
tures where the attack is more severe. The data on 
fuel—cladding reaction have been summarized earlier, 
and attention has been drawn to the importance of the 
O/M ratio of the fuel and to the urgent need for more 
data in this area. For conservative design, allowances 
are made in terms of decreased effective cladding 
thickness. As reliable kinetic equations are developed 
for these reactions, they can be incorporated into the 
codes. 


Fuel. The low thermal conductivity of the oxide 
causes steep temperature gradients in the fuel. Conse- 
quently there are great differences in many of the 
properties of the fuel as a function of radial position. 
The changes that occur in the fuel are primarily 
changes in dimension due to swelling, creep, hot 
pressing, thermal expansion, and fuel cracking. The 
rates of all these kinetic processes change greatly with 
temperature. 


Swelling and Gas Release. 1n the phenomenological 
codes the fuel swelling and gas release due to gaseous 
fission products are calculated using empirical rela- 
tions. For the advanced codes the development of a 
more-detailed swelling and gas-release model is being 
actively pursued. Two approaches are presently being 
adopted. The first, a statistical approach,’® uses the 
Monte Carlo technique to simulate the behavior of gas 
bubbles. The second solves the differential equations of 
the evolution of bubble-size distributions by numerical 
approximations.2® The success of this development 
depends on the ability of these codes to calculate the 
effect of re-solution and various interactions between 
bubbles and microstructural defects. Solid fission 
products also will contribute to fuel swelling, but these 
are less important than gaseous fission products.*? 


Creep. The creep deformation of the fuel is con- 
sidered only in codes where detailed mechanical 
analysis is performed. The creep strain of fuel consists 
of a thermal component and a fission-fragment-induced 
component, similar to cladding. Thermal creep is 
strongly temperature dependent with an activation 
energy®? of ~100 kcal/mole. A review of the relatively 
few data on fission-induced creep®? indicates that at 
temperatures below 1000°C (1830°F) the process is 
essentially athermal. Low-temperature fission-induced 
creep rates are many orders of magnitude greater than 
thermal-creep rates. At higher temperatures, fission- 
induced creep shows a similar temperature dependence 
to thermal creep, but at rates approximately an order 
of magnitude higher. 
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Restructuring and Composition Changes. The 
chemical and the physical states of the mixed-oxide 
fuel will change as a function of burnup. The phe- 
nomena of particular interest include sintering (hot 
pressing), grain growth, uranium and plutonium re- 


distribution, and changes in the O/M. At present the 


treatment of restructuring is adequate in the advanced 
codes,”75 but further development is required to treat 
the kinetics of plutonium and oxygen redistribution. 


Cracking. A realistic fuel-element model must in- 
clude the effects of fuel cracking and requires the 
conditions (temperature and stress) of fuel cracking 
and its geometry. The mechanisms and the rates of 
crack healing must also be included. Attempts are 
being made to predict fuel cracking based on the 
mechanical data of the fuel.*4 


Thermal Analysis. In the codes for normal opera- 
tion, the methods of calculating temperature distribu- 
tion are similar. There are four major areas where 
uncertainties are involved, and the treatment differs in 
various models. These areas include: 

1. The effects of fuel cracking, restructuring, and 
composition changes on thermal performance. 

2.The values for the conductance of the fuel- 
cladding gap. 

3.The treatment of fuel melting and its conse- 
quences. 

4. The effect of as-fabricated porosity on thermal 
conductivity. 


In treating severe transients, calculations of fuel 
response are generally based on calculations of tem- 
perature changes. Simple R-Z geometry, as in 
ARGUS®* and SAS,°° is most common, although a 
more general, three-dimensional capability, such as that 
in TOSS,87 the THT series,2® or TRUMP,®? is 
required for many cases. Phenomena modeled include 
changes in the fuel—cladding gap conductance as the 
fuel expands,®® a hysteresis effect in the gap conduc- 
tance due to cracking,®® spatial variation in power, and 
thermal expansion. Because severe temperature excur- 
sions are involved, fuel-pin temperature calculations are 
coupled to coolant boiling, and transient heat transfer 
is coupled to structure. 


Mechanical Analysis. The mechanical analysis for 
normal operation calculates the deformation of the 
fuel element and the stress level in the cladding and the 
fuel. In all the present models, the axial-symmetric 
boundary conditions are used. Attempts are being 
made to treat the three-dimensional case. In the 
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more-detailed codes,”7*"75 the capabilities to follow 
the power history exist. In LIFE the fuel pin is divided 
into a limited number of concentric rings. The differ- 
ential equations for each ring are solved algebraically. 
This allows a much faster iteration procedure than is 
possible with finite-element codes. In CYGRO a 
finite-element method is used, but no iteration is 
performed. If the fuel is cracked, the mechanical 
analysis becomes more complicated. The cracked fuel 
will cause local stress and strain concentrations in the 
cladding and alter the constraints for fuel swelling. This 
is an area of active research at the present. 

The transient codes are similar, the differences 
being mainly due to the high rate of loading. The 
driving force for mechanical deformation comes from 
the transient pressure produced by fission-gas release, 
differential thermal expansion, and volume change 
resulting from fuel melting. The advanced version of 
SAS®® attempts to treat these cases. The PECT code?! 
is a finite-element analysis of the cladding response to a 
given localized load and is useful for certain situations. 
The materials behavior, including that of both the 
cladding and the fuel under high rates of deformation, 
is important to the mechanical analysis in the transient 
codes. 


Failure Criteria. A fuel element is considered to fail 
when the cladding leaks or ruptures. For normal 
reactor operation the goal is to operate at a burnup 
limit that leads to very few failed elements—thus the 
ability to predict a failure rate of roughly 10° is 
needed. This would require tests of many elements 
under each set of conditions. One possible alternative is 
to find some reasonable parameter, such as cladding 
thermal-creep strain, that is proportional to failure 
probability in this range.?? This is the approach being 
used as the design basis for the FFTF fuel element. 
Radiation-induced creep is similar to the superplastic 
behavior exhibited by many materials. Thus there 
would be no ductility limit for the irradiation-induced 
component. This has led Boltax*? to postulate that 
fuel-element failure will occur after the thermal-creep 
component reaches a certain limit. Li and Yaggee?? 
have considered the various modes of cladding failure 
and their relation to the method of loading and have 
postulated other failure criteria. 


For transient behavior, cladding failure is predicted 
when transient strain reaches a predetermined level. 
The stress needed to give this level is strongly tempera- 
ture dependent, and failure is induced more as a result 
of a thermal excursion than of a pressure excursion. 
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Applications to Normal Conditions 


Most of the applications to date have compared 
code predictions with experimental observations, and 
discrepancies between prediction and observation have 
led to new experiments and revised models. Such codes 
are now being applied for sensitivity and design studies. 

In mixed-oxide fuel elements, the simplest ele- 
ments to treat unambiguously are those with relatively 
low smeared density and moderate burnup. Jackson 
et al.7* have examined the diameter profiles of two 
mixed-oxide fuel elements (85 to 88% smeared den- 
sity) irradiated at about 10 kW/ft peak power in the 
EBR-II to a peak burnup of 46,000 MWd/metric ton. 
They found that the profile can be described satisfac- 
torily by the predicted swelling of the stainless steel, 
indicating that the force exerted on the cladding by the 
fuel is not adequate to cause creep. These authors also 
find that the observed diameter changes peak more 
sharply at the core midplane than the calculated curve. 
Honekamp and Golden?* have observed the same 
effect for unstressed capsules irradiated in the EBR-II. 
Their calculations show that the effect can be ac- 
counted for by the softening of the neutron spectrum 
in moving away from the core midplane. 

Even though mixed-oxide fuel may not initially 
deform the cladding, it undergoes restructuring by 
migration of pores up the temperature gradient toward 
the center of the fuel. A treatment of this dynamic 
restructuring and its effect on thermal conductivity is 
essential in considering questions such as the power to 
cause incipient melting, especially in fresh fuel. Fig- 
ure 14 shows a comparison between predicted and 
measured radial porosity distributions for a model 
developed by Homan and Cox.7° The data are for a 
Sphere-Pac fuel pin irradiated at 13.6 kW/ft to about 
0.7% fission per initial heavy-metal atom in a thermal 
flux. There is good agreement between the experi- 
mental points and the calculated curve except at the 
center of the pin, where more restructuring appears to 
have occurred than was predicted by the model. 
Another check on the model is the comparison of the 
predicted and measured void diameters in the columnar 
and equiaxed grain regions. This comparison has been 
made for several pins from the General Electric F2 
series of irradiation experiments,?* and agreement is 
good for the diameters of the central-void, columnar, 
and equiaxed regions.?° 

Most of the early codes in use assume that the 
element operates under steady-state conditions. Actu- 
ally the tests are carried out in the EBR-II, where the 
reactor is shut down for fuel handling about once a 
month. 
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Fig. 14 Comparison of the calculated and measured radial 
variation in porosity for (Uo.ss, Puo.:5)O. irradiated at 
13.6 kW/ft to 0.7% fission per initial heavy-metal atom. 


The LIFE or CYGRO codes allow the calculation 
of the effect of operating history on the cladding 
deformation. Table 13 presents the results of a sample 
calculation for an actual fuel element in the EBR-II for 
three modes to give the same burnup: the actual 
history, steady state at the peak power, and steady 
state for the time the reactor was operating. The code 
was calibrated to make the predictions for the actual 
history agree with the experimental values. The clad- 
ding creep strain predicted for steady-state operation at 
either peak power or average operating power is below 
that for the actual operating history. This results from 
the stress transients accompanying startup and shut- 
down during power cycles. 

The most-frequent and most-fruitful experimental 
test of the codes now being performed involves a 
comparison of the predicted and experimental profiles 
for well-characterized irradiated fuel elements. When a 
code has been “fine tuned” to fit the results of a series 
of fuel-element experiments, it can be used to predict 
the performance of similar elements over a wider range 
of conditions. 

The OLYMPUS code has been used to correlate the 
behavior of the General Electric F2 series and the 
NUMEC high-burnup pins.*® Figure 15 illustrates the 
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Table 13 LIFE-1 MiArlane Results* for Fuel Element XG-5—F2H 








Actual Steady Steady operat- 
history} peak power ing power 
(911 days) (284.5 days) (369 days) 
Experimental (542 days 45 MW or 35.96 MW or 
values at 0 kW/ft) 14.5 kW/ft 11.6 kW/ft 
Total fluence, 
107? neutrons/cm? ~0.47 0.47 0.47 0.48 
Cladding swelling 
strain, % 0.72 0.72 0.78 
Cladding creep 
strain, % 0.46 0.21 0.0 
AD/D, % 
(final) PAZ 112 0.92 0.78 
Outside radius of 
columnar zone, in. Varied 0.072 0.069 0.039 
Fission gas 
released, % 15 75 15 74 





*Fuel, mixed-oxide pellets (20 wt.% plutonium); cladding, solution-treated 316 stainless 
steel; burnup, 6.75 at.%. 
+Actual history predictions calibrated to agree with experimental values shown. 





predictions that can be made. The data points are the 10 ST ea: eee or 
final observed cladding strain for elements C-15, B-2, Pin ly rp Thickness, in. % AD/D7 ,C-1 
and C-1, The curves are calculated. Element B-2A is a C151 83 0.015 DB 
hypothetical element. The rapid rise in strain with B-2A| 83 0.030 a 
burnup stems from the elimination of available poros- B-2 | 88 0.030 a) y 
ity within the fuel. The continuing fuel swelling then 8)c-1 | 88 0.015 ° 1 / ” 
raises the stress in the cladding, and thermal creep f+ 70-15 
becomes more important than radiation-induced creep 
(lower curves in Fig. 15). Two effects are of special 
interest. The cladding strain rate in this range is 
influenced more by a 5% decrease in smeared density 
than by doubling the cladding thickness. In addition, 
the thermal-creep component of the strain rises much 
more rapidly for the higher-smeared-density elements. 
Since the thermal-creep-ductility limit is only 1 to 2%, 
failure of elements C-1 and B-2 would be predicted at 
about 14% burnup and satisfactory performance of 
C-15 to a significantly higher level. a 
As another example ; of ~ analysis * me a 7 78? 

elements, the results obtained at Argonne National A rie 
Laboratory on an encapsulated NUMEC element C-11 i iii: 7 
can be considered.’7* The profiles of the restructuring aac | > 

zones in the fuel could be fitted well with the LIFE a yy /B-2A _ 


code, but failure to fit the cladding strain profile 4 1 OC-15 
proved to be more informative. The element was clad ‘ j i ae” \ 
in solution-treated 316L stainless steel. The initial 9 10 11 12 13 14 15 
fuel—cladding diametral gap was 1.2 mils in the lower BURNUP, % 

half of the element and 2.5 mils in the upper half. The 

smeared fuel density in the two halves was 88.0 and Fig. 15 OLYMPUS analysis of high-burnup fuel pins. 
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88.6%, respectively. Figure 16 shows the observed and 
predicted profiles for the unstressed capsule and the 
fuel element, where the contribution due to fuel 
swelling has been increased in the code so that the 
deformation may be predicted using a steady-state 
analysis. The flux-induced swelling of the stainless steel 
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Fig. 16 LIFE analysis of high-burnup pin, NUMEC element 
C-11. 


in the capsule and in the element cladding is about the 
same, and so there is clearly about 1% cladding strain 
due to the fuel pushing on the cladding. The predicted 
diameter shows a discontinuity 7 in. from the bottom, 
where the initial gap changed. The observed profile, 
however, shows no such discontinuity. The LIFE code 
at that time did not allow for fuel cracking, and thus it 
predicted that the cool outer portion of the fuel pellet 
developed a tensile stress and restrained the fuel 
expansion. In fact, very early in life the fuel cracks, 
relieving these tensile stresses, and thus increases the 
fuel swelling rate, so that the fuel moves out to contact 
the cladding. This has led to the inclusion of fuel 
cracking in the model. 

Thus to date the application of analytical models 
for steady-state conditions has led to an appreciation 
of the importance of stainless-steel swelling and irradia- 
tion creep in fuel-pin performance and is providing 
detailed mechanistic evidence for the reliability and 
endurance to high burnup of the mixed oxide/stainless 
steel system. 


Application to Transient Conditions 


At present the main effort in developing fuel- 
element-performance codes has aimed at interpreting 
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results from elements run under steady-state con- 
ditions. This will continue to be important, but even 
more valuable will be predictions on the effect of 
sudden transients or off-normal operation at various 
points in the life of the fuel element. It is impossible to 
test fuel-element performance under all these con- 
ditions, yet the accepted fuel-element design must be 
one that can withstand reasonable excursions at any 
time in its expected life. 

In the study of transient behavior, the objective 
changes to reactor safety. Power transients are uninten- 
tional, and, if they occur, the fuel-element and subas- 
sembly designs must be able to contain or localize the 
fuel. In addition, the consequences of failures must be 
characterized well enough so that possible conse- 
quences can be calculated realistically. 

The first major questions in transient experiments 
are thus: “When will the fuel cladding rupture, and 
what is the failure threshold?” All indications to date 
are that the cladding will not rupture in the absence of 
fuel melting. For example, Culley et al.?7 subjected a 
mixed-oxide fuel element (preirradiated to 
8300 MWd/metric ton in the EBR-II) to a transient 
that calculations indicate brought 50% of the fuel to 
within 500°F of the assumed 5000°F fuel-melting 
temperature. Ceramography confirmed that no melting 
occurred. The cladding diameter increased from 0.2 to 
0.5%, but the cladding did not crack. Differential 
thermal expansion between fuel and cladding is be- 
lieved to be responsible for the cladding deformation. 

When the transient leads to fuel melting,?*® clad- 
ding strain is greater, and cladding rupture becomes 
more probable, although the use of low smear densities 
in the LMFBR fuel confers an advantage here. An 
example of tests in which fuel melting occurs is given 
by the first tests conducted on oxide-fuel meltdown in 
the TREAT sodium loop, tests El and E2.?? These 
tests were the first conducted on failure-threshold 
phenomena and postfailure events for LMFBR-type 
oxide fuel under flowing-sodium conditions. They 
involved a power excursion, with a flowing coolant to 
produce typical axial temperature rises and a flow- 
channel geometry for coolant expulsion and reentry. 
The El pin contained a 27.9-cm-long stack of UO, 
pellets of 90% smeared density inside a 304 stainless- 
steel jacket. The E2 pin was similar to El, except that 
it contained a depleted oxide pellet and a 5-cm-long 
steel rod above the fuel stack. Both tests—run with 
fresh fuel and power pulses of about %4-sec duration— 
produced in excess of 80% fuel melting by volume. The 
unfailed El pin had a cladding strain of <0.3%, 
compared with a strain calculated using the DEFORM 
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code in excess of 3%. About 1% of the predicted 
cladding deformation came from differential expansion 
of cladding and fuel before melting and about 2% from 
the volume change resulting from fuel melting. DE- 
FORM is a two-dimensional code that does not allow 
for axial motion of the fuel. Postirradiation inspection 
clearly indicated axial motion of the molten fuel. The 
E2 pin failed at an energy input about 10% greater 
than that giveg the El pin. Measured cladding strain 
from the E2 inspection was 1.4%, in good agreement 
with 1.4% strain calculated up to the onset of fuel 
melting. 

Given experimental data and a basic analytical 
capability, more-detailed modeling can proceed. Fig- 
ure 17 shows the predicted response of an LMFBR 
oxide fuel element to a power transient after extended 
prior steady-state operation, as calculated using the 
improved DEFORM! °° code. Power rises from zero to 
0.48 kW/cm in 100 sec, then to 0.72 kW/cm in another 
3 sec. As the power increases, the gap closes, A, loading 
the cladding elastically. The equivalent stress increases 
rapidly, and the outside of the cladding begins to 
deform plastically, B. By C, all the cladding is plastic, 
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and the fuel—cladding contact pressure begins to 
decrease. At D the innermost fuel starts to melt. This 
relieves the stress on the rest of the fuel so that it 
contracts, unloading the cladding elastically. Gas re- 
leased from melting fuel and the expansion of the 
molten fuel cause the pressure in the cavity to rise 
rapidly, until at E the cavity is almost filled with 
molten fuel. 

In this way, development of analytical capabilities 
to calculate the transient response of mixed-oxide 
LMFBR fuel is proceeding along with the progress in 
experimental work. The first stage can be described as 
a survey phase, in which the capability to calculate the 
conditions of transient exposure up to failure was 
confirmed, and the first experimentally based models 
of behavior were generated. Current emphasis is on 
better simulation of LMFBR conditions, paralleling the 
position for analytical efforts for steady state. The 
next step will be a stronger interaction between 
analysis and experiment to build the basic under- 
standing of the detailed behavior to make such 
analytical and experimental predictions reliable for a 
wide range of conditions. 
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Fig. 17 Predicted response to power transient using the DEFORM code. 
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FUEL-ASSEMBLY PERFORMANCE 


As indicated, the phenomena of cladding swelling 
and irradiation creep improve the fuel-pin performance 
since they extend the critical point (at which fuel 
cladding interaction due to fuel swelling occurs) to 
higher burnup. At the same time, the existence of these 
phenomena has created problems outside the fuel pin, 
and the behavior of the pin-spacer system and the 
interaction of the system with the duct have become 
the important areas in fast reactor fuel elements. The 
up-to-date position on these aspects may be reviewed 
as follows. 

Two systems for the support and spacing of fuel 
pins within a fuel assembly have been adopted: (1) the 
wire-wrap system employs one or more wires wound 
around the pins on a constant-pitch helix and fastened 
to the pin at its upper and lower end caps. In this 
system, each pin is dependent on its neighbors for 
transverse support, with only the outermost pins being 
supported directly by the flow duct; (2) the grid 
pin-support system provides direct support to each 
individual pin. Spacer grids are located at intervals 
along the length of the fuel assembly and are attached 
to the flow duct by means of hanger straps. Both 
spacer systems can impose loads on the fuel pins and 
thus induce cladding strains. 


Wire-Wrap Assemblies 


There is a great deal of irradiation experience on 
wire-wrapped LMFBR fuel assemblies, particularly 
from Rapsodie, as described earlier in the section titled 
“Irradiation Experience in Test Reactors.” In addition, 
the helical wire-wrap spacer has been chosen for the 
first cores for FFTF. After the fuel assemblies are 
exposed to high fluences, the problems arising with the 
wire-wrap spacer concept are: (1) the degree of 
compliancy to accommodate differential swelling be- 
tween the fuel-pin assembly and the duct tube and (2) 
the extent of differential swelling between the cladding 
and the wire since this can produce either wire 
breakage or slackness. 

At any given plane in any fuel assembly, both the 
duct tube and the cladding experience approximately 
the same average flux. If all stainless steel in the 
fuel-assembly cross section has common metallurgical 
characteristics, only the difference in operating tem- 
peratures could cause differences in swelling rates. With 
wire wraps and oxide fuel pins, the temperature 
difference between average duct-tube material and 
average cladding material at the midplane of the central 
fuel assembly will be in the range 100 to 150°F, 
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depending on the design details. At peak fluence in the 
early cores of near-term LMFBRs of 1.5 x 1073 
neutrons/cm*, a 150°F temperature mismatch would 
produce a maximum swelling difference of about 4% 
AV/V in the temperature range of interest. Assuming 
isotropy, this is equivalent to a difference in AL/L of 
about 1.3%. Cladding strain from fuel growth may 
bring the total planar differential expansion between 
cladding and duct tube into the 1.5 to 2.0% range. The 
wire wraps operate at a temperature closer to that of 
the duct tube, however, and suffer less differential 
swelling. Therefore, on the assumption that cladding, 
wire wraps, and duct tube obey the same swelling law, 
it appears that a 1.5% planar—differential—linear ex- 
pansion between a wire-wrapped assembly and its 
wrapper tube is a reasonable design value for near-term 
cores. There is no basic reason why a 1.5% differential 
AL/L cannot be accommodated by clearance in the 
as-fabricated fuel assembly, although it is higher than 
the nominal clearance generally provided to facilitate 
loading of the wire-wrapped pins. Emphasis is therefore 
being placed on the determination of a realistic upper 
limit of allowable stack-up clearances within the fuel 
assembly to establish how much differential swelling 
can be accommodated in this way. In addition, 
pin-assembly compression tests must be carried out. 
Balent and Hirst'®’ have reported the results of 
out-of-reactor tests that measured the inherent com- 
pressibility (springiness) of a wire-wrapped assembly 
and have shown that the assembly exerts some pressure 
on the duct when there is a loose spacing of the 
wire-wrapped pins. This pressure arises from the 
small-amplitude bowed axial shapes of the pins, which, 
when compressed, behave as a multiplicity of low- 
spring-rate members. This effect thus permits provision 
within the assembly of fuel pins to accommodate 
variable amounts of differential swelling between the 
pins and the ducts; it also exerts a positive force to 
reduce the tendencies for vibration and consequent 
fretting and wear. 

To evaluate the question of wire slackness or 
breakage, one must estimate the differential growth of 
the wire and the cladding. The wire is exposed to the 
same axial flux profile as the cladding and can 
presumably be made of the same material and degree 
of cold work. The average temperature of the wire in 
the core region is lower than the average cladding 
temperature by 50°F. Huebotter et al.'°? have esti- 
mated the effects caused by this differential. Thus 
Fig. 18 shows a typical axial trace of helical-wire and 
midwall-cladding temperatures and a typical axial flux 
profile. With present swelling estimates at a fluence of 
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Fig. 18 Typical axial trace of helical-wire and midwall-clad- 
ding temperatures, 


1.5 x 107? neutrons/cm?, the integrated volumetric 
growth of cladding is 7.1%, compared with 6.0% for 
the wire. Again assuming isotropy, the wire will 
elongate 2.0%, and the length to be wrapped will 
increase by 2.4%, giving the differential as 0.4%. 
Additional deficits in wire length are created by 
differential thermal expansion and fuel swelling, which 
raise the total to about 0.5%. Because the wire wraps 
the fuel element along its entire length, of which the 
core accounts for only about 40%, the actual deficit in 
wire length is only about 0.2%. Thus an_ initial 
slackness amounting to 0.2% of the total length of the 
wire should prevent end-of-life strain in the wire. Initial 
slackness is considered preferable to initial tightness 
since the wire may break or cause the fuel pin to 
distort when the wire tends to tighten further by 
differential elongation. On this basis, there should be 
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no assembly that could experience a slackening of the . 
wire wrap or a loosening of the wire-wrapped pins 
within the duct as burnup proceeds in the case 
considered. 

However, these assumptions!®? involve uncer- 
tainties that will only be resolved by irradiation of 
assemblies to high fluences. The helical wire-wrap 
spacer has been chosen for the FFTF fuel assembly; in 
addition, development programs are proceeding on the 
compliant grid-spacer concept. 


Grid Assemblies 


As indicated previously, the honeycomb-grid spacer 
system has been selected for the PFR fuel assembly, 
and a considerable amount of out-of-reactor develop- 
ment, together with irradiation of grid-type assemblies 
in the DFR, has been carried out in the United 
Kingdom. 

Out-of-reactor testing has developed data on mix- 
ing with different spacings of the honeycomb grids and 
has provided indications of the relative mixing level 
with and without grids. Data have also been obtained 
on the correlation factor between mixing in air and in 
sodium, on the pressure-drop characteristics of grid 
assemblies, and on the flow distribution in the sub- 
channels of assemblies. 

Mark 1, the first oxide-grid test assembly, was 
placed in the DFR in late 1963 and reached a burnup 
of 5.5%. There has been a series of subsequent 
assembly tests of the 77-pin type, using honeycomb 
grids. These tests have incorporated a large variety of 
pin—grid variables and reached burnups of up to 7.3% 
at a midwall cladding temperature of up to 1290°F. 
Experience to date has therefore been satisfactory in 
reactor with the grid-spacer concept. 

Significant out-of-pile grid tests have also been 
performed elsewhere.'°* These tests have yielded 
extensive data on flow distribution, fretting wear, 
vibration, and dimensional stability at high tempera- 
tures at test periods up to 4000hr. Vibration and 
fretting-wear data are available over a wide range in 
fuel pin—grid cells of different clearances. 

In the United States the program! °* 
series of 37-pin grid-type assemblies for irradiation in 
the EBR-II, leading to the irradiation of two full-size 
217-pin grid-type assemblies for irradiation in the 
FFTF. In 1971, the first 37-pin grid subassembly, 
Westinghouse Subassembly No. 3 (WSA-3), was under 
irradiation in the EBR-II; and, in January 1972, 
General Electric!®* followed with a similar design 
concept in its F9E Subassembly. Two specific require- 
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ments had to be met prior to the WSA-3 irradiation 
test. 

1. Satisfactory development of the grid and subas- 
sembly hardware. 

2. Demonstration of the pressure-drop and flow 
characteristics for this particular subassembly design. 


In addition, development of the grid design for 
later subassemblies leading to the full 217-pin grid 
subassemblies for insertion in the first FFTF core 
required the design and implementation of a series of 
fretting tests to determine the optimum pin—grid gap. 

The first objective of development of grid-and- 
subassembly fabrication techniques has been ac- 
complished, and a wide range of grids to various 
tolerances can now be made cheaply and reliably. 
Inspection techniques have been developed in parallel. 

For the flow-distribution work, an assembly was 
constructed which was identical to the first to be 
irradiated except for substitution of tungsten carbide 
pellets for the fuel. This assembly was instrumented for 
pressure drop and flow-distribution measurement and 
testing in water. The main purpose of the tests was to 
dimension an integral orifice plate so that the available 
head of sodium in the EBR-II reactor will produce the 
design coolant flow rate. A comparison between 
calculated and measured flow coefficients of the 
assembly and of one grid was also made. The flow- 
distribution tests gave checks on the accuracy of the 
subchannel flow-distribution calculations made with 
the COBRA code. Calculations and flow tests show 
that the orifice plate with seven 0.198-in.-diameter 
sharp-edged holes gives the desired assembly pressure 
drop of 41 psi at a sodium flow rate of 33.7 gal/min. 
Tests at different water temperatures confirmed the 
correctness of the extrapolation methods used for fluid 
flow of different viscosities and densities. The flow 
distribution measured by the pitot-tube method at the 
assembly outlet confirmed the subchannel flow- 
distribution calculations made by the COBRA code. In 
general, agreement between calculation and test results 
was good and provided the basis for the acceptance of 
the irradiation. 

The first fretting tests on grid-type subassemblies 
were completed in 1970. The assembly involved a fuel 
assembly utilizing pins and grids in which relative 
motion between the fuel pins and grids was achieved 
by differential thermal expansion. A cycle consisted in 
directing a 1200°F sodium stream up the fuel-pin 
channels while a 500°F stream was directed downward 
around the periphery of the fuel-assembly wrapper can 
and then reversing the flow directions. The results were 
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satisfactory.'°* A similar program using the wave- 
beam type of grid is also in progress.'°* 


Core Structural Components 


In LMFBR cores, swelling will cause the ducts of 
fuel and blanket assemblies and the guide tubes for 
control assemblies to elongate and dilate during service. 
Further, when the rate of swelling of duct and 
guide-tube material is not circumferentially uniform, 
distortions that can take the form of either simple 
bowing (isoplanar) or complex bowing (multiplanar) 
are possible. Circumferential variations in swelling are 
inevitable because of the flux and temperature gradi- 
ents across all but the central assembly. 

In evaluating the effects of swelling, Huebotter! °? 
took the Atomics International 1000-MW(e) LMFBR 
design'°® as a target core and concluded that dilation 
not accommodated by adequate clearances between 
assemblies would have many unacceptable conse- 
quences, including: 

1. The in-core portions of control-rod-guide and 
safety-rod-guide tubes would be displaced from their 
original positions. Long before nominal core burnup is 
reached, these tubes would suffer sufficient misalign- 
ment with their above-core extensions so that conven- 
tionally designed rods could not be inserted into the 
core. 

2. Handling heads of all but the central core 
assembly would move radially outward as the core 
accumulates fast fluence. At nominal discharge burnup, 
the handling heads of fuel assemblies at the core edge 
would be displaced by about 3 in, Fuel handling under 
these circumstances would be very difficult. 


3. The neutronics of the core would be affected 
greatly by what amounts to a 2.5% weighted (flux- 
squared) average growth in the cross-sectional area of 
the core. On the basis of rough estimates of the effects 
of core expansion on reactivity, the radial expansion 
can amount to a 3-dollar loss of reactivity that must be 
compensated by excess fuel and additional control 
rods. Core elongation also produces a reactivity loss, 
estimated at over 2 dollars. 


The most practical way to accommodate dilation 
of both duct and guide tubes is to provide gaps 
between the tubes. The width of the gap is dictated by 
swelling in the worst lattice position—probably the 
fuel assembly at the core center line. In the Atomics 
International LMFBR design, by provision of gaps, this 
increases the diameter of the core 4% and increases the 
sodium-to-fuel ratio from 0.86 to 1.2. This results in an 
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increase in critical mass and a decrease in breeding tive assembly in each of the rows of the Atomics 

ratio. International 1000-MW(e) reactor design.'°? Free 
It is also advisable to move the spacer pads to bowing, i.e., neglecting mechanical interferences with 

regions of the duct that will not suffer dilation. neighbors, is shown in Fig. 19. 

Without a plane of tightness through the core, the If the outward free bowing were suppressed by 

ability to predict changes in radii of fuel- and control- top-end restraint, the bowing profiles would be like 

assembly center lines is diminished. those in Fig. 20. This shows a strong inward bow of the 


fueled length of about 0.46in. in the outer core 
assembly and more in the inner blanket assembly 
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locating the spacer pad, as was done by Atomics 
International, 10% in. above the core midplane. It is 
evident that this is effective in counteracting positive 
reactivity effects. As discussed previously, however, the 
accommodation of dilation is better assured without 
in-core spacer pads. 

Uncertainty in the magnitude of radiation-en- 
hanced creep, as well as the magnitude of swelling, 
greatly influences the confidence in predicting the 
locations of the center lines of the assemblies and the 
control-rod-guide tubes and the locations of the han- 
dling heads when the top-end restraint is removed. 

Except for the influence of creep, a fuel assembly 
removed from the restraining effect of its neighbors 
will tend to assume the free-bowed shape shown in 
Fig. 19. This could create difficulty if subsequent 
operations are designed to accommodate only straight 
components. Therefore the ease of handling of fuel 
assemblies and guide tubes will be sensitive to the 
amount of straightening accomplished by radiation- 
enhanced creep. 

The FFTF is the first reactor to be designed with 
knowledge of the stainless-steel swelling and irradiation 
creep. On this project, engineers from HEDL, WARD, 
and AEC are pioneering the interpretation of the 
design implications, evaluating engineering solutions to 
the problems, and developing specialized analytical 
solutions to the problem. The approaches have been 
described by Pennell,’°’ and the evaluations and 
sensitivity analyses that contribute have been published 
by Jackson and coworkers.!°*%:!°° In this work the 
role of irradiation creep in reducing bending stresses 
effected by core-restraint arrangements, and also in its 
effects on fuel handling when the core is unrestrained, 
has been clearly demonstrated. The general arrange- 
ment of the radial restraint system and the fuel-as- 
sembly arrangement have also been described.' ? 

The problems introduced outside the fuel pin by 
irradiation-induced phenomena, such as swelling and 
irradiation creep, can therefore be minimized in a 
number of ways: (1) by the choice of the spacer 
support system, (2) by providing gaps between the 
fuel-assembly ducts, and (3) by proper design of the 
core-restraint system, taking into account the implica- 
tions of the role of swelling and creep. 

The economic consequences of swelling have to be 
considered in relation to the design solutions to 
accommodate it. Common to most solutions is the 
need to provide interassembly gaps while retaining an 
end-of-life margin for fuel handling. The resulting 
larger core and fuel-fraction loss lead to a penalty on 
energy costs in typical cases'®? of about 0.1 mill/kWh. 
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There is therefore a large incentive in the longer term 
to solve the swelling problem metallurgically, as indi- 
cated in the section titled “Effects of Fast-Flux 
Irradiation on Cladding and Structural Materials” of 
this review. 


ADVANCED FUELS 


Although all fast breeder fuel-development pro- 
grams are concentrated on the mixed-oxide system, the 
mixed carbides and nitrides offer potential for the 
longer term, and Washburn and Scott'!° have recently 
reviewed their performance capabilities. 

The carbides and nitrides of uranium and pluto- 
nium have a thermal conductivity about five times that 
of the oxide. This high thermal conductivity permits 
operation of fuel pins at higher linear and specific 
powers and enables design flexibility for optimization 
of both the cost of energy and the fissile doubling 
time. The marked advantages of advanced fuels were 
brought out in the recent 1000-MW(e) LMFBR follow- 
on studies.'°® The three oxide-fuel designs each had 
linear heat rates of 16 to 18kW/ft peak and 9 to 
10 kW/ft average. The three sodium-bonded carbide 
designs from these studies had linear heat rates of 27 to 
43.5 kW/ft peak and 16 to 28.5 kW/ft average. Of the 
three oxide-fueled designs, the lowest doubling time 
was 7 years, with a fuel-cycle cost of 0.38 mill/kWh. 
The sodium-bended, carbide-fueled replacement core 
for this design reduced the doubling time to 4.4 years 
and the fuel-cycle cost to 0.11 mill/kWh. Assuming 
similar burnup capabilities, the carbides have a clear 
potential economic advantage relative to oxides. 
Another study!!! has indicated that the economic 
potential of the nitrides would be similar to the 
carbides. 

There are a number of barriers that must be 
overcome prior to commercial use of carbides and 
nitrides. Since these fuels react with air, they must be 
fabricated in an inert atmosphere of high quality. To 
achieve this end, one must take greater care during the 
design and installation of fabrication facilities for 
carbides or nitrides, but subsequent operation is not 
greatly different from that for mixed oxides. The 
stoichiometry of the mixed carbide must be carefully 
controlled; hypostoichiometric carbide swells exces- 
sively, whereas hyperstoichiometric carbide supplies 
carbon to stainless-steel cladding, and sodium bonding 
accelerates the rate of carbon transport. The effects of 
stoichiometry on the performance of the nitrides have 
not been clearly established, but in any event the 
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stoichiometry of the nitride is more easily controlled 
than that of the carbides. Vaporization rates of the 
nitride have been found to be lower than those of 
carbides or oxides, which abates earlier postulations 
relative to severe swelling from free metal or attack of 
the cladding by free nitrogen in the irradiated fuel. 

Irradiation tests on mixed carbides and nitrides 
have shown that swelling rates are acceptably 
low,!!?"1!4 especially when sodium bonding is used 
to keep fuel temperatures below ~2370°F or operating 
at linear heat rates to ~40kW/ft. Problems from 
fission-product redistribution or stoichiometry changes 
due to higher burnups have not been observed. 
Cladding failures have been observed with helium- 
bonded pins when the initial fuel—cladding gaps were 
small. In sodium-bonded capsules with large internal 
gaps, other problems have been observed from fuel 
breakups producing chips, which in turn produce 
fuel—cladding mechanical interactions. Overall results 
are promising, but more steady-state and transient 
irradiation tests are needed for both fuels and, in 
particular, for nitrides, before large-scale development 
programs can be undertaken. 


SUMMARY 


1. To date, mixed-oxide fuel-pin irradiations in test 
reactors have demonstrated the high-endurance capa- 
bilities of the mixed oxide/stainless steel system. 
Fission-gas release and fuel swelling can be accom- 
modated, and the benefits arising from the radiation- 
induced effects of swelling and creep in the stainless- 
steel cladding give confidence that the high-burnup 
targets for LMFBR fuel can be achieved. 

2. Fuel—cladding reaction has emerged as a poten- 
tial problem that may limit the peak cladding tempera- 
tures to which mixed-oxide fuel pins are designed and 
operated. Further detailed results on these limits are 
needed. Fuel—sodium reactions do not appear to 
present a problem under dynamic coolant conditions. 

3. Definition of the thermal-performance models 
requires more data, and statistical testing is needed to 
differentiate between the variations in design judgment 
that exist today to optimize for the reliable economic 
performance required. 

4. Stainless-steel swelling and irradiation creep 
have, in effect, focused attention on the assembly, and 
the behavior of the pin spacer system and the duct are 
the areas in which more-precise materials-irradiation 
data and interpretation for design are needed. 

5. The uncertainties introduced can be overcome 
by design, given more-precise data out to the high 
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fluences required for the LMFBR application. Accom-. 
modation of swelling does, however, introduce an 
economic penalty that highlights the importance of 
metallurgical solutions to the swelling problem, 


THE SEVEN AGES OF FUEL 
(with apologies to W. S.) 


All the world’s a stage 
And all the fuels and materials merely players 
They have their exits and their entrances 
And one fuel element in its life plays many parts 
His acts being seven ages. At first the powder 
Mewling and turning in the mixer’s arms 
And then the whining pellet from the press 
With shining morning surface creeping like snail 
Unwillingly to furnace. And then the sinter 
Leaving the furnace with a black countenance 
Made to his mistress’ eyebrow. Then a 
Fuel element full of strange fixtures and 
Clad like the pard, 
Jealous in honour, sudden and quick to 
fracture 
Seeking the bubble reputation 
Even in the reactor’s mouth. And then the 
reactor fuel 
In fair round belly with good burnup lined 
Full of fissions and low cycle costs 
And so he plays his part. The sixth age shifts 
Into the lean and burned up stage 
With fixtures dangling and pouch on side 
His youthful clad well saved, a world too wide. 
For his shrunk fuel; and his big enriched voice 
Turning again towards natural treble, pipes 
And whistles in his sound. Last scene of all 
That ends this strange eventful history 
Is processing and mere oblivion 
Sans grid, sans clad, sans shape, sans everything. 
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ENDF/B-2: Is It Adequate for the 
Design and Evaluation Needs 
of the LMFBR Program? 


By Harry Alter* 


Abstract: /ntegral parameters for a@ series of fast reactor 
bench-mark assemblies, covering a wide range of energy 
spectra, have been calculated with the reference cross-section 
set ENDF/B-2. Multigroup constants, relative to each assembly, 
were generated and used in Sy transport and/or diffusion- 
theory calculations to obtain the parameters. These included 
kere, spectral indices (activation ratios), and central reactivity 
worths. Results of the several calculations indicate that the 
neutron cross-section data for such fuel-material isotopes as 
235U, 238U and ?3° Pu on the reference set ENDF/B-2 do not 
give satisfactory agreement with measured values of the 
parameters. 


Basic to the design and evaluation of LMFBR concepts 
and proposals are the nuclear data and calculational 
techniques used in the analyses. Valid comparisons 
between a variety of proposals and concepts depend 
largely on minimizing differences in analysis methods 
and providing reference sets of neutron cross-section 
data which may be used as standards. Accurate 
predictions of fast breeder reactor design parameters, 
such as critical mass, material reactivity worths, and 
spectral indices, require the development, maintenance, 
and evaluation of up-to-date sets of basic microscopic 
nuclear data files. The implications of data uncer- 
tainties on the predictions of fast breeder reactor 
design parameters have been reported in detail by 
Greebler et al.' 

The term “data testing,” as used here, is meant to 
represent the comparison of important integral parame- 
ters measured in a reactor environment with calculated 





*Atomics International, a division of North American 
Rockwell Corporation, Canoga Park, Calif. 91304. 


values using the reference data set ENDF/B-2. Correla- 
tion of calculated quantities against measured values 
should demonstrate how well the reference data set, 
together with the processing codes and computational 
procedures, predicts significant reactor parameters. 
Although results of such a testing program will seldom 
demonstrate an exacting relation between specific 
microscopic data and integral values, a systematic 
selection of test cases can show which material data 
consistently give poor correlation of computed and 
measured integral values. Analyses of such correlations 
can identify nuclides in need of reevaluation, perhaps 
with some indications of reaction type and energ 
range. Thus the testing indicates where the suspect data 
are in need of closer examination and perhaps revision. 


FAST REACTOR BENCH-MARK 
SELECTION AND SPECIFICATION 


Testing of the ENDF/B-2 was performed by the 
Cross Section Evaluation Working Group (CSEWG) 
against a set of 10 integral experiments. These experi- 
ments were chosen on the basis of system simplicity 
(insofar as isolating the three principal isotopes, 
235,238) and 739Pu), extent of experiments carried 
out, and an adjudged precision of experimental results. 
The choices, which emphasized three spectral regions 
(hard, medium, and soft) for each of the isotopes, con- 
sisted of four plutonium-fueled criticals, three high- 
enrichment uranium-fueled criticals, and three low- 
enrichment uranium-fueled criticals and are listed in 
Table 1. 
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Table 1 Fast Reactor Bench Marks by the 
CSEWG for Integral Testing of the ENDF/B 





Bench 
mark 
No. Bench-mark assembly 





l Jezebel?”® (bare plutonium sphere) 

2 Vera-11A°~'°® (plutonium-plus-graphite assembly) 
3 ZPR-3-48' ' (Pu-U—C-—Na system) 

4 Zebra-3°"!? (9:1 U/Pu metal assembly) 

5 Godiva?~® ** (bare enriched-uranium sphere) 


6 Vera-1B'* (enriched uranium-plus-graphite system) 
7 ZPR-3-6F°’'* (dilute 1:1 uranium system) 

8 ZPR-3-11' 5"! (7:1 uranium-metal system) 

9 ZPR-3-12'4 (4:1 uranium-plus-graphite system) 
10 Zebra-2' ®?° (6:1 uranium-plus-graphite system) 





In establishing bench-mark specifications, the 
CSEWG included adjustments for differences between 
the actual experimental assembly systems and the 
simplified models. These adjustments involve correc- 
tions for heterogeneity effects, worths of interface 
gaps, edge irregularities, and the cylinder-to-sphere 


Harry Alter 


studies. The bench-mark models thus describe the- 
experimentally determined critical dimensions for the 
homogeneous spheres with the quoted compositions to 
within the specified errors on the experimental eigen- 
value. 

Spectral indices measured at the core centers of the 
bench-mark source experiments are given for each 
assembly, with the values taken from most recently 
published evaluations of the measurements and ad- 
justed with known corrections for effects, such as flux 
depression by fission-chamber walls, and recalibrations. 
Measured values of central reactivity coefficients are 
also given, having been adjusted for sample size effects 
with empirical or theoretical corrections. For many of 
the material worths, large relative errors were assigned 
to account for unknown shielding or multiple scatter- 
ing effects in the large samples. In addition to 
references to the data sources,”-*° discussions of the 
errors On material worths and spectral indices are 
included with each case description.” ' 

Initial specifications for each of the 10 bench-mark 
assemblies were given by Davey and Hess,?' later 


modified where needed by Alter,?? and are sum- 
marized in Tables 2 and 3. 


volume shape factors. The adjustments were adopted 
on the basis of experimental data and/or semiempirical 


Table 2 Fast Reactor Bench-Mark Core Specifications (107? atoms/cm? ) 








Isotope/ 

nuclide Jezebel Vera-11A  ZPR-3-48 Zebra-3 Godiva Vera-1B ZPR-3-6F ZPR-3-11 ZPR-3-12 Zebra-2 

233 Py 3.7050 0.7213 0.1645 0.3465 

aiid 1 0.1751 0.0370 0.0107 0.0183 

241 Py 0.0117 0.0028 0.0011 0.0016 

334 0.0492 0.0092 0.0069 0.0046 0.0046 

2351 0.0016 0.0228 4.5000 0.7349 0.6727 0.4567 0.4516 0.2523 

Ses) 0.7405 3.1560 0.2498 0.0469* 0.75764 3.4392+¢ 1.69488 1.5550 

C 4.6200 2.0770 5.7540 2.6762 3.7950 

Na 0.6231 

Cr 0.1580 0.2531 0.0818 0.1635 0.1918 0.1486 0.1419 0.1008 

Fe 0.6090 1.0180 0.4559 0.6283 OF 712 0.5681 0.5704 0.3856 

Ni 0.0660 0.1119 0.0321 0.0689 0.0839 0.0718 0.0621 0.0397 

Al 0.02334 1.9019 0.0069 

Mn 0.0106 0.0080 0.0208 0.0059 

Cu 0.8050** 0.4794++ 

Mo 0.0206 

H 0.0058 

Radius,tt 6.385 13.985 46.46 23.67 8.741 19.138 22.995 31.61 28.76 45.593 
cm 





*Includes 0.0014 for ?77°U. 
+Includes 0.0029 for 73°U. 
tIncludes 0.0019 for 77°U. 
§ Includes 0.0019 for 77°U. 
q Includes 0.0124 of silicon in steel. 
**Includes 0.0450 for gallium and 0.0010 for aluminum. 
+7 Includes 0.0220 for gallium. 
£+Two-region spherical model based on as-built critical configuration plus adjustment. 
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Table 3 Fast Reactor Bench-Mark Blanket Specifications 
(10? atoms/cm?) 








Isotope/ 

nuclide Vera-L1A  ZPR-3-48 Zebra-3 Vera-IB ZPR-3-6F ZPR-3-11 ZPR-3-12 Zebra-2 
235 0.0250 0.0083 0.0301 0.0250 0.0089 0.0089 0.0089 0.0301 
cased | 3.4400 3.9976 4.0990 3.4400 4.0026 4.0025 4.0026 4.0990 
Cr 0.1680 0.1196 0.0913 0.1682 0.1129 0.1196 0.1237 0.0913 
Fe 0.6460 0.4925 0.3496 0.6464 0.4539 0.4925 0.4971 0.3496 
Ni 0.0710 0.0536 0.0360 0.0708 0.0494 0.0536 0.0541 0.0360 
Al 0.0060 0.1359 0.0060 
Mn 0.0051 0.0047 0.0111 0.0052 

Blanket thick- 
ness, cm 43.0 30.0 30.50 30.50 30.0 30.50 31.70 





FAST-NEUTRON DATA TESTING 
PARTICIPANTS 


The Data Testing Subcommittee of the CSEWG 
consists of scientists with interests in the areas of 
thermal data testing, fast data testing, and shielding 
data testing. Table 4 presents those subcommittee 
members, and their affiliations, involved in the testing 
of fast-neutron cross-section data. 


MULTIGROUP LIBRARY PREPARATION 
AND CALCULATION PROCEDURES 


The microscopic neutron cross-section data on the 
reference ENDF/B data set?* consist of resonance 
parameters, point data at many energy points, angular 
distribution data, and secondary energy-distribution 
data. To calculate the required parameters for the 
bench-mark assemblies, one must obtain cross sections 
for a limited number of energy groups. Since many 
materials have cross sections that vary rapidly with 
energy, knowledge of the flux distribution within a 
group is necessary to produce realistically weighted 
group constants. The computer codes used by the 
various participants to accomplish this task are sum- 
marized in Table 5. 

The computer codes used by the data testing 
participants to calculate Keg, spectral indices, and 
material central reactivity worths are summarized in 
Table 6. 

The general calculational procedure was to use a 
one-dimensional transport-theory code to calculate the 
eigenvalue, spectral indices, flux, and adjoint for each 
of the given assemblies. The flux and adjoint were then 
fed into a one-dimensional transport-perturbation- 
theory code for computing the material central reac- 
tivity worths. 


Table 4 Data Testing Subcommittee 
of the CSEWG 





Member Affiliation 





Argonne National Laboratory (ANLY 
Atomics International (Al) 


E. M. Pennington 

H. Alter and 
R.S. Hubner 

B. A. Hutchins and 
C. L. Cowan 

D. R. Mathews and 
J. M. Neill 

R. A. Grimesey and 
A.W. Brown 

R. J. LaBauve and Los Alamos Scientific Laboratory 
M.E. Battat (LASL) 

J. D. Jenkins Oak Ridge National Laboratory 
(ORNL) 

Westinghouse (WADCO) 


General Electric, Sunnyvale, 
Calif. (GE) 
Gulf Radiation Technology (GRT) 


Idaho Nuclear Corporation (INC)* 


W. W. Little, R. E. 
Schenter, and 
L. D. O'Dell 

T. A. Pitterle and 
N.C. Paik 


Westinghouse (WARD) 





*Presently at Aerojet Nuclear Corporation (ANC). 


Standard specifications for the calculations were 
provided and summarized in Tables 2, 3, and 7. 


RESULTS OF BENCH-MARK 
CALCULATIONS 


The values of Kerp calculated for the series of 
bench-mark assemblies by the data testing participants 
are summarized in Table 8. Tables 9 to 12 give calcu- 
lated results for spectral indices for four assemblies. 
The results in these tables are representative of results 
reported for all 10 assemblies in terms of the range of 
spectra and fuel isotopes covered; i.e., Jezebel has a 
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Table 5 Codes To Generate Multigroup Table 6 Codes To Calculate Bench-Mark 
Constants Assembly Parameters 

Participant Computer codes Participant Computer codes 

Al ETOE,** MC2 (Ref. 25) Al ANISN?* and COMPOZ?® 
ANL ETOE, MC2 ANL ARC ID TRANSPORT?” 

GE ENDRUN,?° TDOWN?’ GE SNID?® 

GRT GAND2?® and GGC5?° GRT IDF*? 

INC ETOP?° and PHROG?! INC SCAMP?° and MONA*! 
LASL ETOE, MC2 LASL DIF 1V** and DAC 1** 
ORNL SUPERTOG?? ORNI XSDRN*? and CITATION? 5 
WADCO ETOX?? and IDX** WADCO DTF IV and PERT V*° 
WARD ETOX and 1DX WARD ANISN and PERT V 








Table 7 Bench-Mark Assembly Calculational Specifications 





No. of intervals 








Calculation 

Assembly code Sy Core mesh Reflector mesh 
Jezebel ID TRANSPORT 16 40 

Vera-11A = ID TRANSPORT 8 40 40 
ZPR-3-48 IDTRANSPORT  4or8 30 20 
Zebra-3 ID TRANSPORT 8 40 30 
Godiva ID TRANSPORT 16 40 

Vera-1B ID TRANSPORT 8 40 40 
ZPR-3-6F ID TRANSPORT 8 40 30 
ZPR-3-11 ID TRANSPORT 4o0r8 30 20 
ZPR-3-12 1D TRANSPORT 5 40 30 
Zebra-2 ID TRANSPORT 4 40 30 





Table 8 Summary of Calculated ke¢¢ for CSEWG Bench-Mark Assemblies 











Laboratory 
Assembly Al ANL GE GRT INC LASL ORNL WADCO WARD Experimental 
Jezebel 0.9894 0.9902 0.9817 0.9880 1.000 + 0.003 
Vera-11A 0.9687 0.9682 0.9685 1.000 + 0.004 
ZPR-3-48 0.9681 0.9668 0.9688 0.9725 0.9880 0.9654 0.9600 0.9649 0.9650 1.000 + 0.004 
Zebra-3 0.9540 0.9480 0.9497 0.9494 0.9450 0.9430 1.000 + 0.003 
Godiva 0.9984 0.9991 0.9962 1.000 + 0.003 
Vera-1B 0.9893* 0.9974* 0.9897+ 1.0000 + 0.0028 
ZPR-3-6F 0.9916 0.9937 0.9910 0.9910 1.0000 + 0.0015 
ZPR-3-11 0.9640 0.9660 0.9724 0.9604 1.0000 + 0.0025 
ZPR-3-12 0.9830 0.9853 0.9763 0.9803 1.0000 + 0.0023 
Zebra-2 0.9600 0.9779 0.9593 1.000 + 0.004 





*S, calculation. 
+S, calculation. 
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hard flux spectrum and is 7*°Pu fueled; ZPR-3-48 has 
a medium flux spectrum and is 77° Pu—??U fueled; 
Godiva has a hard flux spectrum and is principally 
*35U fueled; ZPR-3-12 has a medium flux spectrum 
and is 735:?38U fueled. Tables 13 to 16 give calculated 
results for central reactivity worths for the same 
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In the tables of spectral indices and central 
reactivity worths, the results as reported by the 
participants have been averaged, and the ratio of the 
average calculated parameter to the measured value 
(when a measured value was available) is also given. All 
participant results were factored into the averaging 


representative set of assemblies. process. 


Table 9 Spectral Indices at Jezebel Core Center 
from S, 4 Transport Codes, 0,./a¢ (7°* U) 











Laboratory Average 
Oy LASL ORNL WADCO _ Experimental  calc./exp. 
of? *° Pu) 1.38 1.38 1.37 1.49 + 0.03 0.924 
of(?*7Np)* 0.94 0.99 + 0.05 0.944 
ofc "UP 0.184 0.177 0.183 0.205 + 0.008 0.885 
of(?>3U)* 1.56 1.56 1.61 + 0.10 0.971 
of? *° Pu) 0.868 0.863 
de? *°U)* 0.0818 0.0687 





*Infinitely dilute calculation. 


Table 10 Spectral Indices at ZPR-3-48 Core Center, 0/07 "et 





Laboratory 








Se en ee Average 
Ox ANL Al GE GGA INC LASL ORNL WADCO WARD _ Experimental calc./exp. 
of? 4° Pu) 0.2286 0.229 0.226 0.225 0.226 0.228 0.226 0.226 0.227 0.243 + 0.002 0.934 
of??? Pu) 0.9164 0.916 0.907 0.920 0.915 0.915 0.926 0.921 0.921 0.976 + 0.010 0.940 
of(??*U) 0.02980 0.0301 0.0297 0.0293 0.0293 0.0298 0.0283 0.0297 0.0300 0.0307 + 0.0004 0.963 
Gel *U) 0.141 0.139 0.139 0.132 0.137 0.140 0.146 0.140 0.139 0.138 + 0.004 1.009 
oe" >) 1.440 1.440 1.442 1.435 1.447 1.447 1.480 + 0.015 0.974 
of(?*" Pu) 1.269 1.262* 1.2692 
of(23 7Np)* 0.217 0.216 
aC" °U)* 0.0680 0.0675 0.0664 
af > Uy 0.197 0.197 0.196 
o¢(? 4" Pu) 0.214 0.210* 0.211 
o¢(7*°Pu) 0.2043 0.2057 0.2041 
o-(7°° Pu) 0.220 0.219 0.235* 0.218 0.223 
o¢(237 Np)* 0.8861 0.8705 
ag?" U)* 0.278 0.2717 
o-(?*5U) 0.2799 0.2786 0.2760* 0.2750 
op(??*U)* 0.2895 0.2817 
ae UO) 0.1647 0.1611 
o¢(Au)* 0.192 0.337 
o¢(! °B)? #32 1.34 1.32 
6¢(Mo) 0.0693 0.115* 0.0704 
o¢(Ni) 0.0196 0.0186* 0.0188 
oc(Ke) 0.0070 0.0077* 0.0076 
o¢(Mn) 0.0169 0.0328 
o(Cr) 0.0077 0.0083* 0.0083 
o(Al) 0.0016 0.0015* 0.0015 
o¢(Na) 0.0010 0.0012* 0.0009 
o¢(' ©O) 0.0006 0.0006 0.0005 





*Infinitely dilute calculation. 
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Table 11 Spectral Indices at Godiva Core Center from 
Si.¢ Transport Codes 





Laboratory 








Average 
Ox ANL LASL ORNL Experimental calc./exp. 
Indices as oxlof ? >=") 
of??? Pu)* 1.36 1.36 37 1.42 + 0.02 0.960 
of(?*U) 0:55 ‘O155 0.160 0.156 + 0.005 1.004 
of(??3U)* 1.57 1,57 1.63 + 0.10 0.962 
o¢(Au)* 0.105 0.104 0.105 + 0.010 0.994 
Indices redefined as oxlof ? >S 1) 

of(?**U) 4.9 4.8 4.9+ 0.2 0.990 
o¢(?3* U) 0.56 0.55 0.53 0.47 + 0.02 1.160 
of? > 26h) * 0.244 0.243 0.234 + 0.005 1.040 





*Infinitely dilute calculation. 


Table 12 Spectral Indices at ZPR-3-12 Core Center 
from S, Transport Codes, 0,./07 (7?° U) 

















Laboratory 
re ee SURE LA eee oe es Average 
Ox Al GE ORNL WARD Experimental  calc./exp. 
of???Pu)* — 1.07 1.07. 1.07 1.07 1.12 + 0.02 0.955 
of(?** U) 0.043 0.044 0.044 0.043 0.047 + 0.002 0.929 
oc(?? SU) 0.127 0.127 0.134 0.129 0.123 + 0.005 1.051 
of(??4U) 0.301 0.306 0.303 0.305 + 0.012 0.992 
of(??°U)* Lo 1.50 1.48 + 0.03 1.018 
of(?*° Pu)* 0.330 0.335 07333 0.326 
GAC eu) * 0.138 0.138 0.135 
o¢(??*U) 0.219 0.218 
o(?94U) 0.231 0.230 
o¢(Mo)* 0.0460 0.0461 
o-(Ni) 0.0217 0.0216 
o¢(Fe) 0.0068 0.0071 
o¢(Cr) 0.0060 0.0062 
o¢(Al) 0.0019 0.0019 
aC) 0.0001 0.0001 
*Infinitely dilute calculation. 
EVALUATION OF CALCULATION Assembly Pusicipont 
PROCEDURES 
Jezebel ANL, LASL, ORNL, WADCO 
All bench-mark assemblies were calculated by a Vera-11A GE, LASL, WADCO 
minimum of three data testing participants, with the ZPR-3-48 ALL 
distribution of assemblies calculated by the partici- Zebra-3 Al, ANL, GE, INC, WADCO, WARD 
pants shown at the right. ane ANL, LASL, ORNL 
Vera-1B ANL, INC, WADCO 
It would have been more meaningful if all the ZPR-3-6F = ANL, LASL, ORNL, WARD 
participants had calculated all the bench-mark assem- ZPR-3-11 Al, GE, INC, WARD 
blies. ZPR-3-48 was calculated by all participants. The ZPR-3-12 AI, GE, ORNL, WARD 
Zebra-2 AI, INC, WARD 


wealth of results reported permitted a more complete 
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Table 13 Material Worths at Jezebel Core Center and systematic analysis. Thus, for example, the central 
(10° Ak/k/mole) reactivity worths reported by ORNL (Table 14) are 
immediately suspect. Since all participants used the 
Laboratory A 5 y P ; P P 
nner same reference cross-section set, the ORNL results 
Material LASL ORNL WADCO Experimental  calc./exp. ; oa agt . 
point to significant differences between the calcula- 
hd 3 1813 1436 1847 1972 + 190 0.861 tional procedures used by them and by the other data 
*e?Pu_—-- 3098 2438 = 3098 3025 + 30 0.951 testing participants. The reactivity-worth differences 
237Np* 1145 1499 + 40 0.764 
238 [ys 156 97.2 162 21748 0.638 are therefore most probably due to the ORNL pro- 
235U* 1623 1270.3 1628 1528 + 15 0.986 cedures for generating multigroup constants and/or the 
233 * 2204 2822 2582 + 20 0.973 ORNL codes that calculate the integral parameters. 
Th* ~131.7 -124+4 1.062 iia wee F i 
po -17718 oor 1.078 Several fundamentally different approaches were 
Ta* “4186. =2050 <—190%6 1.009 used by the data testing participants to generate 
—— ~451.2  -550.8 = —477 + 20 1.050 effective core-and-blanket multigroup cross sections for 
Mo* —85.3 -89.7 —84 23 1.042 calculating the fast reactor bench-mark assemblies. It 
Zr* —50.1 -~68 +2 0.737 aay 
; as r é - e nt Sonstan e 
real ai, ns tees ans was required that region dep ndent group constants b 
Fe* _39 37.0 37.3 4044 0.921 generated for the various core and _ blanket 
= —32.3 -29+1 1.114 compositions. 
Ti* —30.9 -~49+2 0.631 In the method used by the MC? and GGCS5 codes, 
Al 25.2 -245 -26.8+0.8 0.927 A PL Rg pana ; 
° i. 4h cone ee 0.798 the energy range is divided into many narrow groups 
N _85.0 “4341 1.977 (~1000 to 5000), and cross sections for each isotope 
Cc -12.3 -10.7 -13.1+ 0.4 0.878 are computed for each fine group using a flux- 
Be 29.9 29+1 1.031 weighting function that depends on the energy- 
ad = Bee dependent cross sections of all other isotopes in the 
a ms . — 2400 ‘ ° : 
mixture. A  zero-dimensional fundamental mode 
*Infinitely dilute calculation. calculation is then performed with the effective broad- 


Table 14 Material Worths at ZPR-3-48 Core Center (10°° Ak/k/mole) 





Laboratory 








: Average 
Material ANL* Al GE* LASL ORNL WADCO WARD Experimental  calc./exp. 
ted a) 23.0 22.1 21.8 22.0 19.2 20.3 22.2 19.4 + 5.0 1.109 
cate) 136.9 132.0 132.4 132.7 118.4 130.2 135.0 106.1 + 0.9 1.235 
tel) TS 7.4 7.6 LE. —7,1 —8.1 7.8 —5.9+ 0.1 1.282 
Stool) 106.2 102.3 104.5 101 91.5 99.9 104.4 79.3 + 0.9 1279 
Tay —27.2 Bo 35.6 -29.6 + 0.9 1.102 
Mo -6.38 —6.24 10.2 -6.22 -5.6 -4,71 4.8 —4.15 + 0.05 1.520 
Ni —1.76 —1,74 1.76 1.72 1.4 Las, 1.73 1.07 + 0.04 1.580 
Fe 0.96 0.97 0.97 0.96 0.8 1.00 —1.01 0.69 + 0.04 1.381 
Mn 2.12 3.17 mF 2.60 1.18 + 0.04 2.200 
Cr —1.01 —1,02 ~1.06 —1.01 0.8 —1.09 1.09 -0.49 + 0.04 2.064 
Al 0.52 —0.53 0.60 0.51 0.4 —0.61 -0.61  -—0.42 + 0.03 1.287 
Na 0.14 0.16 ~0,10 -0.15 0.0 —0.07 0.10 ~0.14 + 0.03 0.744 
O 0.22 —~0.34 0.2 0.32 0.36 —0.11 + 0.03 2.615 
C 0.20 —0,22 -0.26 -0.21 0.1 0.19 -0.23  -0.05 + 0.03 4.014 
sl 3 100.6 97.1 101.2 87.1 102 FES 89.1 + 0.9 1.095 
sect J 3 176.1 155.7 172.7 
42405 2.9 Za 
229 Ue 13.3 12 
ser EU 193.1 190.2 BP235 188.8 





*Diffusion theory. 
tInfinitely dilute calculation. 
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Table 15 Material Worths at Godiva Core Center 
(10° Ak/k/mole) 

, Laboratory Average 
Material ANL LASL ORNL Experimental calc./exp. 
ae Put 1107 910.1 1422 + 130 0.899 
289 Py* 1905 1899 1567 1881 + 13 0.952 
ssid: 145 143 98.9 160+ 2 0.805 
wiring.) 1004 996 810 983+7 0.953 
Th* 34,3 52.6 2 ae | 4.825 
Au* —63.6 PS 49+2 1.440 
ee -319.9 -365 +7 0.876 
Cu* —24.4 -12+1 2.033 
Ni* 34 39.0 29 +1 1.259 
Fe* See, 11.9 344.3 6.577 
Al Sot 33% 1:2 Ae a | 
; 4.3 15,8:4.1,3 0.272 
Be 33.1 48+2 0.690 





*Infinitely dilute calculation. 
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performs a group collapse over space and energy to - 
produce the standard-group structures library. 

In another approach, used by the ETOX—IDX, and 
ENDRUN-—TDOWN codes, few-group (~100) infinite 
dilution cross sections and self-shielding factors for 
each isotope are first generated for a number of 
discrete temperatures and d9’s (where Oo is the total 
cross section per atom of all other isotopes in the 
mixture). Effective multigroup cross sections are then 
computed by averaging these group constants in space 
and energy to obtain the standard-group structure. 

Still another approach was that used by the 
SUPERTOG—XSDRN codes of ORNL. SUPERTOG 
accepts the nuclear cross-section data specified by the 
ENDF/B. These data are averaged over each specified 
group width with a suitable flux-weighting function, 
1/E or that supplied by the user. The code then 


Table 16 Material Worths at ZPR-3-12 Core Center (10° Ak/k/mole) 











Laboratory 
— —_—— —_—_____— Average 
Material Al GE ORNL WARD _ Experimental  calc./exp. 

Put4s%7""Pa)* 232 233 192 249 244 +7 0.919 
sin) 2 -7.9 —2.6 ~9.0 —6.7+ 1.0 1.032 
gid 155 156 121 160 S723 0.942 
ied Ue 219 210 283 269 + 7 0.951 
Ta* 30 14 32 -31+1 0.813 
Mo* 7.3 8.9 4.3 LEC.S 0.969 
Nb? 2 a7 Ot] 1.145 
Ni 552 2.8 -2.5 3:3 21-205 1.092 
re 7 1.4 3 1.6 1.5 + 0.4 0.987 
Al —0.6 0.3 —0.2 -0.5 —0.4 + 0.1 0.989 
Cc 0.9 0.9 £2 0.9 +2.0'2'0/3 0.494 
cs ic 63.5 64.7 64.3 

222 eu* 240.4 241.0 197.8 

a 29.6 30.4 38.1 

Cr LD RS 1.6 7 

*Infinitely dilute calculation. 
group cross sections being obtained by flux weighting computes 100-group cross sections for input to 


the narrow-group cross sections. A consistent PI 
option was employed by MC? users to generate the 
core group constants while the ordinary P1 option was 
used in the blanket region. GGC5 employed the BI 
approximation to compute the group constants in a 
99-group structure prior to collapsing these data to the 
standard 26 groups. MC? collapses the ultrafine groups 
(~1500) and the fine groups (54) to the standard- 
group structure. 

The INC code ETOP processes the ENDF/B data 
into 68 groups for input to PHROG. PHROG in turn 
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XSDRN. Nuclides with resonance parameters are pro- 
cessed twice, once to yield infinite dilution cross 
sections and resonance parameters which that code 
then uses to perform appropriate self-shielding calcula- 
tions. Effective multigroup constants are then calcu- 
lated in a 100-group format for use with XSDRN, 
which in turn calculates spatially and energy averaged 
broad-group constants for the reference-group struc- 
ture. 

The treatment of the resolved resonance region in 
the several codes differs substantially. MC2, ETOX, and 
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Table 17 ENDF/B-2 Microscopic Group Constants for the 
ZPR-3-48 *?°Pu Fission Cross Section Normalized to ANL Data 





Energy-group 








Group lower ey 
index boundary, eV LASL Al ORNL GE INC 
1 6.0650E 06 1.0000 1.0000 1.0173 0.9847 0.9830 
2 3.6790E 06 1.0000 1.0000 1.0027 1.0017 1.0017 
3 2.23 10E 06 1.0000 1.0000 1.0019 0.9998 1.0010 
4 1.3530E 06 1.0001 1.0000 0.9984 1.0008 0.9999 
5 8.2080E 05 1.0001 1.0000 0.9991 1.0012 0.9995 
6 4.9790E 05 1.0000 1.0000 0.9975 1.0015 1.0006 
7 3.0200E 05 1.0000 1.0000 0.9999 1.0001 1.0000 
8 1.8320E 05 1.0000 1.0000 1.0000 1.0000 1.0000 
9 1.1110E 05 1.0000 1.0000 1.0000 1.0000 1.0000 
10 6.7380E 04 1.0000 1.0000 1.0001 1.0001 1.0000 
1] 4.0870E 04 1.0000 = 1.0000 1.0014 1.0005 1.0015 
t2 2.4790E 04 1.0001 1.0000 0.9986 0.9984 0.9991 
13 1.5030E 04 1.0007. 1.0000 1.0124 0.9421 0.9559 
14 9.1190E 03 0.9996 1.0000 1.0184 0.9514 0.9755 
15 5.5310E 03 0.9999 1.0000 1.0301 0.9577 0.9907 
16 3.3550E 03 0.9991 0.9999 1.0215 0.9746 1.0014 
17 2.0350E 03 0.9974 0.9981 1.0328 1.0096 0.9961 
18 1.2340E 03 0.9970 0.9983 1.4385 0.9919 1.0647 
19 7.4850E 02 0.9970 0.9992 1.1226 0.9695 1.0267 
20 4.5400E 02 0.9952 0.9992 1.2223 1.0269 1.0852 
21 2.7540E 02 0.9693 0.9763 1.2323 0.9996 1.0724 
22 1.6700E 02 1.0119 1.0165 0.8902 0.9171 0.8523 
23 1.0130E 02 1.0034 0.9993 0.9465 1.0312 0.8804 
24 6.1440E 01 0.9930 0.9997 1.0288 1.1327 0.9048 
25 3.7270E 01 0.9676 0.9998 0.7658 0.5758 0.5895 
26 2.2600E 01 0.9867 0.9996 0.9560 1.1351 0.7232 





ENDRUN used the narrow resonance (NR) approxima- 
tion. PHROG as well as the combined SUPERTOG 
XSDRN codes used a Nordheim integral treatment. 
GGCS used an “exact” integral slowing-down treat- 
ment with 9327 energy points below 3.355 keV. For 
?39Pu, GRT extended the resolved region to 
3.355 keV with a ladder approach. 

All data testing participants were requested to use a 
standard 26-group structure in generating the effective 
multigroup libraries. Two, WADCO and WARD, pro- 
duced multigroup libraries in a 29-group structure. 
Although their calculated values for keg, spectral 
indices, and central reactivity worths were meaningful 
for intercomparison with similar results calculated by 
the other participants, direct intercomparisons of 
calculated group constants and group fluxes were not. 
Thus the WADCO and WARD data were omitted from 
the following intercomparison of effective multigroup 
constants. 


Examples of differences in microscopic group 
constants calculated for the ZPR-3-48 assembly are 
demonstrated in Tables 17 to 19. Group values for the 
fission and absorption cross section for 77?P and for 
the transport cross section for 7**U, as calculated by 
the several participants, are intercompared. 

These group constants, calculated for ZPR-3-48, 
were chosen for two reasons: (1) the data presented are 
representative of the fuel-material multigroup con- 
stants calculated for the other bench-mark assemblies 
and (2) this information was calculated by all the data 
testing participants. The values presented in the tables 
were normalized to those calculated by ANL. In some 
instances the reasons for disagreement among the 
multigroup cross-section sets have been ascertained. 
For example, in the first few groups, the (n,2n) cross 
section is substracted from the absorption cross section 
by some laboratories (AI, GE, INC, and LASL) but not 
others (ANL and ORNL). ORNL has added the (n,2n) 
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Table 18 ENDF/B-2 Microscopic Group Constants for the ZPR-3-48 
?3°Pu Absorption Cross Section Normalized to ANL Data 





Energy-group 








Laboratory 
Group lower 
index boundary, eV LASL Al ORNL GE INC 
1 6.0650E 06 0.9574 0.9574 1.0173 0.9847 0.9448 
és 3.6790E 06 0.9997 0.9998 1.0027 1.0018 1.0015 
3 2.2310E 06 1.0000 1.0000 1.0020 0.9997 1.0010 
a 1.3530E 06 1.0000 1.0000 0.9985 1.0008 0.9999 
5 8.2080E 05 1.0001 1.0000 0.9993 1.0011 0.9995 
6 4.9790E 05 1.0000 1.0000 0.9982 1.0010 1.0004 
7 3.0200E 05 1.0000 1.0000 1.0003 0.9999 0.9999 
8 1.8320E 05 1.0000 1.0000 1.0003 1.0004 0.9997 
9 1.1110E 05 1.0000 1.0000 1.0000 1.0000 1.0000 
10 6.7380E 04 1.0000 1.0000 1.0004 1.0004 0.9998 
11 4.0870E 04 1.0001 1.0000 1.0032 1.0011 1.0033 
lp 2.4790E 04 1.0002 1.0000 0.9983 0.9972 0.9991 
13 1.5030E 04 1.0016 1.0000 1.0050 0.9513 0.9605 
14 9.1190E 03 0.9991 1.0000 1.0096 0.9618 0.9769 
15 5.5310E 03 1.0000 0.9999 1.0149 0.9691 0.9881 
16 3.3550E 03 0.9990 0.9999 1.0132 0.9806 0.9937 
17 2.0350E 03 0.9987 0.9987 1.0397 0.9630 0.9980 
18 1.2340E 03 0.9975 0.9988 1.0885 0.9759 1.0215 
19 7.4850E 02 0.9965 0.9989 1.1256 0.9695 1.0163 
20 4.5400E 02 0:9951 . 09991 1.2237 0.9882 1.0521 
21 2.7540E 02 0.9634 0.9717 1.2536 1.0268 1.0478 
22 1.6700E 02 1.0059 1.0118 0.8811 0.9066 0.8457 
23 1.0130E 02 0.9940 0.9992 0.9300 1.0156 0.8741 
24 6.1440E 01 0.9909 0.9997 1.0426 1.2010 0.9179 
25 3.7270E 01 0.9734 0.9997 0.8017 0.7219 0.5945 
26 2.2600E 01 0.9859 0.9996 0.9261 1.1088 0.7285 








cross section to the fission cross section. The reason tor 
disagreement among the cross sections in the resolved 
(~1 to 300 eV) and unresolved (~300 eV to 25 keV) 
energy ranges has not been fully determined; however, 
suspicion has centered on the various methods used by 
ie participants to self-shield and Doppler broaden the 
effective multigroup constants. This point needs 
further investigation. 

One of the more obvious disagreements is in the 
transport cross section for 72*U (Table 19). The 
ETOE—MC2 path gave considerably lower values for 
certain groups. An investigation showed that incorrect 
transport cross sections for 72*U have been calculated 
by MC2 using all versions of ENDF/B data tapes. The 
reason for these incorrect (too low) values was traced 
to the resolved resonance-region specifications. In this 
region the scattering cross section is to be calculated 
using single-level Breit—Wigner parameters; however, 
the File 3 smooth cross-section background is zero over 
much of the region. In the energy range 5.0 to 50.0 eV, 
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the background is nonzero. From 50.0 to 3970.0 eV, 
the background is zero. In the latter range the ETOE 
program calculates negative scattering cross sections 
over small intervals. In particular, the cross section is 
negative near the energies 65.0, 100.43, 1964.5, 
1971.0, and 2542.5 eV; in addition, it is much too 
small over wider intervals. This constitutes an error in 
the ENDF/B ?°*°U evaluation for the resolved reso- 
nance region. 

The correction to this error is simple. In this energy 
range the ENDF/B2 data file must provide the correct 
matching of background cross section to the scattering 
cross-section method specification (i.e., single-level or 
multilevel Breit—Wigner parameters) to yield the cor- 
rect evaluated cross section. Thus the original data 
evaluator must be responsible for specifying the 
method for calculating the cross sections as well as the 
required matching background. 

Although such disparity among the effective multi- 
group constants is hardly calculated to promote con- 
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Table 19 ENDF/B-2 Microscopic Group Constants for the ZPR-3-48 
*38U Transport Cross Section Normalized to ANL Data 





Energy-group 








- Laboratory 
Group lower 
index boundary, eV LASL Al ORNL GE INC 
| 6.0650E 06 1.0000 1.0000 1.0140 1.0078 1.0112 
2 3.6790E 06 0.9999 1.0000 1.0112 1.0068 1.0111 
3 2.2310E 06 0.9999 1.0000 1.0047 0.9990 1.0035 
+ 1.3530E 06 1.0000 1.0000 1.0061 0.9996 1.0058 
5 8.2080E 05 0.9999 1.0000 1.0230 1.0010 1.0231 
6 4.9790E 05 1.0001 1.0001 1.0079 0.9997 1.0047 
7 3.0200E 05 1.0003 1.0000 1.0059 1.0009 1.0047 
8 1.8320E 05 0.9999 1.0000 1.0084 1.0087 1.0044 
9 1.1110E 05 0.9997 1.0000 1.0016 1.0022 1.0023 
10 6.7380E 04 1.0002 1.0000 1.0022 1.0019 0.9991 
11 4.0870E 04 0.9969 0.9971 1.0229 1.0115 1.0202 
12 2.4790E 04 1.0001 1.0000 1.0850 1.0257 1.0647 
13 1.5030E 04 1.0013 0.9999 1.0662 0.9751 1.0463 
14 9.1190E 03 0.9988 1.0000 1.0464 0.9265 1.0291 
15 5.5310E 03 1.0000 0.9998 1.0376 0.8815 1.0193 
16 3.3550E 03 0.6977 0.6896 1.1124 0.9165 1.1198 
17 2.0350E 03 0.1503 0.1509 1.5112 1.2849 1.7011 
18 1.2340E 03 0.9059 0.9128 3.1931 2.1777 3.8526 
19 7.4850E 02 0.9580 0.9594 1.4920 0.9766 1.8824 
20 4.5400E 02 0.9235 0.9255 1.2458 0.8953 1.4385 
21 2.7540E 02 0.9903 0.9920 1.1233 0.9427 1.2768 
22 1.6700E 02 0.9943 1.0037 4.2815 1.3409 5.6072 
23 1.0130E 02 0.9917 0.9963 2.3207 0.9630 2.9181 
24 6.1440E 01 1.0360 1.0473 7.5577 2.7971 7.8016 
25 3.7270E O1 1.0167 0.9990 0.9310 0.9066 0.9456 
26 2.2600E O01 0.9860 0.9889 14.6590 0.9611 16.9873 





fidence that the integral parameters calculated by the 
data testing participants can be meaningfully related to 
each other, one must take note of the following: 
Generally discrepancies among the multigroup con- 
stants appear to be significant in the first two groups, 
i.e., for energies greater than ~4 MeV and in the groups 
below ~15 keV. For the 10 bench-mark assemblies 
calculated, the mean fission-energy ranges from a low 
of ~200 keV to a high of ~1.06 MeV. For ZPR-3-48, 
the softest spectrum assembly, ~17% of the core flux 
falls below 15 keV and ~2% of the core flux falls 
above 4 MeV. It is thus reasonable to conclude that 
differences among effective multigroup constants, 
produced by the several computer-code processing 
techniques, are indeed related to the calculated integral 
parameters, but not to a significant degree. It is the 
author’s contention that this conclusion is borne out 
by the majority of the results presented; that minor 
differences are acceptable in that the several results are 


generally consistent relative to each other; and that, 
where this is not the case, the cause is yet to be 
identified. 

The computer codes summarized in Table 6 were 
then used to calculate the required integral parameters. 
Here the standardization of calculation methods by the 
participant was quite good. All reported values for Keg¢ 
were calculated with one-dimensional transport-theory 
codes. The computer codes used were basically minor 
variants of one another. The reported values for the 
spectral indices and central reactivity worths were 
calculated using both transport-theory and diffusion- 
theory methods. For the small bare assemblies, Jezebel 
and Godiva, only transport-theory calculations were 
performed; for the softer ZPR-3-48 assembly, both 
transport and diffusion results were reported. 

Comparison of results for the ZPR-3-48 spectral 
indices calculated with both diffusion-theory 
(CAESAR) and transport-theory (ANISN) codes is 
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presented in Table 20 and reveals insignificant differ- 
ences. 

A similar comparison for central reactivity worths 
(Table 21) shows that, although absorber nuclides 
undergo little change, scattering nuclides do differ 


Table 20 Comparison of CAESAR and ANISN 
Calculations for ZPR-3-48 Using ENDF/B-2 
Cross Sections 





Spectral indices at core center, 0,/af (7° * U) 








Ox CAESAR ANISN Experimental 
of (7°*U) 0.02896 0.02883 0.0307 + 0.0004 
of (?°°Pu) = 0.9147 0.9144 0.976 + 0.010 
of (?*°Pu) —-0.2293 0.2287 0.243 + 0.002 
of(?*'Pu) 1.280 1.280 
Gc (775 U) 0.2832 0.2834 
oc (73 *®U) 0.1316 0.1316 0.138 + 0.004 
oc (7°9Pu) ~=—:0..2259 0.2261 
oc (74°Pu) =: 0.2104 0.2105 
oc (7**Pu) =: 0.2195 0.2196 
Og (Al) 0.001560 0.001560 
Og (Cr) 0.006837 0.006841 
Og (Fe) 0.006434 0.006435 
Gq (Ni) 0.01622 0.01619 
Og (Mo) 0.07065 0.07067 
gq (Na) 0.001502 0.001503 
Og (C) 0.00006954 0.00006950 
0q ('°B) 1.3381 1.3387 
Gq (O) 0.0006269 0.0006250 





Table 21 Comparison of PERT and COMPOZ 
Calculations for ZPR-3-48 Using ENDF/B-2 
Cross Sections 





Material worths at core center 
(10° Ak/k/mole) 








Material PERT COMPOZ Experimental 
gatea § 96.97 97.35 79.3 + 0.9 

iO 7.166 7.36 5.9 + 0.1 

stat 127.4 127.9 106.1 + 0.9 

arma 20.66 20.50 19.4 + 5.0 

shi 2 188.43 189.4 

Al 0.402 -0.430 0.42 + 0.03 
Cr 0.944 0.980 0.49 + 0.04 
Ke 0.687 0,724 0.69 + 0.04 
Ni 1.516 1.556 1.07 + 0.04 
Mo 6.49 6.60 4.15 + 0.05 
Na —0.254 ~0,.290 0.14 + 0.03 
& 0.1087 0.0807 0.05 + 0.03 
IB 103.8 -104.8 89.1 + 0.9 

O 0.0497 0.0772 0.11 + 0.03 
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appreciably when their worths are calculated by the 
different methods. 

There appears to be a systematic irregularity in the 
ORNL reported results for central reactivity worths. 
The reasons for these significant discrepancies have not 
yet been uncovered; however, these results are so far 
out of line with the results of the other participants 
that they were discarded for the final evaluation of 
results. 

Again, from the standpoint of developing a more 
meaningful basis on which to judge and evaluate the 
accuracy and validity of the reference ENDF/B data 
set, it would have been preferable if all the data testing 
participants had been able to calculate the parameters 
as specified. However, with the exception of the ORNL 
worth calculations, one may conclude that differences 
among the various methods used to calculate the 
integral parameters perturb the consistency among the 
calculated results to a lesser degree than do differences 
among the several effective multigroup cross-section 
sets. 


EVALUATION OF CALCULATED RESULTS 


The reactor calculations for kerf were performed 
with one-dimensional transport-theory codes using the 
S, order given in Table 7 for the several bench-mark 
assemblies. The Ker¢ from the S,, calculation should, in 
general, differ by less than ~0.002 from the S. 
eigenvalue. For Jezebel, this difference was ~0.001. 
The anisotropy of the scattering was represented by 
the transport approximation, which introduces an error 
in the calculated Kee. For Godiva, this error is on the 
order of —0.003. For larger assemblies, this error 
becomes negligible. Thus the error introduced in Kerf 
through use of the transport approximation is small. 

The spectral indices and central reaction rates for 
2351) 238 and 73° Pu fission and 77° U capture were 
calculated using both transport and diffusion theory. 
Although the homogeneous bench-mark models have 
been corrected, where necessary, for the influence of 
heterogeneity on Kerr, this effect on the calculated 
spectrum has not been taken into account. Calculations 
by Kallfelz et al.4”7 on ZPR-3-11 for the heterogeneous 
cell showed practically no spectrum change from the 
homogeneous case. Calculations by Meneghetti*® 
agreed for ZPR-3-11; however, he then demonstrated 
that, for several other assemblies, the inclusion of 
heterogeneity tends to reduce the cell-averaged fission 
ratio of 738 U/*35U. The reduction in this ratio tends 
to increase the discrepancy between calculation and 
experiment. Calculations by Pitterle et al.4° confirm 
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that the inclusion of heterogeneity by using cell- Table 22 Ratios to **° U oy of 
homogenized cross sections gives, in general, softer 73°Pu of, 7°*U of, and 7**U 0, 





spectra than homogeneous calculations, thus worsening 





the agreement between calculation and experiment. Assembly (Calculated value) Experiment —Cale./exp. 
Although these effects lessen the value of com- henatal 1.37 1.49 + 0.03 0.92 
paring measured spectral indices with results calculated 0.181 0.205 + 0.008 0.88 
with homogeneous models, such comparisons still Vera-11A 1.12 1.07 + 0.02 1.05 
provide indications of general trends of the adequacy pegrenet pail iiipepesin 1.08 
OE TORR ZPR-3-48 0.917 0.976 + 0.010 0.94 
oe ae 0.0294 0.0307 + 0.0004 —-0.96 
The central reactivity worths were calculated using 0.139 0.138 + 0.004 1.01 
both transport- and diffusion-theory methods. The Zebra-3 1.13 1.19 + 0.014 0.95 
very large disparity between the measured and calcu- na airs 0.0461 2 0.0008 pag 
lated sis seals ‘Caabia tai ful san tina Godiva 1.36 1.42 + 0.02 0.96 
ated worth values imits reir usefulness as an ic enti- 0.157 0.156 + 0.005 101 

fier of errors in the differential cross-section data. The 

reported results have been used chiefly to help identif scat —— apa ap 1 
_ — y P y 0.076 0.067 + 0,001 1.13 
differences in calculational methods. 0.128 0.131 + 0.006 0.98 
Results for the calculated eigenvalues are sum- ZPR-3-61 1.21 1.22 + 0.03 0.99 
marized in Table 8. The results per assembly show 0.071 0.078 + 0.002 0.91 
riebemehn consistence uh differences : 2 0.110 0.104 + 0.003 1.05 
asonable consistency with differences among the ZPR-2-1] 1.13 1.19 + 0.02 0.95 
calculated quantities attributable principally to the 0.034 0.038 + 0.001 0.89 
several methods used to generate effective multigroup 0.119 0.112 + 0.005 1.06 
sets. The calculated Aegp values for the plutonium- ZPR-3-12 1.07 1,12 + 0.02 0.96 
fueled assemblies (ZPR-3-48, Vera-11A, Zebra-3, and oon ers ae = 
; are 0.129 0.123 + 0.005 1.05 
Jezebel) are low on the average from ~1% for Jezebel Fetus. 0.946 0.987 + 0.01 0.96 
nae gee ve j ‘ebra-2 .946 -987 + 0. .96 
to >5% for Zebra-3. For the ?°* U-fueled assemblies, 0.0294 0.0320 + 0.0005 0.92 
the calculated Kerp values are low by <1% (Godiva) to 0.139 0.136 + 0.001 1.02 





>3% (Zebra-2). 

The calculational methods used in the data testing 
exercises are, to a large extent, those currently ac- 
cepted for liquid-metal-cooled fast breeder reactor 1.10 
(LMFBR) conceptual design studies. Although these 
methods do not represent perfection in their applica- 
tion, they do represent an accepted mechanism 
through which reference data sets will be evaluated and 
applied. It is thus quite reasonable to conclude that, 
for the currently accepted reactor-design methods, the 
onus for the discrepant Keg¢ values falls primarily on 
the data, and particularly with the *7°U, ?°8U, and 
239 Du files on the reference ENDF/B data set. 

In Table 22 the calculated ratios of of 7*° Pu, of 
?°5U, and o¢ 7°*U relative to of *°*U are compared 
with the measured values. Here the calculated values 
represent an “average” of values calculated by all the — 
data testing participants. The individually reported 
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, 0.90 
values are reasonably consistent, and the average value i oe oe ee oe oe oe ¢ 
ia i le Seale Raa Diial o> om £€§ f¢469 © o & 
is representative of the individual values. ude SSBeteete sé 
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In Fig. 1 the ratios of of 7*°Pu/or 7°°U plotted Su NN 8 > N NN N 
for the several bench-mark assemblies are generally 1 to CSEWG BENCH-MARK ASSEMBLIES 
8% low. The exceptions, Vera-11A and Vera-1B, are 
currently being reviewed prior to remeasuring several Fig. 1 Ratios of calculated to experimental values of of 
of the integral parameters.*° Current belief is that the 23 °Pu/of??5U. 
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1.10 





reported of ?°°Pu/o¢ 7°°U values are too high. The 
results here tend to confirm the keg¢ results for the 
plutonium-fueled assemblies, i.e., an increase in the 
239Dy fission cross section would tend to reduce the 
disagreement between the calculated and measured 
values of Kege and the *?°Pu fission ratio. In Fig. 2 the 
general trend of the ratios of of 7°? jos SORE) de 1.05 -— = 




















predominantly low. Again, the Vera assemblies are z 
exceptions, but the reported measured values for these > 
assemblies are currently suspect. The other exception S peneneet 
to the trend is Godiva, the small bare ?**5U-fueled & 
assembly. A decrease in the 775 U fission cross sections x 
is not acceptable since the 775 U-fueled assemblies give S 7 
lower-than-desired calculated values of Ker. Thus one 2 
would propose to increase the 77*U fission cross s —— 
section to produce better agreement with measured 3 
values of x¢ ?°*U/ap ?°°U. The increase in this cross < 
section would also help alleviate the low Kerr values oes — _ 
calculated for the majority of the test assemblies. 
1.15 
0.90 
00 00 Le = N oy 
1.10 -— oe . 4 8 % 3 < 
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Fig. 3 Ratios of calculated to experimental values of o¢ 
ronag A, 





In Fig. 3 the general trend of the ratios of o¢ 
8U/or *35UJ is predominantly high. Although there 
is no reported ratio result for Vera-11A, the Vera-1B 
result is assumed to be suspect. 
aupete In this case an increase in the 7°°U fission cross 
section should be matched by a decrease in the capture 


CALCULATED/EXPERIMENTAL 
° sa 
a 8 


cross section of 73*U. It is seen that, as the 73° U atom 
0.90 — saa = density decreases for assemblies ZPR-3-11 to ZPR-3-12 
— to ZPR-3-6F to Vera-1B and for a much smaller 


increase in the *2°U atom density, respectively, for 
these assemblies, the calculated Ker values increase 
0.85 from about 0.967 to 0.981 to 0.992 to 0.991. Thus it 











a = a : & : = 2 appears that a reasonable decrease in the capture cross 
N @ « ca o£ «£ & «& Fs section of **%U would tend to bring both the 
fe 8 7 oF SS SD calculated ep values and the ratio of o¢ 7°*U/o¢ 

CSEWG BENCH-MARK ASSEMBLIES 2351 into better agreement with experiment. 
The average calculated central reactivity worths 
Fig. 2 Ratios of calculated to experimental values of of relative to the experimental values are given in 
298 U/ap?35U, Table 23. The average calculated values do not include 
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the ORNL contributions. It is evident from a review of 
the ORNL results that some systematic error was 
introduced into their calculations. The remaining re- 
sults, as reported by the data testing participants, are 
self-consistent, and the average value is representative 
of the individual values. Again, because of the large 


Table 23 Ratios of Calculated-to 
Experimental Central Reactivity 
Worths of Fuel Isotopes 











(Calculations)/ 
experiment 
Assembly ctl Cecil) Mas! 
Jezebel 1.06 0.73 1.02 
Vera-11A 1.38 0.92 
ZPR-3-48 1.28 1.28 1.24 
Zebra-3 1.40 1.49 | Pe 9 
Godiva 1.02 0.90 1.01 
Vera-1B 1.07 E72 1.08 
ZPR-3-6F 1.19 1.02 1.10 
ZPR-3-11 1.16 1.20 1.06 
ZPR-3-12 1.00 P25 0.96 
Zebra-2 1.28 1.29 1.20 





differences between the calculated and measured 
values, it is difficult and perhaps not even meaningful 
to attempt to relate data differences by the compari- 
sons of this integral parameter. As was determined by 
many independent calculations for these test assem- 
blies as well as others, calculated values of central 
reactivity worths for the fuel isotopes are generally 
quite high when compared to measured values. The 
results of these calculations (Fig. 3), do nothing to 
dispel this observation. 


CONCLUSIONS 


The reference ENDF/B-2 set of cross sections 
produced values of kere which are considerably low, 
especially for the ?*°Pu-fueled assemblies. This sug- 
gests that the cross-section data for the important fuel 
materials (73°Pu, ?38U, and 73°U) are not adequate 
for the needs of the LMFBR program. 


Not all the differences between the experimental 
and calculated quantities can be blamed on inaccurate 
cross sections. The transformation of the actual hetero- 
geneous three-dimensional bench-mark assemblies to 
homogeneous spherical geometry can introduce some 
error; a major area of uncertainty associated with the 
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ZPR-3-48 assembly bench-mark specifications is the 
corrections to be made for the effects of heterogeneity, 
not only on critical dimensions but also on spectral 
indices and reactivity worths. The present state of the 
art of heterogeneity calculations on ZPR cores is such 
that a certain amount of subjectivity was involved in 
simplifying the model to fit current codes. As a result, 
there may be about a 0.5% uncertainty in the 
calculated Kerp of ZPR-3-48 due to heterogeneity 
models. 


Computational methods also contain defects. The 
use of a single fission spectrum for both the core and 
blanket regions and the approximate treatment of 
elastic slowing down by heavy materials in some of the 
data-processing codes introduce errors, the magnitude 
of which is yet to be determined. Significant dif- 
ferences in effective multigroup constants have been 
reported. The reasons for some of these differences are 
understood, however, additional effort is required to 
resolve the remaining differences. 


Changes in basic neutron cross-section data within 
the limits of experimental error can make rather large 
differences in integral quantities. Although the integral 
results must be used with caution since they depend on 
many energy-dependent cross sections (i.e., various 
cross-section changes could combine to give about the 
same changes in integral quantities), it is clear that it is 
not possible to compile sufficiently accurate data for 
reactor calculations by evaluating experimental dif- 
ferential data without reference to results of integral 
measurements. 


We conclude that integral bench-mark testing of 
the reference cross-section set ENDF/B-2 has yielded 
results that can be meaningfully interpreted and 
applied in determining the accuracy of the cross- 
section data for the fuel isotopes **°Pu, **°U, and 
235. 
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FULL-LENGTH EMERGENCY COOLING HEAT 

TRANSFER PROGRAM TEST PLAN 

Author: R. G. Bock, General Electric Company, Atomic Power 
Equipment Department, San Jose, Calif. (for Phillips 
Petroleum Co., Idaho Falls, Idaho). 

Abstract: The test program plan for the BWR portion of the 
Full-Length Emergency Cooling Heat Transfer (FLECHT) 
Program is presented. The project organization, estimated 
program costs, and a test program schedule are included. 
(auth) 

Available as GEAP-5676 from National Technical Information 
Service, U.S. Department of Commerce, Springfield, Va. 
22151. 


BWR STANDBY COOLING HEAT-TRANSFER 

PERFORMANCE UNDER SIMULATED 

LOSS-OF-COOLANT CONDITIONS BETWEEN 

15 AND 300 PSIA 

Authors: J. D. Duncan and J. E. Leonard, General Electric 
Company, Atomic Power Equipment Department, San Jose, 
Calif. (for Idaho Nuclear Corp., Idaho Falls, Idaho), 

Abstract: Tests were conducted to determine the heat-transfer 
effectiveness of the emergency cooling systems at pressures 
up to 300 psia. A full-scale, electrically heated, stainless-steel- 
clad heater bundle was used to simulate a fuel bundle. 
Bottom-flooding and top-spray transient tests were per- 
formed under loss-of-coolant conditions. The thermal re- 
sponse of the cladding is reported. Individual local heat- 
transfer coefficients were calculated as a function of 
location, time, and the various test parameters. The variation 
of the coefficients provides insight into the heat-transfer 
mechanisms significant during the emergency cooling phase 
of the postulated loss-of-coolant accident. Increasing the 
system pressure above atmospheric pressure improved the 
performance of both the bottom-flooding and the top-spray 
modes of emergency cooling. The time required for the 
cladding to cool to saturation temperature in the flooding 
mode and the maximum cladding temperatures in the spray 


*Compiled by the Reactor Technology Section of the AEC 
Technical Information Center, Oak Ridge, Tenn. 37830. 
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mode were both significantly reduced as pressure was 
increased to 200 psia. 18 references. (auth) 

Available as GEAP-13190 from National Technical Informa- 
tion Service, U.S. Department of Commerce, Springfield, 
V4; 22951. 


REMOVAL OF RADIOACTIVE PARTICLES 
BY SPRAYS 


Author: L. F. Parsly, Oak Ridge National Laboratory, Oak 
Ridge, Tenn. 

Abstract: The theory of particle removal by sprays is reviewed, 
and the results of a number of recent experiments by others 
and 15 experiments in the ORNL Nuclear Safety Pilot Plant 
are discussed. The experimental results show that particle 
removal is not significant when no steam is condensing on 
the spray drops but is greatly enhanced when condensation 
occurs, A major part of the removal by condensation is due 
to drop nucleation on the particles and the combined effects 
of diffusiophoresis and Stephan flow. Drop nucleation is a 
transient phenomenon and appears unpredictable, but the 
rate of removal due to the other phenomena can be predicted 
within a factor of 2 or better. 43 references, (auth) 

Available as ORNL-4671 from National Technical Information 
Service, U. S. Department of Commerce, Springfield, Va. 
22051. 


EMERGENCY COOLING IN BOILING-WATER 

REACTORS UNDER SIMULATED 

LOSS-OF-COOLANT CONDITIONS. 

BWR-FLECHT Final Report 

Authors: J. D. Duncan and J. E. Leonard, General Electric 
Company, Atomic Power Equipment Department, San Jose, 
Calif. (for Idaho Nuclear Corp., Idaho Falls, Idaho). 

Abstract: The Boiling Water Reactor Full-Length Emergency 
Cooling Heat Transfer (BWR-FLECHT) test program has 
been completed. The results of the program are summarized. 
Five Zircaloy-clad and five stainless-steel-clad full-scale BWR 
bundle mockups were tested under simulated loss-of-coolant 
conditions, Approximately 150 transient and steady-state 
tests were conducted from 1968 to 1970, The test program 


SELECTED RECENT OR RECENTLY RELEASED PUBLICATIONS 77 


has significantly increased the understanding of the heat- 
transfer performance during the emergency cooling phase of 
the postulated BWR loss-of-coolant accident. The test data 
and the analysis indicate that both the top-spray and 
bottom-flooding modes of emergency cooling will terminate 
the cladding temperature increase and minimize the con- 
sequences of the postulated loss-of-coolant accident. The 
models developed for spray cooling heat transfer allow a 
conservative prediction of the consequences of the accident. 
Specifically the local cladding temperature and local cladding 
strength can be calculated as a function of time during the 
accident. 17 references, (auth) 

Available as GEAP-13197 from National Technical Informa- 
tion Service, U.S. Department of Commerce, Springfield, 
Va. 22151. 


ANALYSIS OF AEROSOL BEHAVIOR FOR FFTF 

HYPOTHETICAL ACCIDENTS 

Authors: R. L. Koontz and L. Baurmash, Atomics Interna- 
tional, Canoga Park, Calif. 

Abstract: The distributions of the released mass of fuel, fission 
products, fuel cladding, and sodium coolant within the FTI 
containment volumes following several postulated accidents 
are estimated. The mass leaking from the inner and outer 
vessels is computed assuming various leak rates from each 
compartment of the containment system, 15 references. 
(auth) 

Available as Al-69-Memo-95 from National Technical Informa- 
tion Service, U.S. Department of Commerce, Springfield, 
Vas 22158; 


PHYSICS DESIGN OF THE GAS-COOLED FAST 

BREEDER REACTOR DEMONSTRATION PLANT 

Author: Bruno Pellaud, Gulf Oil Corporation, San Diego, Calif. 

Abstract: The reactor physics work related to the Gas Cooled 
Fast Reactor Demonstration Plant is reviewed. Special 
attention is devoted to control rods, fuel management, and 
blanket design. A complete description of studies related to 
static and burnup behavior, reactivity coefficients, and 
reactor kinetics is supplied. A major section on basic data 
gives special attention to the cross sections of 77° Pu,?*° Pu, 
and ?3%U. Methods for calculating neutron spectra and 
streaming effects are described. 72 references. (auth) 

Available as GA-10509 from National Technical Information 
Service, U.S. Department of Commerce, Springfield, Va. 
PLEYEL. 


MEASURED LOSS OF SELECTED FISSION 
PRODUCTS FROM HIGH-BURNUP MIXED-OXIDE 
FUEL IN A MINIATURE PUMPED SODIUM LOOP 


Authors: D. E. Plumlee and P. E. Novak, General Electric 
Company, Breeder Reactor Department, Sunnyvale, Calif. 

Abstract: A defected section of mixed-oxide fuel preirradiated 
to ~41,000 MWd/Te was exposed to ~1150°I flowing 
sodium for 18 hr in a miniature pumped loop in a hot cell. 
Fuel-pin diametral increase and the transfer of fuel and 
fission-product species from the fuel specimen to different 
portions of the loop were measured to provide an improved 
design basis for the accommodation of failed fuel in LMF BR 


systems. Posttest examination of the fuel specimen showed: 
an ~3% increase in. the average diameter of the fuel section, 
less than 1% loss of fuel from the specimen with uranium 
preferentially leached compared with plutonium, and sub- 
stantial loss of fuel-specimen fission products ranging from 
29 to 85% for the isotopes investigated. Posttest examination 
of the loop indicated plateout of fission products in the 
colder regions, with the exception of cesium, which remained 
in solution in the sodium. 8 references. (auth) 

Available as GEAP-13731 from National Technical Informa- 
tion Service, U.S. Department of Commerce, Springfield, 
Va. 22151. 


MULTINODE ANALYSIS OF B&W’s 2568-MW(t) 
NUCLEAR PLANTS DURING A LOSS-OF-COOLANT 
ACCIDENT 

Authors: C. E. Parks, B. M. Dunn, and W. J. Schermer, 
Babcock & Wilcox Company, Nuclear Power Generation 
Department, Lynchburg, Va. 

Abstract: Multinode analyses were conducted to evaluate 
several loss-of-coolant accidents for 2568-MW(t) pressurized- 
water reactors. Time-step and noding studies are presented, 
basic assumptions are listed, and the selection of base 
parameters is discussed. The multinode model was used to 
perform a core cooling analysis for a spectrum of break sizes. 
Curves showing the parameters of interest are presented. 13 
references. (auth) 

Available as BAW-10034 from National Technical Information 
Service, U.S. Department of Commerce, Springfield, Va. 
ZEUS: 


COMPILATION OF MOST RECENT 

EXPERIMENTAL RESULTS ON FUEL-PIN 

FAILURE AND FUEL—SODIUM INTERACTION 

Author: G. Kessler, Kernforschungszentrum, Karlsruhe, West 
Germany (translated for Argonne National Laboratory, 
Argonne, Ill., from Report PSB-389/1971). 

Abstract: The results of the TREAT experiment with series S, 
E, and H published up to March 1971 have been compiled 
and compared with the meltdown experiment of Atomics 
International. A critical analysis shows that the mechanical- 
energy-conversion factor is <0.2% in all experiments. A 
comparison of the particle-size distributions occurring during 
the fuel—sodium interaction shows that the mean particle 
size falls in the range 0.2 to 0.5 mm in all experiments. 
Similar results were also obtained in the SPERT experiments. 
rhe significance of these results for an estimate of the effects 
of severe hypothetical accidents in fast sodium-cooled 
reactors is discussed. 20 references. (auth) 

Available as ANL-trans-903 from National Technical Informa- 
tion Service, U.S. Department of Commerce, Springfield, 
Va. 22151. 


ECCSA AND MUCHA: COMPUTER CODES FOR 

THE ANALYSIS OF EMERGENCY-CORE-COOLING 

SYSTEMS. Topical Report, Task 18 

Author: Ronald A. Cudnik, Battelle Memorial Institute, 
Columbus Laboratories, Columbus, Ohio 
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Abstract: Several computer codes have been developed for the 
analysis of a reactor loss-of-coolant accident and subsequent 
emergency-core-cooling injection. Included are the ECCSA4 
code, which was developed to predict the thermal and 
hydraulic behavior of a single reactor fuel rod and its 
associated core channel; the MUCHA2 and MUCHAI codes, 
which are pressure-coupled, parallel core channel codes 
initially developed to investigate the possibility of preferen- 
tial flow due to thermally dissimilar channels; and the 
expanded version of the MUCHA2 code, which was devel- 
oped to provide the capability of analyzing the complete 
primary system with the core heat-transfer effects intimately 
coupled to the system analysis. The mathematical models 
that have been developed to describe the thermal and 
hydraulic processes are discussed in detail. These include 
models for heat conduction, heat generation, metal—water 
reaction, core fluid behavior, and plenum behavior. Source 
decks and sample problems for ECCSA4, MUCHA2, and 
MUCHAI will be available through the Argonne Code Center 
for the CDC-6400 computer. 32 references. (auth) 

Available as BMI-1916 from National Technical Information 
Service, U.S. Department of Commerce, Springfield, Va. 
22051. 
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INTERIM STRESS DESIGN CRITERIA FOR 
LMFBR VESSEL AND CORE STRUCTURES 


Authors: E. Y. W. Tsui, J. V. Marron, and A. T. Brown, 


General Electric Company, Breeder Reactor Department, 
Sunnyvale, Calif. (for Oak Ridge National Laboratory, Oak 
Ridge, Tenn.) 

Abstract: Interim stress criteria for LMFBR vessel and core 
structural design were developed. The components to which 
the criteria apply are classified into two groups: (1) long-life 
components and (2) core components. Examples of long-life 
components are: the reactor container, primary coolant 
piping, reactor closure assembly, core clamping mechanisms, 
etc. Examples of short-life components are: fuel-assembly 
hexagonal channels and internals, blanket assemblies, and 
control-assembly channels. The scope involves the establish- 
ment of basic stress limits and recommends appropriate 
accumulative damage criteria considering the effects of hold 
time (frequency) and mean strain on low-cycle fatigue at 
temperatures in the creep range. 24 references. (auth) 

Available as GEAP-13719 from National Technical Informa- 
tion Service, U.S. Department of Commerce, Springfield, 
Va, 22151. 
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Letters 


The following exchange of letters was prompted by the R. B. Briggs article, “Tritium in 
Molten-Salt Reactors,” which appeared in Vol. 14, No. 4, of this journal. 


Dear Mr. Briggs: 


We have found a flaw in your model depicting tritium transport in the molten-salt reactor 
system. Whenever a diatomic gas dissolves in a liquid as a monatomic species, Henry’s law does not 
hold. In your case, species in the liquid include TF. For this, Sievert’s law (pe = HsC) would be 
expected to describe equilibrium conditions better. Some justification for this is the fact that the 
hydrogen—liquid sodium system follows Sievert’s law down to very low concentrations. 

Substitution of Sievert’s law into the transport equations gives a flux which is directly 
proportional to tritium concentration in the molten salt. This change would provide better 
agreement between calculated and measured results, which you indicated on page 338 could be 
brought about instead by decreasing the permeability or raising the exponent on the diffusion 
equation. We think our change would be the more reasonable one. 


Paul Barton 
Visiting Scientist 


John Holmes 
Associate Chemical Engineer 


Argonne National Laboratory 
Idaho Facilities 


Gentlemen: 


Thank you for your letter of February 4 and your thoughtful interest in my paper on tritium 
in molten-salt reactors. The model used in the calculations does indeed have flaws, so I have given 
considerable thought to your suggestion that use of Henry’s law rather than Sievert’s law to 
describe the relationship between the partial pressure of hydrogen and the concentration in the salt 
is a flaw not previously recognized. I conclude that we should continue to use Henry’s law in the 
way that we use it. The space allotted to methods of calculation in my paper was not sufficient to 
be very informative. I hope that the following discussion will provide a better understanding of the 
basis for my conclusion, 

According to our chemists, chemical considerations alone lead to the expectation that 
hydrogen should dissolve in the 2Lif—BelF, fuel carrier and coolant salt of the MSRE as a 
diatomic rather than a monatomic species. This thesis is supported by recent, not-yet-published 
measurements of the solubility of hydrogen in very pure 2LiF—BeF,, salt at 600°C. The solubility 
was found to vary linearly with pressure over the range from | to 2 atm, and we see no reason for 
the relationship to change at lower pressure. The Henry’s law constant has about the same value as 
for neon, which, although a monatomic gas, has about the same molecular radius as hydrogen. This 
finding is consistent with the theory of solubility of unreactive gases in molten salt. Preliminary 
measurements at 700°C show that the solubility of hydrogen, like that of neon and the other noble 
gases, increases with temperature, but the temperature dependence does not seem to be as great. 

The solubility of hydrogen fluoride in 2Li—BeF, has also been measured* and again Henry’s 
law applies. Hydrogen fluoride, however, dissociates and forms complex ions in the salt, and the 
solubility decreases with increasing temperature. 





*MSR Program Quarterly Progress Report for Period Ending October 31, 1958, Report 
ORNL-2626, p. 85ff, Oak Ridge National Laboratory. 
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With pure 2Lil'—Bel’,, then, one can dissolve hydrogen in the salt, no significant amount of 
hydrogen fluoride is produced, and the concentration of hydrogen is proportional to the partial 
pressure, Also, one can dissolve hydrogen fluoride in the salt, no significant amount of hydrogen is 
produced, and the concentration of hydrogen fluoride-related species, which I will refer to simply 
as hydrogen fluoride, is proportional to the partial pressure. When both hydrogen and hydrogen 
fluoride are present, we relate their individual concentrations in the pure salt and their partial 
pressures, appropriately, we believe, by the Henry’s law constants that were measured separately. 
Phe fluid in the secondary or coolant system of the MSRE was 2LiF—BeF, , no hydrogen fluoride 
was added, so I find no alternative to the use of Henry’s law to relate the concentration of 
hydrogen in the salt to the partial pressure in that system. 

The fluid in the primary or fuel system of the MSRE was 2LiF—BeF, containing about 
5 mole % Zr, and 0,2 mole % UF, and UF. The zirconium fluoride should have little effect, but 
hydrogen reduces UI‘, to produce hydrogen fluoride by the reaction 


UP, + 7h, SUF, Fel (1) 


The equilibrium quotient for this reaction has been measured by Long and Blankenship* and is 
of the form 


_ CUR, x PHI 
Cur, x (PH, )% 


where the symbols are those defined in the nomenclature attached to this letter. 

In the MSRE the concentration of total uranium did not change significantly with time. 
Beryllium metal or other reducing agent was added periodically to counteract the oxidizing effect 
of fission and keep the concentration of UF, at the desired level. In our calculations, then, we 
relate the partial pressures of hydrogen and hydrogen fluoride by 


% | CUE (3) 
PH, = + ~=—  }) PHE - 
H, K Cur, 


Because we believe that hydrogen molecules and hydrogen fluoride species coexist in the salt, 
we relate their concentrations to the partial pressures by the Henry’s law constants to obtain 


Lf Gen 
1 UE ‘ 
(Hy Cy)? =—(-— ) (AnFCur) (4) 
H,*H, z(cor) 


Assuming the validity of the discussion above, let us go on with the calculation of the 
distribution of tritium in a reactor. Consider a simple case in which the tritium is produced 
uniformly in a fuel salt that circulates through a piping system at 650°C. Tritium escapes from the 
salt by diffusing through the pipe walls and by transfer to a gas stream that is equilibrated with the 
salt and discharged to an off-gas system, The rate of transport of tritium molecules to the pipe 


walls is 
Ow =hA (Co Cw) (S) 
Assuming that the concentration of tritium outside the pipe wall is zero,+ 


PA % _PA we 
Ow=— Pw=—> (Hy,Cw)” (6) 


The amount of tritium leaving in the off-gas, which is at 650°C and | atm, is 


OG = BGHy, Co (7) 


*G, Long and F. F. Blankenship, The Stability of Uranium Trifluoride. Part II, Stability in 
Molten Fluoride Solution, Report ORNL-TM-2065, Part II], Oak Ridge National Laboratory, 
November 1969, 

+Although we suggested, on the basis of Webb’s data (Report NAA-SR-10462) and many other 
references, that the exponent on the pressure might change from 2 to 1 at very low pressure, 
recent measurements of our own indicate that such a change is highly unlikely. 
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The salt contains uranium, so it also contains tritium fluoride in a concentration in the bulk 
salt that is related to Co by Eq. 4. Then 


K (Cu 
Cr=zs>— ca) (Ay Co)” (8) 
Hur \Cur, 2 


At this point, we introduce two assumptions that remain to be substantiated and are major 
faults if they are seriously in error. 

One way for the tritium in tritium fluoride to be made available to pass through the pipe walls 
is for it to be reduced by the metal of the wall in a reaction of the type 


M° + 2H* = M?* + 2H (9) 


We assume that these reactions are relatively slow and a negligible source of tritium for 
diffusion through the metal. 

The second assumption is related to the equilibrium between tritium, tritium fluoride, and the 
uranium species in the salt. We assume that equilibrium exists in the bulk of the salt but not in the 
salt film through which the tritium passes to the metal surface. The conversion of tritium fluoride 
to tritium at the metal surface would be controlled by the rate of reaction of tritium fluoride with 
UF, in salt near the surface, We have not considered the kinetics of the reaction yet and usually 
assume that the rate is too slow to be a source of tritium for diffusion through the metal. 

Because of these two assumptions, the tritium in tritium fluoride cannot pass through the pipe 
walls, but it is free to pass into the gas stream. The rate of release to the off-gas is 





3 


CUF 
OGFr = /, BGHyFCr = 16K (ce ) (Hy, Co)" (10) 


The sum of the removal rates must equal the production rates, so 
R=QwtQcGtQcGr (11) 


Equations 5, 6, 7, 10, and 11 are solved simultaneously to obtain a common value for Cog, and this 
value is used to calculate Cy and the fraction of the tritium that escapes to each of the sinks. 

Calculation of the distribution in the MSRE is more complicated because tritium was sorbed by 
the graphite in the reactor core, and tritium that passed through the walls of the heat-exchanger 
tubes dissolved in the coolant salt and escaped by several paths. Also, hydrogen was introduced 
into the primary system by decomposition of oil in the pump bowl. Nevertheless, the principles are 
the same, so I have calculated the flow through the pipe walls and to the off-gas for the simple 
model by use of the constants shown in the nomenclature, They are the same as those used for 
calculating the MSRE reference case, except that the surface area in the simple model was 
increased to account for the tritium transport to the secondary system in the MSRE and the 
tritium production rate was reduced to roughly compensate for sorption by the graphite. 

The equations, after substitution of the numerical values for the constants, are 


Ow = 2.9 x 103(Cg — Cw) = 2.4 x 10! ‘CH 
QG=5.7x10Co Ogr=13x 10°C? 
Cr=1.4 x 10°CR 

PH, =9.6x10°7'CH = pHER=0.2 x 10'*CuE 


Solving the equations gives 


Ow= 1.5 x 10'5 = 13 Ci/day 


Og = 3.0 x 10'5 = 25 Ci/day 

Ocr = 0.1 x 10'* = 1 Ci/day 
€n=5.25010"" Po=six 10° 
Cy = 0.4 x 108 Py=2.4'xi107° 
Cr=1.0x10'? — pp=2.9 x 107 
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On examination of the results, one sees immediately that the major barrier to the transport of 
tritium through the pipe walls is the salt film. The metal is transparent. All other conditions 
remaining constant, one must decrease the mass-transfer coefficient to the wall or greatly reduce 
the permeation coefficient of the metal if he wishes to reduce the flow through the metal. There is 
no basis for choice in this example, but a combination of the two appeared to give the best 
agreement between the calculated and the measured distributions in the much more complicated 
MSRE cases. Because of the many uncertainties, this merely tells us that we should test oxidized 
Hastelloy tubing to determine whether the permeability is greatly less than that of clean tubing. It 
also points up the need for good values of the mass-transfer coefficients. 

Although I have stated the conviction that hydrogen molecules and hydrogen fluoride species 
coexist in the salt, what if this were not true? Suppose that only hydrogen fluoride species exist in 
the fuel salt (in effect that molecular hydrogen is insoluble in the salt), but, in the presence of UF, 
and UF,, the hydrogen is released from the salt as an equilibrium mixture of hydrogen and 
hydrogen fluoride. What difference would this make in the results of the calculations? I find little 
difference for the MSRE conditions as they are represented by the example calculated above. 

First, we would have to assume that only tritium fluoride is produced when tritium is formed 
in the salt. The partial pressure of tritium fluoride would be 


Pe= Hyper (12) 


and that of tritiur would be obtained from Sievert’s law relationship 


Hr (Cur 
% 2p UF, : 
= c 13 
Po ra Cy ss) F (13) 


If the calculated concentration of tritium fluoride in the bulk salt, Cp, turns out to be about 
the same as we calculated for the example, then the assumption that equilibrium exists in the bulk 
salt and the mathematical relationships used would cause po and pf to have values equal to those 
we obtained from our calculations with tritium molecules in the salt. The rate of removal of 
tritium and tritium fluoride by the off-gas would be unchanged. 

The flow of tritium through the salt film to the pipe wall would be in the form of tritium 
fluoride species. In our sample calculation, Cf ~ 2Cg, but the number of tritium atoms per 
molecule is only half as great in the fluoride, For lack of better information, we might assume that 
the mass-transfer coefficients for tritium molecules and tritium fluoride species are about equal. 
Then the calculated rates of tritium flow to the wall will be the same whether transport is by 
tritium fluoride species or by tritium molecules, if their concentrations at the wall are much 
smaller than in the bulk salt, The concentration of tritium fluoride in salt at the wall that would 
make py = 2.4 x 10?° atm is 8.5 x 10° molecules/cm?. This is indeed small in comparison with 
the concentration in the bulk salt, so the flow through the pipe wall would remain unchanged. 

Although the results are the same for both models when UF ,/UF, = 0.01, they differ when 
UF,/UF, is larger or smaller. With larger UF,/UF, the “molecular hydrogen” transport model 
predicts little change in the distribution, The “hydrogen fluoride” transport model predicts that 
the flow through the metal wall should decrease and the flow to the off-gas should increase. With 
smaller UF ,/UF, the molecular hydrogen transport model predicts lower escape through the pipe 
walls, while the hydrogen fluoride transport model predicts the opposite effect. The effects are 
shown in Table | by the results of calculations for the simple reactor system. 

The model that best represents what I believe to be the actual physical situation is the one that 
allows for the coexistence of hydrogen molecules and hydrogen fluoride in the salt and transfer of 
hydrogen through salt films to metal surfaces by both. We probably will adopt such a model when 
we have better information about the reaction kinetics and the mass-transfer coefficients for the 
hydrogen fluoride species. Results of calculations based on this model and the assumption that 
chemical equilibrium exists every where in the salt are also shown in Table 1. 

On the basis of all that has gone above, I conclude: 

1. The model that was used in the calculations for my paper properly describes hydrogen 
transport in the secondary system of the MSRE and of the MSBRs considered there. 

2. For the present, the model adequately describes hydrogen transport in the primary system if 
the calculations show half or more of the hydrogen in the bulk salt to be in the form of molecular 
hydrogen. This is generally the case when UF, /UF, > 0.01, the reference condition for the MSRE 
and MSBR calculations, 

3. Changes in the mass-transfer and permeation coefficients would have been required to obtain 
best agreement between the measured and calculated tritium distributions in the MSRE regardless 
of the model used to describe the transport in the primary system. 
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4. Until assumptions are validated concerning rates of reaction of hydrogen fluoride species in 
the salt with the metal wall and with reduced metal ions in the salt film, distributions calculated by 
use of the model must be viewed with more than ordinary skepticism when most of the hydrogen 
in the bulk salt is in the form of hydrogen fluoride species. Transport through the metal walls tends 
to be underestimated. This is generally the condition when UF,/UF, < 0.001, a condition that 
was shown in my paper to be favorable for retaining tritium in an MSBR. 


Table 1 Effect of Model on Calculated 
Distribution of Tritium in a Simple Reactor System 





Tritium distribution 





Through a ——— 
Transport walls, as TF, 
model UF, /UF, %o %o % 
H, 107 35 65 0 
HF 4 96 0 
H, + HF 35 65 0 
H, 167 33 65 2 
HF 33 65 2 
H, + HF 50 48 2 
H, 10° 26 54 20 
HF 87 6.5 6.5 
H, + HF 87 6.5 6.5 
H, oe 5 11 84 
HF 89 0 11 
H, + HF 89 0 11 
H, 10° 0 0 100 
HI 63 0 37 
H, + HF 63 0 37 





Hopefully my explanation has not been too lengthy or raised more questions than it has 
answered. Please let me know if you have problems with it. I try to keep up with work on 
hydrogen and tritium in EBR-II and LMFBRs through the ANL Reactor Development Program 
progress reports and Sodium Technology quarterly reports but have not read of any calculations 
being done on the behavior of tritium in LMFBR systems, If you know of such work, I would 
appreciate hearing about it. 


NOMENCLATURE 


A = area of metal surface, cm? = 4.8 x 10° 
B = conversion factor, molecules/(cm* )(sec) atm = 7.9 x 1018 
Cy, = concentration of H, in salt, molecules H, /cm* 
CyF = concentration of HF in salt, molecules HF/cm? 
Co = concentration of tritium in bulk salt, molecules T, /em? 
Cr = concentration of tritium fluoride in bulk salt, molecules TF/cm? 
Cy = concentration of tritium in salt at wall, molecules T, /cm? 
Cy F,/Cur, = ratio of concentrations of uranium species in salt = 0.01 
G = gas flow rate, cm?/sec = 120 
h = mass-transfer coefficient for H, in salt, cm/sec = 0.006 
Hy, = Henry’s law constant for H, , atm/(molecule H, /cm*) = 0.6 x 107” 
Hyp = Henry’s law constant for HF, atm/(molecule HF/cm?) = 0.2 x 107 
K = equilibrium quotient for reduction of UF, by hydrogen, atm # = 1.1 x 10° 
P = permeation coefficient of metal, (molecules H, )(mm)/(cm? )(sec)(atm) =2x10'° 
PH, = partial pressure of H, , atm 
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PHF = Partial pressure of HF, atm 
Po = partial pressure of T, over bulk salt, atm 
PF = partial pressure of TF over bulk salt, atm 
Pw = Partial pressure of T, over salt at wall, atm 
Qg = flow of tritium to off-gas, molecules T, /sec 
Oar = flow of tritium fluoride to off-gas, molecules T, /sec 
Qw = flow of tritium through wall, molecules T, /sec 
R = rate of production of tritium, molecules/sec = 4.6 x 10’* 
t = metal thickness, mm = 10 


R. B. Briggs, 

Associate Director 

Molten-Salt Reactor Program 
Oak Ridge National Laboratory 





SYMPOSIUM ON THE MANAGEMENT OF RADIOACTIVE 
WASTES FROM FUEL REPROCESSING 


Paris, France, Nov. 27—Dec. 1, 1972 


The European Nuclear Energy Agency and the International Atomic Energy Agency are jointly 
sponsoring a symposium on the management of radioactive wastes. It will be held at OECD Headquarters, 
2, rue Andre Pascal, Paris. 

Topics to be covered include: sources of radioactive wastes, waste production and management, 
routine releases from reprocessing plants, storage of fission-product solutions, solidification processes for 
fission-product solutions, ultimate storage or disposal, decommissioning of waste-storage facilities and 
reprocessing plants, and economic aspects of waste management. 

As a part of the symposium, technical visits will be arranged (at the participants’ expense) to plants in 
France, Germany, and possibly other European countries. Participation will be limited to persons 
nominated by the government of a Member State of the ENEA, IAEA, or an international organization 
invited to participate. 

For details about submission of papers to be presented, registration, accommodations, etc., persons 
outside the United States should write 


ENEA/IAEA Symposium on the Management of 
Radioactive Wastes from Fuel Reprocessing 
c/o European Nuclear Energy Agency 
38, Boulevard Suchet 
F-75 Paris 16e, France 
and persons in the United States should write 
Mr. John H. Kane 
Special Assistant for Conferences 


U. S. Atomic Energy Commission 
Washington, D. C. 20545 
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A USAEC Publication of Special Interest 


Nuclear Reactors Built, Being Built, or Planned 
in the United States as of Dec. 31, 1971 


This compilation contains current information about facilities for domestic use or export which 
are capable of sustaining a nuclear chain reaction. Civilian, production, and military reactors are 
listed, as are reactors for export and critical assembly facilities. 40 pages, 8 x 10%, paperback. 


Available as TID-8200(25th Rev.) for $3.00 from 


National Technical Information Service 
U. S. Department of Commerce 
Springfield, Virginia 22151 





Recent USAEC publication of special interest 


A SELECTED, ANNOTATED BIBLIOGRAPHY OF THE 
CIVIL, INDUSTRIAL, AND SCIENTIFIC USES 
FOR NUCLEAR EXPLOSIONS 


Compiled by ROBERT G. WEST, Division of Peaceful Nuclear Explosives 
and 
ROBERT C. KELLY, Division of Technical Information Extension 
U.S. Atomic Energy Commission 


This bibliography contains 1594 selected references to reports and published literature on the 
USAEC’s Plowshare Program which were published since the inception of the program in 1957 to 
Dec. 31, 1970. The Plowshare Program was begun to study the feasibility of using nuclear explosions 
for peaceful purposes. Studies have been conducted on many potential applications in such areas as 
excavation of harbors, canals, and mountain passes; production of transplutonium isotopes; neutron 
physics research; gas and oil recovery and storage; waste disposal; water resources management; and 
geothermal power production. 


Author, experiment, foreign literature, and report number and availability indexes are included. (8 by 
10¥,, 359 pages) 


Available as TID-3522(9th Rev.) for $6.00 from: 
National Technical Information Service 
U.S. Department of Commerce 
Springfield, Virginia 22151 
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nuclear-explosion seismology Chitioal 


Howard C. Rodean 
Seismic/Evasion Research & ; 


Lawrence Livermore Laboratory 


September 1971 
166 pages, 49 illustrations, 7 x 9°/, in. 
Library of Congress Card Number: 73-170333 


Mr. Rodean investigates the generation of seismic signals by underground nuclear explosions, including 
decoupling and low-coupling phenomena that minimize these signals. Careful throughout to define terms 
and avoid ambiguities, he describes the effects of rock-property variations on seismic-signal strength and 
summarizes principles of elastic-wave theory for compressional, shear, and Rayleigh waves. Theoretical 
compressional wave solutions are given for both ground motion and seismometer response to this motion. 
Available seismic-magnitude data are presented as a function of explosion yield and environment. Mr. 
Rodean concludes with a summary of what is and is not known about seismic-signal generation by 
underground nuclear explosions. 


CONTENTS 


Inelastic Response to Nuclear Explosions 

Distribution of Energy from Underground Nuclear Explosions 
Linear Elastic (Seismic) Waves from Explosions 

Anelastic Attenuation of Seismic Signals 

Seismometer Characteristics and Operating Environment 
Distant Seismic Signals from Underground Nuclear Explosions 
Techniques to Minimize Seismic Signals 

References 

Author and Subject Indexes 


This book is available as TID-25572 for $3.00 from 
National Technical Information Service, U. S. Department of Commerce, 
Springfield, Virginia 22151 
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NUCLEAR SCIENCE ABSTRACTS 


The U. S. Atomic Energy Commission publishes Nuclear Science Abstracts (NSA), a semimonthly journal 
containing abstracts of the literature of nuclear science and engineering. 


NSA covers (1) research reports of the U. S. Atomic Energy Commission and its contractors; (2) research 
reports of government agencies, universities, and industrial research organizations on a worldwide basis; 
and (3) translations, patents, books, and articles appearing in technical and scientific journals. 


Complete indexes covering subject, author, source, and report number are included in each issue. These 
indexes are cumulated and sold separately. 


Availability 


SALE WNSA is available on subscription from the Superintendent of Documents, U.S. Government 
Printing Office, Washington, D.C. 20402, at $42.00 per year for the semimonthly abstract issues and 
$34.50 per year for the cumulated-index issues. Subscriptions are postpaid within the United States, 
Canada, Mexico, and all Central and South American countries, except Argentina, Brazil, Guyana, French 
Guiana, Surinam, and British Honduras. Subscribers in these Central and South American countries, and 
in all other countries throughout the world, should remit $52.50 per year for subscriptions to 
semimonthly abstract issues and $43.00 per year for the cumulated-index issues. The single-copy price for 
the abstract issues is $1.75 postpaid, with this exception: Add one-fourth of $1.75 for mailing to the 
countries to which the $52.50 subscription rate applies. 


EXCHANGE NSA is also available on an exchange basis to universities, research institutions, industrial 
firms, and publishers of scientific information. Inquiries should be directed to the USAEC Technical 
Information Center, P. O. Box 62, Oak Ridge, Tennessee 37830. 
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